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Preface 


Beginning with Volume XX, the Deep Space Network Progress Report changed 
from the Technical Report 32- series to the Progress Report 42- series. The volume 
number continues the sequence of the preceding issues. Thus, Progress Report 
42-20 is the twentieth volume of the Deep Space Network series, and is an uninter- 
rupted follow-on to Technical Report 32-1526, Voliune XIX. 

This report presents DSN progress in flight project support, tracking and data 
acquisition (TDA) research and technology, network engineering, hardware and 
software implementation, and operations. Each issue presents material in some, 
but not all, of the following categories in the order indicated. 

Description of the DSN 

Mission Support 

Ongoing Planetary/Interplanetary FUght Projects 
Advanced Flight Projects 

Radio Science 

Special Projects 

Supporting Research and Technology 
Tracking and Groimd-Based Navigation 
Communications — Spacecraft/Ground 
Station Control and Operations Technology 
Network Control and Data Processing 

Network and Facihty Engineering and Implementation 
Network 

Network Operations Control Center 
Ground Communications 
Deep Space Stations 

Operations 
Network Operations 
Network Operations Control Center 
Groimd Communications 
Deep Space Stations 

Program Planning 
TDA Planning 
Quality Assurance 

In each issue, the part entitled “Description of the DSN” describes the functions 
and facilities of the DSN and may report the current configuration of one of the 
five DSN systems (Tracking, Telemetry, Command, Monitor & Control, and Test 
& Training). 

The work described in this report series is either performed or managed by the 
Tracking and Data Acquisition organization of JPL for NASA. 
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DSN Functions and Facilities 

N. A. Renzetti 

Office of Tracking and Data Acquisition 


The objectives, functions, and organization of the Deep Space Network are 
summarized. Deep space station, ground communication, and network operations 
control capabilities are described. 


The Deep Space Network (DSN), established by the 
National Aeronautics and Space Administration (NASA) 
Office of Tracking and Data Acquisition (OTDA) under 
the system management and technical direction of the 
Jet Propulsion Laboratory (JPL), is designed for two-way 
communications with unmanned spacecraft traveling ap- 
proximately 16,000 km (10,000 mi) from Earth to the 
farthest planets of our solar system. It has provided track- 
ing and data acquisition support for the following NASA 
deep space exploration projects, for which JPL has been 
responsible for the project management, development of 
the spacecraft, and conduct of mission operations; 

(1) Ranger. 

(2) Surveyor. 

(3) Mariner Venus 1962. 
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(4) Mariner Mars 1964. 

(5) Mariner Venus 1967. 

(6) Mariner Mars 1969. 

(7) Mariner Mars 1971. 

(8) Mariner Venus/Mercury 1973. 

The DSN has also provided tracking and data acquisi- 
tion support for the following projects: 

(1) Lunar Orbiter, for which the Langley Research 
Center carried out the project management, space- 
craft development, and mission operations func- 
tions. 
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(2) Pioneer, for which the Ames Research Center car- 
ried out the project management, spacecraft devel- 
opment, and mission operations functions. 

(3) Apollo, for which the Lyndon B. Johnson Space 
Center was the project center and the Deep Space 
Network supplemented the Spaceflight Tracking 
and Data Network (STDN), which is managed by 
the Goddard Space Flight Center (GSFC). 

(4) Helios, a joint United States/^^'est Germany project. 

(5) Viking, for which the Langley Research Center pro- 
vides the project management and Lander space- 
craft, and conducts mission operations, and for 
which JPL provides the Orbiter spacecraft. 

The Deep Space Network is one of two NASA net- 
works. The other, the Spaceflight Tracking and Data 
Network, is under the system management and technical 
direction of the Goddard Space Flight Center. Its function 
is to support manned and unmanned Earth-orbiting and 
lunar scientific and advanced technology satellites. Al- 
though the DSN was concerned with unmanned lunar 
spacecraft in its early years, its primary objective now and 
into the future is to continue its support of planetary and 
interplanetary flight projects. 

A development objective has been to keep the network 
capability at the state of the art of telecommunications 
and data handling and to support as many flight projects 
as possible with a minimum of mission-dependent hard- 
«'are and software. The DSN provides direct support to 
each flight project through that project’s tracking and 
data systems. This management element is responsible for 
the design and operation of the hardware and software in 
the DSN which are required for the conduct of flight 
operations. 

As of July 1972, NASA undertook a change in the inter- 
face between the network and the flight projects. Since 
January 1, 1964, the network, in addition to consisting of 
the Deep Space Stations and the Ground Communications 
Facility, had also included the Mission Control and Com- 
puting Facility and had provided the equipment in the 
mission support areas for the conduct of mission opera- 
tions. The latter facilities were housed in a building at 
JPL known as the Space Flight Operations Facility 
(SFOF). The interface change was to accommodate a 
hardware interface between the network operations con- 
trol functions and the mission control and computing 
functions. This resulted in the flight project’s picking up 
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the cognizance of the large general-purpose digital com- 
puters, which were used for network processing as well as 
mission data processing. It also assumed cognizance of 
all of the equipment in the flight operations facility for 
display and communications necessary for the conduct of 
mission operations. ’The network has already undertaken 
the development of hardware and computer software 
necessary to do its network operations control and monitor 
functions in separate computers. ’This activity became 
known as the Network Control System implementation. 
A characteristic of the new interface is that the network 
provides direct data flow to and from the stations via 
appropriate ground communications equipment to Mission 
Operations Centers, wherever they may be; namely, 
metric data, science and engineering telemetry, and such 
network monitor data as are useful to the flight project. 
It accepts command data from the flight project directly 
into the ground communications equipment for trans- 
mission to the station and thence to the spacecraft in a 
standardized format. 

In carrying out its functions, the network activities can 
be divided into two general areas. The first includes those 
functions which are associated with the in-flight support 
and in tracking the spacecraft; its configuration can be 
characterized as follows: 

(1) DSN Tracking System. Generates radio metric data; 
i.e., angles, one- and two-way doppler and range, 
and transmits raw data to mission control. 

(2) DSN Telemetry System. Receives, decodes, records, 
and retransmits engineering and scientific data 
generated in the spacecraft to Mission Control. 

(3) DSN Command System. Accepts coded signals from 
Mission Control via the Ground Communications 
Facility (GCF) and transmits them to the space- 
craft in order to initiate spacecraft functions in 
flight. 

The second category of activity supports testing, train- 
ing, and network operations control functions and is con- 
figured as follows: 

(1) DSN Monitor and Control System. Instruments, 
transmits, records, and displays those parameters of 
the DSN necessary to verify configuration and 
validate the network. Provides operational direction 
and configuration control of the network and 
primary interface with flight project mission control 
personnel. 
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(2) DSN Test and Training System. Generates and con- 
trols simulated data to support development, test, 
training, and fault isolation within the DSN. Partici- 
pates in mission simulation with flight projects. 

The capabilities needed to carry out the above func- 
tions have evolved in two technical areas: 

(1) The Deep Space Stations that are distributed 
around Earth and which, prior to 1964, formed 
part of the Deep Space Instmmentation Facility. 
The technology involved in equipping these sta- 
tions is strongly related to the state of the art of 
telecommunications and flight/ground design con- 
siderations and is almost completely multimission in 
character. Table 1 gives a description of the Deep 
Space Stations and the Deep Space Communica- 
tions Complexes (DSCCs) they comprise. 

(2) Ground communications. This technology supports 
the Earth-based point-to-point voice and data com- 
munications from the stations to the Network 
Operations Control Area at JPL, Pasadena, and to 
the Mission Operations Centers, wherever they may 
be. It is based largely on the capabilities of the 
common carriers throughout the world which are 
engineered into an integrated system by the 
Goddard Space Flight Center for support of all 
NASA programs. The term “Ground Communica- 
tions Facility” is used for the sets of hardware and 
software needed to carry out the functions. 

The Network Operations Control Center is the func- 
tional entity for centralized operational control of the 
network and interfaces witli the users. It has two 
separable functional elements; namely, Network Opera- 
tions Control and Network Data Processing. 


The functions of the Network Operations Control Cen- 
ter are: 

(1) Control and coordination of network support to 
meet commitments to network users. 

(2) Utilization of the network data processing com- 
puting capability to generate all standards and 
limits required for network operations. 

(3) Utilization of network data processing computing 
capability to analyze and validate the performance 
of all network systems. 

The personnel who carry out the above functions are on 
the first floor of Building 230, wherein mission operations 
functicn.s are carried out by certain flight projects. Net- 
work personnel are directed by an Operations Control 
Chief. The functions of the Network Data Processing are: 

(1) Processing of data used by Network Operations 
Control for the control and analysis of the network. 

(2) Display in Network Operations Control Area of 
data processed in Network Data Processing Area. 

(3) Interface with communications circuits for input to 
and output from Network Data Processing Area. 

(4) Data logging and production of the intermediate 
data records. 

The personnel who carry out these functions are lo- 
cated in Building 202, which is approximately 200 m from 
Building 230. The equipment consists of minicomputers 
for real-time data system monitoring, two XDS Sigma 5’s, 
display, magnetic tape recorders, and appropriate inter- 
face equipment with the ground data communications. 
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Table 1. Tracking and data acquisition stations of the DSN 


DSCC 

LiK-atinn 

DSS 

DSS serial 
designation 

Antenna 

Year of initial 
operation 

Diameter, 

m(ft) 

Type of 
mounting 

Cnldstnne 


Pioneer 

11 

26(85) 

Polar 

1958 



Echo 

12 

26(85) 

Polar 

1962 



( Wnus )•■> 

13 

26(85) 

Az-EI 

1962 



.Mars 

14 

64(210) 

Az-El 

1966 

Ticihinliilla 

Australia 

Weemala 

42 

26(S5) 

Polar 

1965 



Balliuia 

43 

64(210) 

Az-EI 

1973 

- 

Australia 

1 loneysuekle Creek 44 

26(85) 

X-Y 

1973 

Madrid 

Spain 

Uobledo 

61 

26(85) 

Polar 

1965 



Cebreros 

62 

26(85) 

Polar 

1967 



Uobledo 

63 

64(210) 

Az-EI 

1973 

"A maintenance- 

facility. Besides the 26-m (85-ft) diam Az-EI mounted 

antenna, DSS 13 

I has a 9-m (30-ft) diam Az-EI 

mounted antenna that is used for interstation time correlation using lunar reflection techniques, for testing the design of new 

equipment, and for support of ground-bitsed radio science. 
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Compatibility Test System for Use With 
the Mark III DSN Data Subsystems 
Implementation 

A. I. Bryan 

DSN Systems Engineering Office 

G. H. Winn 

DSN Facility Operations Office 


This article reports on the Compatibility Test System (CTS) that will be used 
at the Compatibility Test Area, JPL, Pasadena (CTA 21), and the Spacecraft 
Compatibility/ Monitor Station, Merritt Island, Florida (STDN MIL 71), to 
perform telecommunications compatibility tests with the Mariner lupiter-Satum 
1977 and the Pioneer Venus 1978 spacecraft. The functional design of the new 
system utilizes capabilities provided by a Compatibility Test Assembly and the 
Mark 111 DSN Data Subsystems (MDS) configuration. A discussion of the 
Compatibility Test System implementation identifies the special-purpose 
equipment which comprises the Compatibility Test Assembly. 


I. Introduction 

Implementation of the Mark III Deep Space Network 
Data Subsystems (MDS) at the Compatibility Test Area, 
JPL, Pasadena (CTA 21), and the Spacecraft Coropatibil' 
ity/ Monitor Station, Merritt Island, Florida (STDN MIL 
71), necessitates a redesign of the Compatibility Test 
System (CTS) used to verify telecommunications compati- 
bility between a spacecraft and the Deep Space Network 
(DSN). The new CTS, as designed by the DSN Systems 
Engineering Office, will be developed and implemented at 
CTA 21 to support the Mariner Jupiter-Satum 1977 


Spacecraft-DSN Compatibility Test Program, whidi starts 
during the second quarter of calendar year 1976. The test 
system will provide capabilities to satisfy both Mariner 
Jupiter-Satum 1977 and Pioneer Venus 1978 Compatibil- 
ity Test Program requirements. 


II. Functional Design 

The functional design of the Compatibility Test System 
is based on the following criteria: 
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(1) Implement a test capability for use at CTA 21 and 
MIL 71 to provide test monitor, control and 
simulation capabilities which will satisfy Mariner 
Jupiter-Satum 1977 and Pioneer Venus 1978 Com- 
patibility Test Program requirements. 

(2) Provide a cost-effective implementation of the 
Compatibility Test System through utilization of the 
capabilities provided by the Mark III DSN Data 
Subsystems for test control and performance moni- 
toring. 

The Compatibility Test System hardware for Mariner 
Jupiter-Satum and Pioneer Venus 1978 will consist of a 
Compatibility Test Assembly interfaced to a standard 
Mark III DSN Data Subsystems configuration at CTA 21 
or MIL 71. The Compatibility Test Assembly contains 
special-purpose equipment that provides the capability to 
mea^e and analyze analog signals. Other compatibility 
test requirements will be met by utilization of the 
standard MDS capabilities for control, monitoring and 
simulation. 

The functions of measurement, processing and reporting 
of analog signals will be performed by specialized 
applications software. Elach applications program will 
perform as a high-priority task under the control of 
Operating System software. The standard Modular 
Computer Systems Corporation (MODCOMP) Operating 
System (MAX III) software will be extended to provide an 
analog input capability and a plotting capability to 
support the applications programs. 

The Compatibility Test Assembly Controller will 
process analog information to perform two functions: 

(1) Measure the phase stability of a spacecraft’s 
transmitted carrier and telemetry subcarrier(s). 
Phase stability measurements are performed in real 
time by concurrently sampling the phase detector 
error signals of two ground receivers which are 
phase-locked with the spacecraft’s downlink carrier. 
’The root-mean-square (rms) voltage from either 
receiver indicates the phase noise of the receiver 
plus that of the spacecraft’s transmitter. Cross- 
correlation of these two signals, displayed in the 
frequency and time domains, provides a measure- 

^ ment of the phase noise contributed by the 
spacecraft’s transmitter. Subcarrier phase stability 
measurements utilize the same technique by sam- 
pling two subcarrier demodulators. 

(2) Measure the relative power of the components of a 
spacecraft’s transmitted radio-frequency signal and 
provide a frequency domain plot of the spectral 


components. Measurements of a spacecraft’s trans- 
mitted spectrum components are performed by 
digitizing an analog recording of a spacecraft’s 
transmitted RF sign^. Time samples of the radio- 
frequency signal are transformed to the frequency 
domain by a Fast Fourier Transform algorithm. The 
resulting phase components of the frequency 
spectrum are processed to provide discrete spectral 
components of power. All power components are 
converted to relative values and referenced to the 
carrier. ’The resulting plot contains frequency and 
relative power values for each discrete spectral 
component. 

’The Compatibility Test System, by utilizing capabilities 
provided by the Mark III DSN Data Subsystems imple- 
mentation, will also provide the following: 

(1) Control and monitor of tests to determine the 
performance of a spacecraft’s Command, Telemetry 
and Radio Metric Subsystems under simulated 
nominal and worst-case flight conditions. Test 
control and monitoring will be performed from a 
central input/output terminal. 

(2) Spacecraft Radio Frequency Subsystem (RFS) tests 
utilize the control capabilities of the Block FV 
Receiver-Exciter and the special CTA 21 /MIL 71 
Microwave Subsystems. ’These subsystems provide 
capabilities to simulate doppler and to measure the 
acquisition characteristics and tracking performance 
of the spacecraft receiver. 

III. CTS Implementation 

Implementation of the Compatibility Test System 
requires an extension of the Mark III DSN Data 
Subsystems capabilities. ’This extension is provided by the 
Compatibility Test Assembly interfaced to the Communi- 
cations Monitor Formatter (CMF) Assembly and special- 
purpose test software, which resides on the CMF 
Assembly Disk. ’The CMF redundant MODCOMP Model 
D-25 minicomputer is utilized as the Compatibility Test 
Assembly Controller. 

Special-purpose equipment required to implement the 
Compatibility Test Assembly will be housed in one 
standard DSN 48-centimeter (19-inch) roll-aroimd rack and 
will interface with the redundant CMF minicomputer 
through a single cable coimected to an existing plug on 
the CMF Peripheral Controller Interface (PCI) Connector 
Panel. Only one set of Compatibility Test Assem'Diy 
equipment will be procured for use at both CTA 21 and 
MIL 71. Upon completion of testing at JPL and prior to 
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initiation of testing in Florida, the Compatibility Test 
Assembly equipment will be shipped from CTA 21 to 
MIL 71. 

It should be noted that when the Compatibility Test 
Assembly is interfaced to the CMF redundant minicom- 
puter, it will in no way disturb the hardware integrity of 
the standard CMF Assembly. The Compatibility Test 
Assembly interface merely extends the input /output bus 
capabilities of the standard CMF Assembly configuration. 

The new Compatibility Test Assembly consists of the 
following equipment: 

(1) A twelve-bit Analog-to-Digital (A-D) Converter and 
Power Supply. 

(2) A twelve-bit Digital-to-Analog (D-A) Converter and 
Power Supply. 

(3) A Peripheral Controller Interface Assembly, which 
contains the A-D and D-A Interface Logic. 

Standard MDS configuration equipment that will be 
used in conjunction with the Compatibility Test Assembly 
to provide the CTS capability consists of: 

(1) The CMF redundant minicomputer. 

(2) A General Electric TenniNet keyboard-printer as 
normally connected in the standard CMF Assembly 
configuration. 


(3) A Varian high-speed printer-plotter as normally 
connected in the standard CMF Assembly configu- 
ration. 

(4) The direct access storage device (Disk) that is part 
of the standard CMF Assembly configuration. 

The block diagram in Figure 1 illustrates the hardware 
configuration of the Compatibility Test System functional 
capability. 

The software packages to be used by the CTS will 
consist of: 

(1) The standard DSN Operational Software Package 
used to process Telemetry, Command, Radio Metric 
and Monitor data. 

(2) A software package consisting of a version of the 
standard MDS Operating System (MAX III) ex- 
tended to handle the interface between the Compat- 
ibility Test Assembly and the CMF’s MODCOMP 
minicomputer, and contaiiting applications programs 
such as Phase Jitter Analysis and Downlink Spec- 
trum Analysis. Applications programs will be 
prepared in Fortran IV (a standard software package 
on the MODCOMP U). 
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Fig. 1. Compatibility Test System configuration for CTA 21 and Mil 71 
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Viking Mission Support 

D. J. Mudgway 

OSN Systems Engineering Office 

D. W. Johnston 
OSN Operations Office 


This article reports DSN activity in support of Viking 1 and 2 cruise 
operations, together with the implementation and new training directed toward 
completing the planetary configuration of the DSN for Viking support by 1 
February 1976., 


I. Background 

As Viking 1 and 2 settle into the cruise phase of the 
mission, DSN Operations continues to support the Viking 
Flight Teams in providing tracking support from both 26- 
meter subnets, as well as the 64-meter subnet. At the same 
time, additions are being made to the Network, particu- 
larly at the 64-meter stations, to bring the cruise 
configuration up to meet the planetary requirements. 

This situation leads to a demand for station time which 
far exceeds the resources available. The situation is further 
aggravated by additional requirements for operational 
tests to train the station crews and to qualify the added 
capabilities for flight support. 

Some degree of resolution has resulted from the 
establishment of an integrated Viking scheduling commit- 
tee, which coordinates all DSN demands for implementa- 
tion testing and training with Project demands for flight 


support, sets priorities within this total body of work, and 
presents a total Viking support request to the scheduling 
ofiice. 

il. Planetary Testing and Training 

A considerable number of Operational Verification 
Tests (OVTs) were carried out during this period. These 
are discussed below. 

A. Operational Verification Tests 

In the previous article, the DSN testing plan for 
planetary operations was outlined. At. this time, the first 
testing phase is progressing smoothly and should be 
completed as forecast (reference Progress Report 42-30). 
The first phase testing at DSSs 11 and 14 (Goldstone, 
California) has been completed, and the program for DSSs 
42 and 43 (Australia) and DSSs 61 and 63 (Spain) is 
proceeding as planned. Listed below are the tests 
completed at each facility and a brief summary of results. 
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DSS 14 64-meter planetary OVTs. 

(1) OVT 1—90% successful, discovered more time 
required for pre-test calibrations. 

(2) OVT 2—95% successful, low SNRs on one of six 
telemetry streams. 

(3) OVT 3—90% successful, block decoder assembly 
down. 

(4) OVT 4—100% successful. 

DSS 11 and 14 combined planetary OVTs 

(1) OVT 1—50% successful, problem with microwave 
link between DSSs 11 and 14. 

(2) OVT 2—50% successful, delay in reconfigurations 
and validation of telemetry data rates. 

(3) OVT 3—100% successful 

(4) OVT 4—95% successful, SSA overheating problem. 

DSS 43 64-meter planetary OVTs 

(1) OVT 1—95% successful, minor operational problems. 

(2) OVT 2—90% successful, no DODR recalls 
performed. 

(3) OVT 3—65% successful, interference from real-time 
operations. 

DSS 42 and 43 conjoint planetary OVTs 

(1) OVT 1—100% successful. 

(2) OVT 2—100% successful. 

(3) OVT 3—100% successful. 

DSS 63 64-meter planetary OVTs 

(1) OVT 1—100% successful. 

(2) OVT 2—70% successful, interface problem between 
Ground Communication Facility (CCF) and Net- 
work Operations Control Center (NOCC). 

(3) OVT 3—90% successful, delay in data validation. 

DSS 61 and 63 conjoint planetary OVTs 

(1) OVT 1—80% successful, test start delayed, and 
equipment problems in the Network Analysis Team 
area. 

(2) OVT 2—100% successful. 

The above takes the DSN testing eflbrt up to December 

8. Phase one will be completed by early January. Phase 


two will begin the week of March 8, 1976. There wiU be 
four OVTs scheduled at each complex, consisting of two 
tests with the 64-meter stations by themselves, and two 
tests involving their 26-meter conjoint stations. This final 
OVT phase is scheduled to be completed during the week 
of April 8, 1976. 

ill. DSN Support of Additional Viking Testing 

In addition to the above DSN tests, the support of 
System Integration Tests (SITs) and Ground Data System 
(CDS) Tests is progressing as forecast in the previous 
article. To date, DSSs 11 and 14 have participated in the 
following tests: 

(1) Planetary SIT— 100% successful. Interface problem 
between Ground Communication Facility and Vi- 
king Mission Computer Control Complex was 
resolved. 

(2) CDS 5.1— Considered successful. No DSN problems 
encountered. 

Following CDS test 5.1, the next CDS test will be 5.31 
to be run December 15. This will involve DSSs 11 and 14 
at Coldstone, California. The test results will be reported 
in the next article of this series. 

In January, SIT and CDS testing will begin overseas at 
Australia and Spain. On January 8, a SIT test is scheduled 
for DSS 43 (Australia), with a retest, if necessary, on 
January 13. On January 15, a SIT is ^heduled for DSS 63 
(Spain), with a retest on January 22. 

CDS 5.32 will be run in January, with the' overse^lS 
portion beginning at DSSs 42 and 43 on January 15. On 
January 16, CDS 5.32 will be performed with DSSs 61 and 
63. In February, CDS test 6.0 is scheduled. This will 
involve DSS 14 at Coldstone and DSS 43 in Australia. The 
remaining CDS test, 11.0, is scheduled for May, but no 
station has been designated at this time. 

IV. Flight Operations Personnel Test and 
Training 

As reported in the last article, the Flight Operations 
Personnel Test and Training follows CDS testing. These 
tests are scheduled to begin in mid-February. They should 
be completed by April. Details concerning these tests are 
not readily available at this time, and will be covered in 
future articles. 
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Pioneer Venus 1978 Mission Support 

R. B. Miller 

DSN Systems Engineering Office 


The differential long baseline interferometry experiment for the purpose of 
measuring the wind velocities in the atmosphere of Venus as a part of the Pioneer 
Venus 1978 multiprobe mission is described. 


I. Introduction 

The Pioneer Venus multiprobe mission includes a dif- 
ferential long baseline interferometry experiment, which 
will attempt to measure the wind velocity in the atmo- 
sphere of Venus as four probes descend through the atmo- 
sphere. Basically, the experiment will be using inter- 
ferometry techniques to measure the components of the 
wind velocity perpendicular to the line of sight (Earth- 
spacecraft direction) and established doppler techniques 
to measure the velocity components along the line of sight. 
As described in previous articles, the bus spacecraft is 
retarded by a trajectory correction after it releases the 
four probes, so that it will enter the Venusian atmosphere 
after all four probes have reached the surface of the 
planet. In this way, the bus serves as a reference signal, 
undisturbed by the Venusian atmosphere. A corrected 
difference is taken between the bus signal and each of 
the probe signals at a particular tracking station to elimi- 
nate ionospheric and interplanetary effects, and a second 


difference is taken between pairs of tracking stations 
which produces a measure of the rate of change of the 
angle subtended by the two stations and a probe. 

Each pair of stations resolves only one component of 
the velocity. In order to resolve both components of the 
wind velocity perpendicular to the line of sight and to 
provide some measure of redundancy, four stations will be 
equipped to support this experiment. The two 64-meter 
DSN stations located at Coldstone, California, and Can- 
berra, Australia, will be utilized, along with 12-meter 
Spaceflight Tracking and Data Network (STDN) stations 
located at Santiago, Chde, and Guam. The Principal In- 
vestigator for this experiment is Dr. Counselman and the 
Co-Investigator is Dr. Pettengill, both of the Massachu- 
setts Institute of Technology. The Tracking and Data 
Acquisition Office of the Jet Propulsion Laboratory is 
responsible for seeing that all four ground stations are 
equipped for this experiment. The experimenters will be 
responsible for all nonreal-time processing of the data. 
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II. Basis of Differential Long Baseline 
Interferometry 

TTie fundamental concept of interferometry and its ap- 
plication in the Pioneer Venus case are illustrated in 
Fig. 1. Pictured are two tracking stations located on the 
surface of Earth and two signal sources located near the 
planet Venus. For simplicity, consider that the two space- 
craft and the two tracking stations are located on lines 
which are perpendicular to the Earth-Venus line of sight. 
Looking first at a single spacecraft, since the distance 
from Earth to Venus, r, is very much greater than the 
displacement of the spacecraft from the line-of-sight, y, 
then using similar triangles, the angles </• and 0 are ap- 
proximately equal. The difference in path length from the 
spacecraft to the two stations is shown as the distance 8; 
8 equals <f>d, where d is the distance separating the two 
stations. Similarly, tj is equal to rd. Since ^ is approxi- 
mately equal to 6, then y is approximately equal to (r/d)S. 
8 represents the phase difference between the single 
spacecraft signal received simultaneously at the two track- 
ing stations. As an indication of the potential power of the 
interferometry technique, if it were possible to measure 8 
to within 1 degree of phase at S-band, using the fact that 
1 Hz at S-band is approximately 13 cm, an Earth-Venus 
distance of 50,000,000 km and a station separation of 
8,000 km in the expression derived above, the displace- 
ment of the spacecraft could be resolved to within 3 
meters at Venus. Unfortunately, there are several sources 
of error which would prevent making such a direct mea- 
surement. The two most significant effects are that the 
signal will have traveled through two completely different 
locations in Earth’s ionosphere and, second, that in order 
to process the received signal, it is necessary to beat the 
signal against a local oscillator at each of the two stations. 
Differences in the local oscillators at the two stations 
would map directly into an error in the determination 
of 8. The differential technique is utilized to virtually 
eliminate both of these error sources. 

Retuniing to Fig. 1, if a second spacecraft located in 
the vicinity of the spacecraft of interest is tracked simul- 
taneously, then the same expression as derived above 
could be used to derive a differential expression: 

y, - 1/2 s -3 (81 - 82) 
ft 

Now, with a differential measurement, two important 
things happen. First, because the two spacecraft are lo- 
cated close to each other compared to the Venus-Earth 
distance, their signals will follow essentially the same ray 


path through Earth’s ionosphere; therefore, differencing 
the signals will cancel out the ionospheric and atmospheric 
effects of Earth. Second, if the two signals are received 
through a single receiver at each of the stations, then at a 
particular station, both signals will have been beat against 
the same oscillator, and therefore, when they are differ- 
enced, the variations in the local oscillator will cancel out. 
This latter point is very important when considering the re- 
quirements on the ground equipment. In looking at the 
error sources of the experiment due to contributions from 
station equipment, because of this differential effect, only 
error sources which introduce a differential phase error 
between the two received signals have a significant effect 
on the experiment. Error sources which cause equal 
changes in the two received signals (such as local oscilla- 
tor drift) have no first-order effect on the experiment. 

In the above simplified discussion, one point which was 
ignored is the fact that 8 in practice contains an unre- 
solvable ambiguity. In practice, 8 is many wavelengths 
long (one S-band wavelength is approximately 13 cm, 
where 8 can be on the order of hundreds of kilometers) 
and should be better represented by the expression 


where X is the signal wavelength, n is an unknown integer, 
and p is the fractional phase difference expressed in 
radians, p is defined as the fringe phase. Since it is not 
possible to determine n to sufficient accuracy, only the 
time variation of p is meaningful, and it is the time varia- 
tion of p which is termed the fringe frequency in radio 
interferometry. However, in the Pioneer Venus case, that 
is exactly what is desired. Determining the time rate of 
change of p, we have the time rate of change of 8 and 
therefore the derivative of y, which represents the velocity 
of the spacecraft perpendicular to the line-of-sight; 

(yi - y?) ^ _r d (Si - 82) _ r, dp 
dt ~ d dt d lit dt 

where dp/dt, the fringe frequency, is the observable. 


III. Differential Long Baseline interferometry 
Requirements 

The key requirements for the Pioneer Venus 1978 dif- 
ferential long baseline interferometry wind measurement 
experiment will be briefly described. 
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The objective of the experiment is to btable to measure 
the wind velocities as the probes fall through the atmo- 
sphere of Venus to about a 10-cm-per-second accuracy 
using a 100-second integration time for the signal-to-noise 
ratios expected at the DSN 64-meter antennas. It is fairly 
easy to show that this requirement translates into the need 
to be able to determine the phase difference between pairs 
of signals when averaged over 100 seconds to within 1 
degree of phase at S-band. For the 9-meter STDN sta- 
tions, which will have a significantly less favorable signal- 
to-noise ratio, it will clearly be necessary for the experi- 
menters to integrate over much longer times in order to 
achieve the same accmacy, therefore sacrificing time reso- 
lution in the rate of change of the velocity. It is this 
1 degree of relative phase error versus time over the 
bandwidth of interest requirement that will be the most 
difficult to meet, and the DSN and the STDN have not 
yet determined what can actually be achieved. Pre- and 
post-experiment calibration of the station equipment 
involved in the experiment will be necessary, as well as 
some form of calibration signals recorded along with the 
actual data. The experimenters are confident that a num- 
ber of that order can be achieved based on similar experi- 
ments that were performed at lunar distances using 
ALSEP signals. 

As was described in previous Progress Report articles, 
the total bandwidth which the five signals from the multi- 
probe mission might occupy (four probes plus the bus) is 
1.7 MHz. It is therefore necessary to have receivers which 
can pass a 2-MHz bandwidth, and open-loop receivers 
will be modified for this purpose. Analog recording is felt 
to be incapable of meeting the differential phase require- 
ment of this experiment, and therefore digital recordings 
will be made in real-time. Three-bit quantization is re- 
quired, and this, together with the 2-MHz bandwidth, 
means that the recorders will have to be able to operate 
at at least 12 megabits per second. These recorders repre- 
sent the most significant implementation for the Pioneer 
Venus mission. 

The experiment also requires that the mean rate frac- 
tional accuracy of the sampling be three parts in 10‘^ and 
that the jitter on the samples be held to 10 nanoseconds 


root-mean-square. Additionally, the calibration tones 
which will be inserted' in real-time should have an abso- 
lute accuracy of three parts in 10’*. This requirement wfll 
be met by hydrogen masers set with cesium standards at 
the DSN stations and cesium standards at the STDN 
stations. 


IV. Ground Station Configuration 

Four stations will be equipped for the Pioneer Venus 
1978 differential long baseline interferometry wind mea- 
surement experiment. Figure 2 is a block diagram of the 
configuration which will be implemented at the DSN 
64-meter Coldstone and Canberra stations and die STDN 
9-meter Santiago and Guam stations. The five spacecraft 
signals will be detected by low-noise amplifiers, which 
at the STDN stations will be parametric amplifiers with a 
total system temperature of 100*, and at the DSN stations 
ruby masers with a total system temperature of less than 
24*. Some form of yet-to-be-determined calibration ref- 
erence frequencies will be inserted at this point in order 
to calibrate out drifts in the system. The five signals plus 
calibration tones will then pass throu^ the open-loop 
receivers, which will pass a 2-MHz bandwidth. Signals 
will then go through an analog-to-digital converter and 
sampler and onto the 12-megabit-per-second digital re- 
corders. The recorders will be redundant at each of the 
stations. The recorders will be able to record the 12- 
megabit-per-second rate at 76 cm (30 inches) per second, 
and therefore 80 minutes of recording will be possible 
between tape changes. 

There are two principal remaining open areas in the 
differential long baseline interferometry wind measure- 
ment experiment. First is the determination by the DSN 
and the STDN of what is the actual differential phase 
error achievable by the ground equipment at a given 
signal-to-noise ratio. The second area concerns the details 
of both the pre- and post-experiment and real-time cali- 
bration. The complexity and sophistication of the required 
calibration will be dictated by the DSN and the STDN 
determined error contributions introduced by each of the 
elements in the ground station configuration. 
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Fig. 2. Differential long baseline interferometry configuration for Pioneer Venus 1978 multiprobe 

wind measurement 
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Helios Mission Support 

P. S. Goodwin 

DSN Systems Engineering Office 

W. G. Meeks, R. E. Morris, and S. E. Reed 

Network Operations Office 


Having observed its first year anniversary in orbit about the Sun, the Helios-1 
spacecraft has performed satisfactorily. During its first year, valuable scientific 
data were obtained about the Sun and our solar system. As Helios 1 enters its 
second year of successful operation, preparations for the launch of Helios B are 
in their final stage. 


I. Introduction 

This is the seventh article in a series that discusses the 
Helios-1 Bight support and Helios-B test and training 
activities by the Deep Space Network (DSN). The 
previous article (Ref. 1) reported the results of the Helios- 
1 second superior conjunction and second perihelion 
operations as well as Helios-B test results. T^ article 
covers Helios-1 cruise operations after 'the second 
perihelion through second aphelion. Also discussed are 
significant spacecraft anomalies, DSN tracking coverage, 
Helios DSN-STDN cross-support, and Helios-B test and 
training results. Helios-B launch is presently scheduled for 
0534Z on January 15, 1976. 

II. Mission Status and Operations 

A. Helios-l Second Aphelion Operations 

One year of successful fii^t operations was completed 
by Helios 1 on 10 December 1975. The solar probe with 
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its attendant scientific instruments has performed well, 
with the only exceptions being the . initial problems 
encountered with Experiments 1 and 5 after launch (Ref. 
2), and the recent Traveling-Wave Tube Amplifier-1 
(TWTA-1) failure on 31 October 1974. 

The failure of TWTA-1 manifested itself when an 
approximate 28-dB drop in the downlink signal level was 
observed during a spacecraft tracking pass over DSS 11 at 
Goldstone. The anomaly occurred at 1834 GMT, and to 
accelerate isolation, Goldstone DSS 12 support was 
requested by the DSN. Both stations verified identical 
downlink signal leveb at -167 dBm. This level was below 
the telemetry threshold for 26-meter stations. The Helios 
Project declared a spacecraft emergency at 2040 GMT 
and requested 64-meter support for the impending 
Australian viewperiod. At 2121 GMT, the Project 
commanded the spacecraft from the failed higbipower 
TWTA-1 mode (20 watts) to the 10-watt medium-power 
mode. This seemingly returned the probe back to a 
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normal operational mode. What had been thought of as 
normal operations was short-lived; a total TWTA-1 failure 
disrupted the DSS 43 downlink at 2253 GMT. The Project 
analyzed this failure, and the spacecraft was commanded 
to the low-power mode (0.5 watt) of operations at 2344 
GMT. Once RF lock was achieved, the spacecraft was 
commanded to 16 bps, and normal operation in this mode 
was achieved. 

With a satisfactory low-power downlink established, the 
Project continued to analyze the anomaly. The decision to 
switch to TWTA-2 medium-power mode was accom- 
plished at 1610 GMT 1 November 1975. The spacecraft 
has operated normally since the switch to TWTA-2. 
Special DSN surveillance was provided by DSSs 11 and 63 
until the spacecraft emergency was terminated during the 
DSS 11 track of the spacecraft on 1 November 1975 at 
2130 GMT. 

A Helios Project-sponsored TWTA Failure Analysis 
Review was conducted at JPL on 12-14 November 1975. 
The outcome of these meetings produced 14 action items 
plus 7 recommendations, directed to various Helios team 
members. The DSN was requested to continue to provide 
operational information on whether automatic gain 
control fluctuations are observable within the downlink 
after the switch to TWTA-2. The Failure Review also 
illustrated that further investigations will have to be 
conducted to reach a final conclusion about this TWTA 
failure. 

B. Spaceflight Tracking and Data Network (STDN) 
Cross-Support 

A DSN-STDN Helios Interface Agreement is now in 
preliminary form and will provide the necessary DSN- 
STDN Operations plans and interfaces which are required 
for Helios-1 and Helios-B STDN cross-support. The DSN- 
STDN cross-support configuration has also been defined, 
and 15 January 1976 is the designated target date to 
commence operations on a regular basis. Initially, only the 
Madrid STDN station will be available for support; the 
Goldstone STDN station, currently down for major 
rework, will become operational and provide Helios 
support after 1 March 1976. 

There is one remaining interface problem concerning 
timing signal processing during tape recording playback. 
The problem has shown up as a time-tag ofiEset between 
the Spacecraft Compatibility/Monitor Station, Merritt 
Island, Florida (STDN (MIL 71)) playback data records, 
produced from the Goldstone STDN tapes, and the 
Mission Control and Computing Center (MCCC) data 
records which were simultaneously recorded in real-time 

16 


from DSS 12. The relative time-tag differences are a 
function of the Data Decoder Assembly (DDA) perform- 
ance and the recorded bit rates. Indications are that the 1 
pulse-per-second and 1000 pulse-per-second timing signals 
from the Time Code Translator to the DDA interface at 
STDN (MIL 71) are incorrect. A modified Time Code 
Translator has been shipped from JPL to STDN (MIL 71). 
When the timing problem is corrected, a test will be 
conducted to re-evaluate the STDN (MIL 71)-MCCC 
interface. The re-evaluation should be completed by the 
first week of January 1976. 

C. Helios-B Test and Training 

As we approach the laimch of Helios B, DSN and 
Ground Data System testing is being successfully com- 
pleted, and emphasis is shifting to Mission Operations 
System testing. Simulation System, Ground Data System, 
and DSN operational testing were completed during this 
reporting period. Only the DSN Configuration Verification 
Tests (C\^s) and the Mission Operations System tests 
remain before the final Operational Readiness Test and 
launch of the Helios-B spacecraft. 

1. Ground Data System Test. One Ground Data System 
test was conducted on 5 December. This was a retest to 
validate a few items which failed to meet specifications 
during other U.S.-German Groimd Data System tests. This 
test was successful, completing Ground Data System 
testing for Helios B. 

2. DSN Testing. Network testing for initial acquisition 
with Australia Deep Space Stations 42 and 44 was 
completed by mid-November. Each station’s operational 
crews participated in at least one Operational Verification 
Test; some in as many as three. Test scheduling conflicts 
with operational commitments, plus some doubt about the 
exact Helios-B launch date (the official date of 15 January 
1976 being aimounced later), prevented the conduct of as 
many initial acquisition Operational Verification Tests 
with the crew selected to support launch activities as 
would have been desired. However, several Mission 
Operations System tests will exercise these procedures 
with the Station launch support crew. 

Step n maneuver Operational Verification Tests began 
with the first test at Goldstone DSS 12 on 7 November. As 
of this writing, six successful tests were conducted, four 
with DSS 12 and two with Goldstone DSS 11. 

The plan for DSN support for the Helios-B Step II 
maneuver will differ slightly from that of Helios 1. This 
difference lies in the antenna polarization configuration 
prior to track. DSS 12 (prime) will configure for linear 
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horizontal polarization. DSS 11 will configure for linear 
vertical polarization. As the spacecraft’s aspect angle is 
changed by attitude commands, the downlink signal 
polarization will go from linear horizontal to linear 
vertical. Thus the signal level received at DSS 12 will 
decrease while that at DSS 11 will increase. When the 
signal strength at DSS II surpasses that at DSS 12, the 
Project will transfer the uplink signal to DSS 11 for the 
remainder of the track. DSS 12 will reconfigure to linear 
vertical polarization for the remainder of the pass. 

3. Mission Operations System Test and Training. 
Mission Operations System test and training activity 
continued at a predetermined pace throughout this 
reporting period. Three Inter-Team Training Tests, 
without DSN participation, were conducted in October. 
The DSN support began with two Combined Test and 
Training exercises in mid-November. Four Operational 
Demonstration Tests are scheduled from 25 November 
through 6 January 1976. These tests are designed to 
demonstrate the combined U.S.-German Networks ability 
to properly support all phases of the mission from launch 
to cruise. To date, all Mission Operations System testing 
has met objectives and has been considered successful. 

The Operational Readiness Test is scheduled to take 
place on 12 January 1976. Launch for Helios B is now 
scheduled for 0534 GMT on 15 January 1976. 

D. Actual Tracking Coverage Versus Scheduled 

Coverage 

HeiioS'l tracking coverage, along with Pioneer 10 and 
11 tracking coverage, is now allocated on an equal priority 
basis for extended mission operations. This report covers a 
56Klay period for Helios-1 tracking coverage, from 17 
October through 11 December 1975. Operations during 
this period consisted of normal cruise phase operations. 

The total tracking coverage shared between the two 
Pioneers and Helios-1 spacecraft was 322 passes, equaling 
2290 hours. The Helios allocation was approximately one- 
third of this total, or 99 passes and 699 hours of tracking 
coverage. 

Due to Viking Sight operations, the number of Helios-1 
passes, the number of hours tracked, and the percentage 
of coverage supported during this period were ail 
decreased from those of the last period. In conjunction 
with the other decreases, the duration of the average pass 
dropped from 7.6 hours to 7.0 hours. Inasmuch as the 26- 
meter subnet can support at a data rate of 2048 bps, its 
support increased 58 percent during this time frame; but 
the 64-meter network only supported Helios 1 for 183 


hours, which was a drop of 71 percent from the last 
reporting period. 

Helios-1 coverage should remain relatively stable at the 
level allocated during this period if the STDN support is 
provided as plaimed. 

III. DSN System Performance for Helios 

A. Command System 

The command activity for Helios 1 continued to rise 
during this reporting period. The spacecraft’s second 
perihelion occurring on 21 September 1975 (coverage 
lasting until 5 October) and the fact that this is the 
spacecraft’s closest approach to Elarth since launch phase 
have contributed to this increase in the spacecraft 
scientific data, and resulted in increased command activity. 
A total of 4669 spacecraft commands were processed 
through the Command System in the months of October 
and November. This is an increase of 1320 commands over 
the last period, and boosted the DSN cumulative 
command total to 23,301. Two Command System aborts 
were experienced in October; none in November. Both 
system aborts were attributed to operator error and 
involved no command equipment malfunctions. The 
Command System aborts cumulative total was increased 
to 5, with the overall abort total at 13 since spacecraft 
launch. 

Command System downtime due to equipment prob- 
lems during this reporting period was only 3.81 hours. 
There were an additional 1.95 hours lost due to high-speed 
data line outages. These figures show a significant decrease 
in downtime over the last period. Three transmitter 
failures were included in the 17 instances reported. Two 
of these occurred at Goldstone DSS 14 and one at 
Australia DSS 42. 

B. Tracking System 

The DSN Tracking System continued to provide very 
good support to the Helios Project during the second 
perihelion and aphelion. Several anomalies were detected 
during this period which required considerable support 
and analysis of the DSN Tracking System data. 

What was later diagnosed as a Helios-1 spacecraft 
failure was noticed on 9 October. Ranging modulation on 
the downlink was just barely detectable, even at 64-meter 
statiorts operating with high power. At the time, it was 
believed that spacecraft ranging was lost forever. How- 
ever, this was later determined to be associated with the 
spacecraft Traveling-Wave Tube Amplifier (TWTA-1) 
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failure of 31 October. Ranging returned after the switch to 
TWTA-2. To utilize Helios-1 spacecraft ranging the 
modrilator-exciter 1 must be configured to TWTA-2 by 
command froiii the Project. The first ranging in this 
maimer was used on 9 November 1975. This operational 
procedure is being used and all spacecraft ranging 
commitments are being met. 

The changeover from the special 2-MHz doppler bias 
for Helios-1 tracks to the standard 1-MHz bias occurred as 
scheduled on ll October. Almost immediately, Australia 
DSS 43 and Spain DSS 63 began to show a very large 
doppler residual. The problem was traced in real-time to 
the Biased Doppler Detector Module of the floppier 
Extractor. When the unit was tested at the Helios 1-MHz 
biased doppler frequency, a harmonic distortion showed 
up. Several replacement modules were tried. Some gave 
good doppler residuals while some did not. It was 
concluded that operation of the module so near its design 
specification is marginal and risky. 

The Helios-B DSN Initial Acquisition Study, detailing 
DSN acquisition strategy at Australia DSS 42, was 
released on 10 November 1975. The report, prepared by 
the Network Operations Analyst for Tracking, was very 
comprehensive. Portions of the study will be published as 
part of the Helios-1 and Helios-B Network Operations 
Plan (Document 613-21), Revision B. 

C. Telemetry System 

Telemetry System analysis, as an on-going efifbrt, was 
initiated during this reporting period as a result of the 
Traveling- Wave Tube Amplifier Failure Study Meeting. 
This action item was to analyze the Helios-1 ACC 
dropouts occurring after the switch to TWTA-2. The 
intent of the study which followed was to determine the 
presence of dropouts and to compare them with those 
dropouts previously found on TWTA-1. ACC strip-chart 
data were collected from several 26-meter stations which 
tracked Helios 1 from 14 to 22 November 1975. 
Examination of the data found that large cyclic variations 
in the ACC clouded any search of dropouts. Automatic 
Gain Control variations ranged from 2 dB to 5 dB in 
amplitude, with a period of approximately 20 seconds. 
This periodic variation seemed to have a secondary 
fluctuation (possibly ACC dropouts) having an amplitude 


of about 0.5 dB with an occasional 1 to 2 dB, and an 
occurrence frequency of approximately 30 per hour. 
Starting on 21 November, the periodic fimction decreased 
in amplitude, and by Pass 348 on 22 November (over DSS 
11) the data showed that the periodic variation bad gone. 
Remaining, however, was the secondary fluctuation which 
appeared to be ACC dropouts. (Figure 1 gives a picture of 
a typical strip-chart record.) The study is continuing. 

There were only six Telemetry System Discrepancy 
Reports recorded during October, with an average mean- 
time-between-failures of 62 hours. The breakdown of 
November Discrepancy Reports is not yet available. 

Helios-1 telecommunication link analysis for October 
and November is still in progress at this writing by the 
DSN Telemetry Analysis Group. 

IV. Conclusions 

Just prior to reaching its second aphelion, the Helios-1 
spacecraft experienced its first major malfunction, a 
Traveling-Wave Tube Amplifier failure. Extra DSN 
tracking coverage and surveillance were given in support 
of this spacecraft emergency. Special recording and 
analysis of telemetry data were provided in support of the 
Helios Spacecraft Team’s investigation of the anomaly. 

A DSN-STDN Helios Interface Agreement was pre- 
pared during this reporting period which defines inter- 
faces, configurations, and dates for STDN support of the 
Helios Project. Minor problems exist, but solutions are 
seen in the near future, and the plan has promise for 
providing future Helios telemetry coverage. 

Helios-B launch preparations are progressing according 
to plan. Several test phases have been successfully 
completed during this reporting period. MOS testing will 
culminate with the Helios-B Operational Readiness Test 
scheduled for 12 January 1976. First opportunity for 
launch is 0534Z on 15 Januaiy 1976, and overall testing 
success has provided reassurance that all elements of the 
DSN are prepared for this opportunity. 

With the launch of Hdios B and increased activities of 
the Viking Project in preparation for the Mars encounter 
and landing, DSN coverage of Helios 1 is likely to 
decrease in favor of Helios B. 
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Very Long Baseline Interferometry 
Covariance 

J. W. Layland 

Communications Systems Research Section 


This article summarizes the methods and results of a covariance calculation 
for determination <)f DSN station location and clock parameters using very long 
baseline interferometry (VLBI) observation of Quasars. Errors of observation are 
assumed to arise from random (thermal) noise at the receiving antenna, and 
from errors in modeling the tropospheric delay. The critical error source is seen 
to be the troposphere: If errors in troposphere modeling are uncorrelated 
between observations, then 30-50 cm station-location determinations can be 
obtained with current modeling accuracies; if, however, these errors are 
correlated, then a factor-of-ten improvement in prediction of the tropospheric 
delay is needed to achieve this accuracy range. 


I. Introduction 

This article summarizes the methods and results of a 
covariance calculation for determination of DSN station 
location and clock parameters using VLBI observation of 
Quasars. A given observation set may currently include up 
to five stations and ten radio sources. Unknown parameters 
assumed to be simultaneously solved include the three- 
dimensional station location, two station clock parameters 
and the location of the Quasars. Both time delay and 
fringe frequency are observed. Errors of observation are 
assumed to arise from random (thermal) noise at the 
receiving antenna and from errors in modeling the 
tropospheric delay. The critical error source is seen to be 
the troposphere: If errors in troposphere modeling are 
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uncorrelated between observations, then 30-50 cm station- 
location determinations can be obtained with current 
modeling accuracies; if, however, these errors are 
correlated, then a factor-of-ten improvement in prediction 
of the tropospheric delay is needed to achieve this 
accuracy range. Thermal noise, which is independent 
between observations, is comparatively insignificant 
relative to current troposphere errors if sufficient data are 
recorded for reliable detection of each Quasar (10^ to 10^ 
bits per observation). 

II. Problem Model 

The radio stars are located in space by flieir right 
ascension (RA) (a) and declination (dec) (S). Tracking 
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id/dt)(,Ap;) = -^Pz «K 7 


stations are assumed to lie on the surface of a spherical 
Earth, and are hence also located by their a and S 
coordinates and by a third coordinate, p, representing 
local variations from the mean Earth radius. (The 
geometry of the problem is described in Ref. 1.) Station 
clock offsets and clock rate offsets are assumed to be 
nominally zero but are treated as solved-for parameters in 
the covariance calculation. An observation is possible 
whenever any radio star is simultaneously in view of two 
or more of the tracking stations. A radio star is in view of 
a station whenever its elevation angle at that station is 
above some prespeciiied limit value. An observation set 
consists of a 24-hour interval during which the set of stars 
to be observed is rotated through the sky in steps of 
(typically) sr/100 radians. At each step, any observation 
which is possible is added to the data set. No attempt is 
currently made to optimize the observing pattern or to 
eliminate conflicts in station usage. Some observations 
may occur simultaneously at three or more stations if their 
locations permit. The number of stations and their 
locations, and the number of stars and their locations are 
all variable data parameters. Most of the calculations to be 
described in the following have been performed using the 
three DSN stations at Coldstone, Calif., Canberra, 
Australia, and Madrid, Spain, for which no three-station 
observations are possible. 


A VLSI observation results in noisy estimates of the 
time delay and doppler frequency difference between the 
observed Quasar noise signab at the two stations. Errors in 
the observation arise from thermal noise, from errors in 
modeling the propagation medium, and various other 
causes not to be considered here. Earth’s troposphere 
consists of a thin shell over its surface. At Zenith elevation 
angle, this shell adds approximately 2 meters of additional 
delay to the path of a radio signal. As the elevation is 
lowered, the signal path travels through an increasingly 
longer segment of this tropospheric shell. If we let Ap^ be 
the standard deviation of the error in modeling the 
tropospheric delay at zenith, and assume that this 
modeling error scales as the length of the in-troposphere 
signal path, then the standard deviation of the tropo- 
spheric delay model error at an elevation y is 

Ap, = Ap,/lA + sin (y)] (1) 

approximately, where A °= 0.02. Furthermore, if we 
assunie that any given observation is taken through a zone 
of troposphere which is reasonably layered and static, then 
the frequency error induced into the measurement results 
directly from Earth rotation, and changes in y. 


X (d/df)(y)/[A + sin yP (2) 

The frequency error inferred in this way is clearly 
correlated with the delay error. Nothing in this formula- 
tion prevents us from assuming that different observations 
are taken through different zones of troposphere, and 
hence that the tropospheric modeling errors for each 
observation are independent. Current modeling capability 
is believed to be accurate to about 3% (Ref. 2), giving us 

Ap, 2 X lO"*® sec (3) 

using only surface weather measurements at the tracking 
stations. It is expected that this figure can be reduced by a 
factor of 3-5 by using water vapor radiometer calibra- 
tions. 

The effect of thermal noise depends upon the amoimt of 
data recorded and the time-and-frequency space spaimed 
by the data-taking operation. Let us suppose that the 
spanned bandwidth of an observation is 20 MHz and the 
spanned time is 200 seconds. Suppose further that enough 
data are recorded to achieve a 6% resolution of the 
resultant ambiguity function; this can, for example, be 
achieved with 10^ recorded bits from a 0.5 flux-unit ratio 
source observed by two 64-m antennas (Ref. 3). With this 
strategy, the standard deviation of the time delay 
measurement error is 

At * 3 X 10"® sec (4) 

Likewise, the frequency measurement error is AVf ■= 3 X 
10^ Hz. Assuming an ^band center frequency, and scaling 
by the Earth rotation rate for numerical convenience, we 
define 

AD = Af/w, « 2 X lO"® sec (5) 

The three parameters Ap^ Ar, and AD are variable data 
for the covariance calculations, with nominal vidues given 
by Eqs. {3)-{5). For any observation, the standard deviation 
of the error of observation is the rms of the thermal noise 
of the observation, and the tropospheric model error at 
each of the two stations. In the covariance calculation, all 
observations may be equally weighted, or weighted by the 
reciprocal of their variance. More sophisticated weight- 
ings dependent upon the errors projected into station 
locations could presumably improve the results. 

Although the observed parameters are the Quasar noise 
signal time delay and doppler, the parameters of interest 
are the geometrical station and star locations. The time 
delay and doppler at a station, S„ relative to a fictional 
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station at the Ceo-Center, are defined in terms of the 
observing geometry by the two nonlinear equations: 

T bZ, sin fi + r, cos 6 cos (a, - a) 

+ C, + Ca, (6) 

D “ -r, cos sin (o, - a) + C, (7) 

where (a, 6) are the Quasar RA and dec, (a„ r^ Z,) the 
cylindrical coordinates of S, at the time of observation, 
and C, and C, are the station clock offset and clock rate 
ofiEset. Using R, as the mean Earth radius, Z, = R, {p, sin 
Br) is the projection of the station position on Earth’s spin 
axis, and r, « R, p, cos £, is its spin radius, or the 
projection on the equatorial plane, where (a„ p,) are 

the previously defined spherical coordinates of *1110 
observed time delay and doppler are of course the 
differences between the geocentric delays and dopplers at 
two observing stations. 

The parameters of interest are derived from the 
observed parameters through a nonlinear least-squares 
estimation procedure. Details of such procedures in a 
radar context may be found in Ref. 4. The parameters to 
be solved for are station locations, rather than baselines, 
for ease in handling multiple baseline data. One real 
station, usually Coldstone, is used as inference for station 
coordinates, and one Quasar RA is also used as reference. 
The solution is assum^ to be performed in a cylindrical 
coordinate system. The covariance of the errors resulting 
from the nonlinear least-squares estimation procedure is 
estimated by pretending that the problem is linear in the 
neighborhood of the solution point and applying the 
covariance result for a linear least-squares problem. 


III. Summary of Calculations and Results 

The first set of calculations, shown in Fig. 1, was 
intended to provide some direct comparison with the 
published results of Thomas and Fanselow (Ref. 5). These 
calculations do not include tropospheric errors. Two 
difiSculties arise with this comparison. First, the coordinate 
systems used are not identical, so direct vector comparison 
is not convenient, although rss error-length comparison is. 
Second, the shape and size of the location error pattern 
has been observed to be a very strong function of the 
observation strategy. The observation set of Ref. 5 was at 
least partially optimized, whereas the observation set used 
here was wholly naive. The rss of the Madrid location 
error is shown in Fig. 1 for two observation sets for At ■= 
1 X 10~^ second, and various AD. Other values of Ar can 
be determined by appropriate scaling of both At and AD. 


The poorer of the two observation sets shown consists of 
118 observations of three stars located at (RA, dec) ■= (0., 
0.8), (1., -0.05), (2.2, 0.1). The better observation set 
consists of 95 observations of these same three stars, plus 5 
observations of a star located at (1.3, 1.5). *rhe Madrid 
location error for both of these sets is noticeably poorer, 
with larger doppler errors than the Coldstone-Madrid rss 
baseline error as estimated from the results reported in 
Ref. 5. 1 believe at this time that this difference results 
wholly from the comparison difficulties cited above, 
principally the difference in observation strategy. 

The remainder of the calculated errors to be presented 
include tropospheric effects, and most are computed at 
the nominal parameter values of £)qs. (3)-(5). For each 
calculated experiment, only the largest vector component 
of the station location error set is graphed; in most cases, 
this largest component is the Z-projection of the 
Canberra, Australia, station. 

Figure 2 shows the one-sigma station location errors for 
a VLBI experiment involving the three DSN 64-meter 
stations and three radio stars at declinations of 0 and ±0.55 
m. The lower set of curves shows the expected resultant 
station location error for this particular experiment set. 
The upper set of curves shows the station location error 
scaled by the square root of the number of observations, 
which itself varies with elevation angle. This serves as a 
figure-of-merit for the observation pattern, and also 
indicates crudely what the deterioration in results would 
be if the errors of observation were correlated rather than 
independent. The improvement in results by wei^ting of 
the input data is also clear from this figure. 

Some very low-elevation observations result in dramatic 
reduction of the location error indicated in Fig. 2. Upon 
inspection, it appears that these critical observations are 
the very few observations possible on the improbably long 
Madrid-Canberra baseline. With the exception of these 
observations, the best elevation limit would appear to be 
in the neighborhood of 10 to 12 degrees, and perhaps such 
a limit is appropriate on baselines with better visibility, 
with observations made at 3 to 4 degrees on the Madrid- 
Canberra baseline only. 

Figure 3 is another look at the Canberra declination 
error for a VLBI experiment involving the three DSN 64- 
meter stations and three radio stars at declinations of 0 
and ±0.35 m. This smaller declination spread allows 
observations to be made on the Canberra-Madrid baseline 
at sli^tly higher elevations, which is, in turn, reflected in 
a different signature for location error vs. elevation. The 
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decreased spread in declination of the stars results in some 
decrease in the attainable station location resolution. 

Looking at Fig. 2, there is some temptation to believe 
that the observations at or near the horizon are the most 
important to reducing the station location errors. Accord- 
ingly. a covariance calculation was performed which 
limited the elevation angle of at least one station of the 
pair to a band within 10 degrees of the selected horizon. 
This is shown in Fig. 4 for &e same data as in Fig. 2. At 
best, only a slight benefit has accrued from this strategy at 
the higher horizon limit values, and a sizeable loss has 
accrued for low horizon limits. 

Figure 5 shows the Canberra Z error obtained when the 
observations are perturbed by the tropospheric modeling 
error only. As in Fig. 2, the stars are located at 0 and 
±0.55 m. This approximates the situation obtained when a 
very large amount of data is recorded for each observa- 
tion. The dependence upon the very low-elevation 
observations is enhanced in this figure. 

There is at least an academic interest in knowing what 
improvements can be achieved by using additional 
tracking stations. This has been investigated for specific 
additions, with the results shown in Figs. 6 and 7. The stars 
observed are at declinations of 0 and tO.55 m. Figure 6 
shows the Canberra Z error calculated for an experiment 
using the DSN plus the Spaceflight Tracking and Data 
Network (STDN) station in Hawaii. Figure 7 shows the 
Canberra Z error for a 4-station VLBl experiment 
including the former DSN station at Johaimesburg, South 
Africa. Roughly speaking, the addition of a properly 
placed station to the DSN results in a sizeable improve- 
ment in the largest station location error component and a 
reduced need for low-elevation observations. 

Figure 8 shows the effect of varying the thermal noise 
contribution to a VLBI experiment involving the three 
DSN stations and four radio stars at declinations of -0.5, 0, 

0.17 and 0.8 radians. The noise contribution is indicated as 
number of bits recorded per station, although such 


variation can also result from varying source strengths. 
The point corresponding to the noise parameters of Eqs. 
(4) and (5), of 10* bits per observation, is shown for 
reference. Even if we assume optimistically that tropo- 
sphere modeling errors are independent, we can operate 
with the total bit count in the neighborhood of 10* bits 
per observation, or 2 X 10® bits per total experiment— 
which fits on one computer tape-^with little loss in 
resolution. If the assumption of independent troposphere 
errors is false, the thermal noise is still more strongly 
dominated by tropospheric effects; then, however, the 
desired 50-cm accuracy becomes unachievable. The truth 
should lie somewhere between these curves. 

IV. Conclusions 

Only a few conclusions arise naturally from the results 
presented here; 

(1) Tropospheric modeling errors represent a dominant 
error source, which masks thermal receiver noise as 
long as enough data is recorded for detection. 

(2) Independence between samples for the errors in the 
troposphere model is crucial to the delivery of 
meaningful results. If it is not achieved, then at least 
an order-of-magnitude improvement in model accu- 
racy will be required. 

(3) Observations on the Madrid-Canberra baseline are 
important when only the DSN stations are utilized. 

V. Future Work 

Several specific actions seem required to complete the 
analysis presented here: 

(1) Add the capability to do directed observations as 
well as clock-driven observation sets. 

(2) Add terms for ionospheric and other error sources. 

(3) Experiment further with observations and weighting 
strategies. 

These actions will be accomplished as time permits. 
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ELEVATION LIMIT 
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Fig. 2. Canberra location error va. horizon with 3 
DSN atationa and 3 atara at 0 and ±0.55 m 
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CANBERRA Z-IOCATION ERROR, cm 
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Fig. 3. Canberra location error vs. horizon with 3 
DSN stations and 3 stars at 0 and ±0.35 m 
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Fig. 4. Canberra location error with horizon 
window for 3 DSN stations and 3 stars 
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CANBERRA Z-LOCATION ERROR 
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Fig. 5. Canberra location error vs. horizon for 3 
DSN stations, 3 stars, troposphere noise only 
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ng. 6. Canberra location error for DSN plus 
Hawaii and 3 stars 
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Fig. 7. Canberra location ernir for DSN plus 
Johannesburg, South Africa, and 3 stars 
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Binary Sequential Ranging With Sine Waves 

W. L. Martin and J. W. Layland 

Communications Systems Research Section 

Current ranging systems of the Deep Space Network estimate the range to a 
spacecraft by measuring the phase of a square-wave range code after its round- 
trip from station to spacecraft and back. Distortion of this waveform, in the form 
of phase-shifts of the harmonics of the square wave, can seriously degrade ranging 
accuracy. In this report, we show that such degradation can be largely eliminated 
by zonal filtering the range code to only its sine-wave fundamental at no net cost 
in ranging signal power or accuracy. 


I. Introduction 

Current ranging systems of the Deep Space Network 
estimate the distance from the DSN Tracking Stations to 
a spacecraft by transmitting to the spacecraft a square- 
wave range code of approximately 2 /is period, and then 
measuring the phase of this square wave as it is returned 
from the spacecraft relative to a time-delayed and 
doppler-shifted copy of transmitted code (Ref. 1). The 
ambiguity inherent in the periodicity of this square wave 
is resolved by using lower frequency codes which are 
coherent with it. There is currently some interest in in- 
creasing the precision to which this measurement is per- 
formed, or at least in ensuring that the degree of precision 
of which the system is theoretically capable is achieved 
in practice. This article proposes a minor change to 
existing-type ranging systems which can help to achieve 
that goal. 

Figure 1 is an overall functional block diagram of the 
ranging operation. The delay through the station alone is 
measured via the zero-delay or calibration subsystem. The 
delay through the spacecraft transponder is separately 


measured prior to flight. The net delay through the trans- 
mission medium, from which is inferred the true range, is 
estimated as the delay value remaining when the station 
and spacecraft delays are subtracted from the overall 
measured delay. We are in all cases not measuring the 
delay per se, but are instead measuring the phase of the 
square-wave range code. The exact relationship between 
the measured phase of any nonsinusoidal waveform and 
its “delay” can be somewhat fuzzy. This fuzziness results 
if waveform dispersion occurs anywhere along the signal 
path. We can combat it by using only a sine wave for the 
highest frequency code component, or by performing all 
measurements using only the sine wave fundamental of 
this component. In the following, we review some of the 
ramifications of this possible change to the binary coded 
sequential ranging systems. 

II. Additive Noise Effects 

Let us presume that we wish to use only the funda- 
mental sine wave component of a range code for delay 
measurement. The first thing to determine is what effect 
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this will have on the measurement errors due to noise. 
The signal received, r(t), is assumed to he an ideal square 
wave, of power s, and delay-shift 8, plus a White noise 
component n(f) with one-sided spectral density N„. The 
reference signal is an ideal unit-amplitude (unit-energy) 
square wave, and correlation detection is employed at 
two reference phases to estimate 6. The two correlations 
obtained are 


^ ~ Ji, 

Y = f‘sQr(^V-j)-rMdv 


( 1 ) 


As an alternate technique, let 'us use unit-energy sinu- 
soids for the reference signals. The two correlations ob- 
tained are 


X' = fJy/2 sin 
Y = J^'v^cos ^— drj 


( 5 ) 


As before, the variance of both X' and Y' is N„ * t/2. The 
mean values of X' and Y' will be determined by expanding 
the received square-wave range code in its Fourier series 
form. 


In E)q. 1, t is the measuring integration time and T is the 
bit period of the code. The expected value of X and Y has 
peak value \/s‘ t and depends upon 6 as shown in Fig. 2. 

The variance of both X and Y is No • The estimated 
delay, ■?'is given by 


r(v) - "(*?) + - 7 sin |"i • (ij + 8)1 

(6) 

The two correlations are then approximately given by 


^ " |X| + |Y| ’ 2 

whenever X is positive. Following Goldstein (Ref. 1) we 
use 


X' = \/Ss • t • — • cos + n, 

Y' = \/^ • t • — * sin + n. 


( 7 ) 



which implies 


( 3 ) 


The estimate of 0 is now given by 
T 

8 = — arctan (Y,X) 

7T 




2 




|y| 

|X1 + ]Y| 


)■ 


+ 


/ -Ysgnx y ) 

V |X| + |Y| / / 


a.’ < 


NoP 

8st 


a.' > 


NoP 

16 St 


( 4 ) 


The estimator variance is at its maximum when either X 
or Y is zero, and at its minimum when |X| = |Y|. Equa- 
tion 4 will be our reference against which the variance of 
the estimated phase of a modified range code will be 
measured. 


Using the same technique as before in £q. 2-4, we find 


a,* = • t • (sin’ ^ + cos* fl) J 


N,P 

16st 


( 8 ) 


This surprising result indicates that by using only the 
fundamental component of the received range code, we 
achieve the same performance with respect to additive 
noise that we could achieve with the entire square wave, 
and we achieve such performance irrespective of the 
actual delay between the received and reference code 
waveforms. Clearly the harmonics of the square-wave 
range code are of no benefit to us for conventional ranging 
operation. They can in fact be detrimental as we shall see 
shortly. 
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III. Harmonic Phase-Shift Effects 

To understand the effects that harmonic phase shifts, 
or waveform dispersion, have upon the ranging operation, 
we need to refer to the overall system diagram of Fig. 1. 
The calibration-subsystem in the station is wideband and 
as close to delay-free as possible. As a result, when range 
measurement is performed in-station to calibrate station 
erjuipment delays, it is almost entirely the station equip- 
ment, the receiver, the exciter, modulator, and transmitter, 
and the ranging subsystem itself that establishes the de- 
tails of the range code waveform, and its harmonic struc- 
ture. The spacecraft transponder, on the other hand, has 
a limited bandwidth that passes only the fundamental and 
its third harmonic; in current transponders, the third har- 
monic i.s .somewhat attenuated, although not phase-shifted. 

The range to the spacecraft is measured as follows; first, 
the station delay is measured, with all station equipment 
in the configuration which will be used with the space- 
craft measurement. Next, the delay is measured through 
the Earth-to-spacecraft link and the spacecraft trans- 
ponder, and then the station delay and the previously 
measured spacecraft delay are subtracted from this value. 
(See Ref. 2 for more detail.) For error-accounting, the 
calibration subsystem is directly in the range-measurement 
path. We envision two cases: first, where the station 
equipment generates and processes a perfect square wave 
range code, and second, where some station equipment 
has changed in such a way that the harmonics of the range 
code are phase-shifted and/or attenuated relative to the 
fundamental. For both cases, the spacecraft transponder 
passes first and third harmonics of the range code with 
no relative phase shift, although with some, say 3 dB, 
attenuation of the third harmonic. Additive noise will be 
ignored. 

The range code received at the ranging subsystem via 
the station’s receiver may be represented as 

-i- -^sin^S^ • t -i- + ••• I (9) 

where we denote by o*/k, the gain, and tf*, the phase shift 
of the kth harmonic relative to the first. The two correla- 
tion values used to compute the delay value are 


- C-< cosf-rj-l- 2^ — COsffc- 

odd 

Y(n) = X(r, - T/2) 


— ' T+6i 




( 10 ) 


During most of its operation, the local code reference of 
the ranging subsystem is constrained to track so that 
X(t) = Y(t), and the real delay measurement is contained 
within the physical shifts accumulated by the local coder. 
This forces 

0 = X(t) -X(t- T/2) 
0^{cos(^-t)-cos(y-'--|)} 

( 11 ) 

kss ' 

odd 

-COs(k.^.r-i-tf»-j)} 

The true tracking point is at t = T/4. If we define ij as 
the offset from the true tracking point, i.e., = t — T/4, 

change variables in Eq. (11), and manipulate it with 
trigonometric identities, the tracking point is defined by: 

fcss.odd 

• sin • ^ + (12) 

or 

O^sin(j) •sin(^i,) -I- 

ksd.edd 

• |sin (k • • cos(»t) + cos(k • 

• sin (flt)| 

Denoting by [’ll, the greatest integer less than x, we can 

represent sin s= sin * (” 1) ^ for all odd 

k. Let us also assume that i] is smaU; that for any offset 
which is of genuine interest with respect to ranging, we 
can approximate 
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Then the offset from the true tracking position is given by 


V 


T ^ 


(14) 


Note that there are two obvious conditions which can 
make 7 = 0; if either Ot — 0, for all k, or if at = 0, for 
all k > 3. In both cases, the delay measured by the 
ranging system is the delay suffered by the range code 
fundamental and there is no fuzziness in the meaning of 
that delay. The second case corresponds to performing 
the range measurement with the code fundamental only. 

The details of the error in the ranging operation can 
now be determined. If the station equipment generates 
and processes a perfect square wave, then there is no 
error caused by the code harmonics in measuring station 
delays, and r/ = 0. Similarly, even though the transponder 
filters out harmonics so that at ^ 0 for k > 5, and a.t < 1, 
we still have = 0 so that 1; = 0 in measuring delays 
through the spacecraft also. If we have also a calibration 
of the spacecraft delays using a perfect square wave, the 
net error in measuring range to that spacecraft is nil. 


In an alternate scenario, let us assume that some har- 
monic distortion arises somewhere within the complex 
equipment of the tracking station. The delay measured 
for the ground system then suffers an error which is well 
described by Eq. (14). Let be the gain factor of the 
spacecraft alone at the third harmonic of the high fre- 
quency range code, and use /S* = 0 for k > 5 to indicate 
that the transponder does not pass the fifth and higher 
harmonics. Assume that the spacecraft adds no phase 
shift to the third harmonic, and let at, Bt as before denote 
station equipment amplitude and phase characteristics. 
Then the delays measured through the spacecraft are in 
error by 


1jt/c — 


— aiP, sin B, 
j* y 

1 + aaPz COS Ba 


(15) 


And if we assume that the spacecraft delay itself was 
calibrated with a reference station without distortion, the 
net error in measuring spacecraft range is i],/c ~ i?> or 



For a concrete example, consider a station system which 
has 0,1 = as = 1, o* = 0 for k > 7, and has 6a = 10 deg, 
Bs = —30 deg, and a spacecraft for which pa = 1/2. 
Here, the resultant range error is 



I 


sin 10 deg H — ^ sin 30 deg 
9 25 

1 + — cos 10 deg — ^ cos 30 deg 
3 5 


sin 10 deg 

1 + cos 10 deg 
6 


A 


RE = .008 * T 


(17) 


In current ranging systems, where T is approximately one 
microsecond, this resultant error is eight nanoseconds, or 
approximately 2-1/2 meters. Larger errors than this are 
quite plausible, considering the amount of equipment in 
the ranging signal path! 


IV. Ameliorating the Waveform Dispersion 
Effects 

There are a number of alternative things which can be 
done in response to the magnitude of the effects evi- 
denced in Eq. (16). First of all, if the overall error budget 
is of the order of 10 to 20 meters, it can simply be ignored, 
as reasonable care with the equipment can restrain RE 
well below this value. If, however, we are asking how 
accurate the spacecraft range can be measured, then more 
care must be used. As a second approach, we could design 
and control the parameters of all tracking stations of the 
Network so that phase-shifts imposed upon the harmonics 
of range code by the station equipment are held within 
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prespecified tolerances. The actual phase-shift tolerances 
allowable may be determined from Eq. (16) once the 
allowable range-error is established. A third alternative 
is to make the calibration subsystem itself mimic the 
transfer characteristics of the spacecraft transponder. In 
this way, the harmonic structure of the range code arriv- 
ing at the ranging subsystem via either the spacecraft or 
the calibration subsystem is identical, and the range error 
equation analogous to Eq. (16) is likewise ideally zero. 
Considering component tolerances, and in particular, 
variations between spacecraft transponders, there would 
in fact be some slight residual range error resulting. 

The fourth alternative is that one which is being pro- 
mulgated by this report, namely the use of the funda- 
mental sine wave component only of the range code. This 
can be accomplished by inserting an appropriate filter 
anywhere on the signal path within the station equip- 
ment that is common to both the spacecraft and calibra- 
tion measurements, or by use of a sine-wave reference 
waveform as was analyzed in part II. Use of the sine-wave 
reference seems impractical from the standpoint of avail- 
able equipment. Appropriate filtering of the received 
signal on the path from the receiver to the ranging sub- 
system can achieve the identical result. Filtering at the 
receive-end of the signal path will remove all coherent 
harmonics of the range code, not only those generated 
intentionally as part of the square-wave range code, but 
also those which arise from nonlinearities along the signal 
path. Filtering could also be applied to the uplink range 
code between the ranging subsystem and the modulator. 
If the average uplink modulation index is maintained, and 
no peak deviation limit is imposed, the ranging signal 
power is in fact slightly increased. Filtering of the uplink, 
however, does not eliminate the desirability of further 
filtering of the downlink, since nonlinearities in the rang- 
ing path beyond the uplink filter can still inject harmonics 
at strange phase-angles. The filter (or filters) are directly 
in the group-delay path of the ranging measurement and 
their delay stability is critical to range accuracy, even 
though the exact delay value is not. 

For the experiments to be described, the filter was 
inserted in the lO-MHz signal line from the receiver as 
shown in Fig. 3. This point is preferable to others since it 
is the last opportunity to remove distortion prior to cor- 
relation and therefore serves to offset any degradation 
occurring elsewhere in the system. In principle, a corre- 
lator using a pure sine-wave reference would be superior. 
However, because the present systems are all mechanized 
using digital techniques the external filter was selected as 
the most practical approach. 


The characteristics of the filter needed to be uniform 
across the band pass vrith a linear phase-frequency slope 
so as to ensure maximum group delay stability. However, 
once the band edge was reached, a rapid attenuation was 
desired to eliminate unwanted harmonics. It was decided 
to design for a cutoff frequency of ±I MHz and a rejec- 
tion of at least 60 dB in the stop band beyond ±1.5 MHz 
where the offending harmonics are located. 

A schematic diagram of the filter is shown in Figure 4 
(component values will be found in Table I). An elliptic 
function design was selected as providing the best band- 
pass relative to band-stop characteristics. While ripples 
exist both within and without the band pass, their magni- 
tude is sufficiently small that the performance of the 
filter should not be affected. Moreover, the steep attenua- 
tion characteristic obtainable with this design makes its 
use in this application attractive. 

The basic design consists of a low-pass filter, 1 MHz in 
bandwidth which has been translated to 10 MHz by 
resonating elements. Because certain, identifiable, com- 
ponents determine the attenuation characteristics while 
others control center frequency and symmetry, it is im- 
portant to select these critical components carefully. 

For example, Li to L4 and C5 to Cll comprise the 
basic lew-pass filter and therefore determine the band- 
pass/band-stop characteristics as well as the transition 
between. Because the placement of poles and zeros is 
critical to this design, it is important that actual com- 
ponent values be selected to accurately reflect the com- 
puted numbers. The procedure followed with inductors 
was to adjust the number .of turns with the aid of aii 
inductance bridge until the closest possible value was 
obtained. Using this method, inductances within one 
percent of the calculated values could be realized. With 
the capacitors, up to three were connected in parallel to 
bring the capacitance to within one percent of the com- 
puted value. 

Because of the high performance requirements and the 
critical characteristics of this filter noted earlier, it was 
necessary to utilize high Q circuit components. Glass 
capacitors were used throughout and the torroidal induc- 
tors were wound with the largest size wire which could 
reasonably be fitted on the ferrite cores. Measured Q’s 
for the finished inductors were in excess of 100. 

After collecting the components comprising the low- 
pass portion of the filter, it was necessary to translate the 
entire spectrum to a 10-MHz center frequency for use 
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with the ranging system. This was accomplished with 
Cl to C4 and L5 to Lll. Since the only purpose of these 
components is to adjust the center frequency, confor- 
mance with calculated values is less important. The pro- 
cedure followed was to pick an element within the low- 
pa.ss filter and adjust the translating component so that 
the pair was resonant at 10 MHz. Thus Cl was trimmed 
until LI and Cl were a .series tuned circuit at 10 MHz. 
Likewise, turns were changed on L8 until it formed a 
parallel tuned circuit with C8. In general, each of the 
11 tuned circuits (Ln, Cn) comprising the filter was 
adjusted individually. 

The element pairs were then assembled into a specially 
made copper module as shown in Figure 5. Because of 
the high attenuation required, it was deemed advisable 
to separate the various sections into four compartments to 
minimize leakage. An effort was made to place the com- 
ponents in the center of the cell to reduce stray capaci- 
tance; the module itself sen'ed as circuit ground. 

Following assembly, tests were made to determine the 
bandwidth of the filter, and group delay and phase delay 
characteristics. Figure 6 is a photograph of the amplitude 
response as measured with a General Radio model 1711 
Sweep Generator and associated Model 1714 Display 
Unit. A sweep speed of 500 kHz per centimeter was used 
so as to provide maximum resolution in the frequency 
band from 8 through 12 MHz. 

The response appears to be reasonably symmetrical 
about the 10-MHz center frequency particularly at ±500 
kHz where the first range code sidebands are located. 
Beyond approximately ±900 kHz the filter begins to 
reduce the signal rapidly. At ±1.5 MHz where the third 
range code harmonic is situated, the attenuation exceeds 
60 dB. While there is a slight asymmetry on the high 
frequency side of the curve, between about 1.2 and 1.4 
MHz, it is considered unimportant since the design goal 
had been met at the ±1.5 MHz frequencies. 

Figure 7 shows the phase delay over the same band of 
frequencies. Because of the measurement technique, the 
figure appears to have a sawtooth pattern. In fact, the 
heavier lines, sloping down to the right, should be treated 
as a single continuous, unsegmented line. Ideally, the 
phase change should be linear with varying frequency. 
Practically, this is difficult to obtain in a multipole filter 
such as this one. In this case, the phase-frequency rela- 
tionship is reasonably linear over a range sufficient to 


contain the first range code sidebands. Since the code 
spectrum is discrete, nonlinearities occurring beyond 10 
MHz ±500 kHz should not materially affect the group 
delay. 

Filter group delay appears in Figure 8. Note that the 
filter has an insertion delay of approximately 600 nano- 
seconds at 10 MHz. This increases to about 700 nano- 
seconds for the lower range code sideband (9.5 MHz) and 
to about 900 nanoseconds for the upper sideband (10.5 
MHz). The average sideband delay compares well with 
the 808 nanosecond delay actually measured by the 
ranging system. 

At this point the intrinsic problem with filters is ob- 
vious. Where large attenuations and steep transitions are 
required, a design meeting these conditions will also 
include a rather substantial group delay. With large 
delays come the concomitant problems with stability 
which may introduce errors considerably in excess of 
those sought to be eliminated. The underlying problem, 
of course, is one of trying to resolve a phase comparison 
to approximately 0.1 percent. Clearly, this defect could 
be substantially reduced by widening the bandwidths and 
raising the code frequency. However, since the intent 
here was merely to investigate the feasibility of a filter as 
a solution to the harmonic phase shift problem, little 
attention was devoted to stabilizing the unit for field 
operations. 


V. Ranging Test Results 

To test the effectiveness of the filter as well as the 
hypothesis that much of the ranging inaccuracy is caused 
by nonuniform phase shifts in the higher order range code 
harmonics, the system was connected as shown in Figure 
3. First, the system delay was measured using a zero delay 
device. Second, a Mariner Venus-Mercury transponder 
was substituted for the zero delay device and the total 
delay remeasured. The difference between this and the 
prior measurement represents the transponder delay. 
Thereafter, the phase modulator in the ground exciter 
system was changed and the entire procedure repeated. 
In principle, the transponder delay should remain un- 
changed. However, because of the way in which the 
modulator and spacecraft ranging channel act upon the 
higher order range code harmonics, the transponder delay 
undergoes an apparent change. Installation of the filter 
served to remove the range code harmonics so that corre- 
lation is only with the code fundamental. 
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Results of these tests are shown in Figure 9. Trans- 
ponder delays were measured at 1-dB and 9-dB carrier 
suppression. Two ground system modulators were em- 
ployed. Number 2002 exhibited fairly normal operating 
characteristics. However, number 2034, while allegedly 
meeting Network test specifications, had been observed 
to exhibit unusual behavior such as low sensitivity. 

It can be seen from the figure that the transponder 
delay varies greatly as the modulation index is changed 
on number 2034. This distortion is substantially greater at 
low modulation levels than it is at high levels. Note that 
with the filter installed, the difference is reduced to less 
than one-third of its former value at small modulation 
angles. As the carrier suppression is increased to 9 dB the 
differences have disappeared within the uncertainty of 
the measurement as indicated by the brackets. However, 
while making a dramatic improvement, the effect is still 
present. Since it seems unlikely that any higher order 
harmonics are passing through the filter, the explanation 
must lie elsewhere. At this juncture it is uncertain why 


divergence exists and further tests will be necessary to 
develop a better understanding of the causal relationships. 


VI. Summary and Future Plans 

Our preliminary analysis, described herein, has con- 
firmed the existence of a real, or at least a potential, limi- 
tation to ranging accuracy resulting from phase shifts of 
the higher harmonics of the square-wave range code. Ex- 
perimental and analytical evidence has been developed 
which shows that this limitation can be largely eliminated 
by filtering the range code to its fundamental component 
only at the receiver. 

Additional testing should be performed to more exactly 
quantify the improvement to be obtained in this manner. 
Also, since that filter is now itself a critical component of 
the station’s group delay path, the stability of that filter 
under temperature and other environmental variables 
must be determined. 
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Table 1. Filter componentt 


Component 

designation 

Computed 

value 

Measured 

value 

Measured 

Q 

LI 

5.11 /<H 

4.95 pH 

185 

L2 

4.17 >>H 

4.27 pH 

195 

L3 

5.09 /iH 

4.90 pH 

195 

L4 

4.75 ;iH 

4.60 pH 

195 

L5 

0.792 

0.790 pH 

185 

L6 

0.259 

0.262 pH 

185 

L7 

0.395 pH 

0.400 pH 

175 

L8 

0.157 pH 

0.157 pH 

140 

L9 

0.U5 pH 

0.116 pH 

130 

LIO 

0.137 pH 

0.136 pH 

130 

Lll 

0.127 pH 

0.126 pH 

140 

Cl 

50 pF 

50 pF 


C2 

60.7 pF 

61 pF 


C3 

49.7 pF 

50 pF 


C4 

S3 pF 

53 pF 


C5 

320 pF 

320 pF 


C6 

977 pF 

980 pF 


C7 

641 pF 

640 pF 


C8 

1612 pF 

1612 pF 


C9 

2198 pF 

2200 pF 


CIO 

1849 pF 

1847 pF 


Cll 

1989 pF 

1980 pF 
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Fig. 2. Ranging eorretetion values for aquan-wave raferenea Fig. 4. Ranging code fitter schematic 
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Fig. 5. Ranging code filter module 



SWEEP SPEED. HORIZ. SOO KHz per cm 
VERT. 50 deg per cie 


SWEEP SPEED: HORIZ. 500 Ic Hz per cm 
VERT. 10 dB per cm 


CENTER FREQUENCY 10 MHz 


CENTER FREQUENCY 10 MHz 


Fig. 6. Filter amplitude response 


Hg. 7. HIter phase delay 
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SWEEP SPEED. HORIZ. 500 KHz per cm 
VERT. 500 ns per cm 

CENTER FREQUENCY 10 MHz 

Fig. 8. Filter group delay 


PT = -II5d8m 



UPUNK CARRIER SUPPRESSION (dB) 


Fig. 9. Filter insertion results 
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Multipath Tests on 64-m Antennas Using 
the Viking Orbiter-1 and -2 Spacecraft as 
Far-Field Illuminators 

T. Y. Otoshi 

Communications Elements Research Section 

0. L Brunn 

R. F. Systems Development Section 


Far-field multipath tests were performed on the 64-m antennas at Goldstone, 
DSS 14, and Madrid, DSS 63, by use of the transponders on the Viking Orbiter-1 
and -2 spacecraft. At the time of the tests. Viking Orbiter-1 and -2 spacecraft 
were in their interplanetary orbits to Mars and were respectively about 21.9 
million and 13.9 million km from Earth. The test results showed that the effects 
of multipath in the far-field of the 64-m antenna were to cause less than a 5-ns 
peak-to-peak variation on two-way range and 0.1-dB peak-to-peak variations on 
received signal level. The multipath signal level was calculated to be 
approximately 40 dB weaker than the primary signal in the far-field main beam 
direction. 


I. Introduction 

In a previous article (Ref. 1). it was shown that 
multipath eflPects on 64-m antennas could be the cause of 
the 15-m range residual observed between DSS 43 and 
DSS 63 during Mariner 10 Mercury Encounter 1 on March 
29, 1974. To support this conclusioi^ however, it was 
necessary to prove that on the 64-m antenna, the 
multipath error occurs only in ground station delay 
calibrations and does not occur in a far-field range 
measurement to a spacecraft. There was, therefore, 
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considerable interest in performing tests to show that 
multipath effects in the far-field were negligibly small. 

In practice it is difficult to perform a far-field ranging 
test on the 64-m antenna because, in order to be in the 
far-field of the 64-m antenna at S-band, the transponder 
must be located at least 62.8 km (39 mi) away and 5.6 km 
(3.5 mi) high for a minimum elevation angle of 5 degrees. 
Due to the non-availability of suitable collimation toyvers, 
radiosonde balloons, or satellites for ranging tests, it was 
necessary to use a spacecraft located in the far-field. A 
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spacecraft t^t is sufficiently far from Elarth will appear 
like !L point source that radiates constant power. Since the 
ground station tracking antenna can accurately follow any 
movement of the spacecraft (due to the spacecraft orbit 
and Elarth’s rotation), the test conditions are nearly 
equivalent to those which would be created by a 
stationary far-held collimation tower in the sky. 

The first known far-field test of a large antenna using a 
spacecraft was reported by Levy et al. (Ref. 2) in 1967. In 
this test the Surveyor spacecraft on the moon was used as 
the far-field illuminator. One of the interesting test results 
obtained was the focus curve shown in Fig. 1. The focus 
curve shows the relationship of far-field gain changes as a 
function of subreflector position relative to nominal. It 
was shown by Potter (Ref. 3) that the ripples observed on 
the focus curve in Fig. 1 were caused by a multipath 
component originating from a reflection from the horn 
aperture. This test result, although restricted to signal 
level data, presented the first known evidence that a small 
multipath effect existed in the far-field of the 64-m 
anteniut. 

In July 1975, a series of far-field multipath tests was 
performed on the 64-m antenna at DSS 14 using the 
ranging transponder on the Helios spacecraft as the far- 
field illuminator. However, due to the long interplanetary 
media of about 2 AU and because Helios was only 7 
degrees from the Sun as seen from Earth, the differenced 
range versus integrated doppler (DRVID) and signal level 
data were excessively noisy. Therefore, no conclusions 
could be made on far-field multipath effects. 

A unique opportunity to perform far-field tests with a 
spacecraft under ideal conditions came about two months 
after the launching of the Viking Orbiter (VO) -1 and -2 
spacecraft. It is the purpose of this article to present the 
far-field results obtained with the Viking spacecraft as well 
as to compare the experimental data with theoretical data. 


II. Theoretical Results 

As was shown by Potter (Ref. 3), the focus curve for a 
Cassegrain antenna has a parabolic shape for defocusing 
losses of less than 1 dB. The approximate expression for 
the focus curve is 

CdB(*) “ (Co)dB “ ■A2(* - *o)* (1) 

where 

G^{x) ^ antenna gain, dB 

(Co)dB = antenna gain at maximum, dB 


A 2 = constant 

* = focus position 

Xq = focus position for maximum gain 

Figure 2 shows the focus curve of the 64-m tricone system 
at 2.295 GHz. This curve as based on Eq. (1) has a value 
for A 2 of 0.02126 dB/cm^ (0.13716 dB/in.2) for the 64-m 
tricone configuration (see Potter, Ref. 3). If multipath 
effects are present in the far-field, they will introduce 
ripples in the focus curve similar to that shown in Fig. 1. 

Although there is a similar focus curve for the uplink 
frequency of 2.113 GHz, the spacecraft radio system 
removes the uplink signal level variations and transmits 
back a constant coherent downlink signal back to Earth. 
Therefore, in a far-field multipath test, one only needs to 
consider defocusing losses on the downlink signal at 
approximately 2.295 GHz. 

To the authors’ knowledge, no published information is 
available concerning the far-field group delay and phase 
delay changes resulting from the defocusing of the 
subreflector on a Cassegrain antenna. A defocused 
subreflector computer program written by Sorensen (Ref. 
4) was modified by Rusch to yield far-field data for a 
defocused Cassegrain antenna. The phase versus frequency 
and phase versus subreflector position data were used to 
compute far-field group and phase delay changes, 
respectively. Figure 3 shows the group delay change' due 
to defocusing of a symmetrical Cassegrain antenna in the 
absence of multipath. The primary feed was assumed to 
have an E- and H-plane power pattern of (43.63) (20 logjo 
cos y), which is nearly equivalent to the main lobe power 
pattern for the Potter horn described in Refs. 5 and 6. The 
angle y is the angle from boresight of the primary feed. 
Although the present 64-m antenna has a corrugated 
feedhorn and a tricone configuration, it is not expected 
that the delay change curve for the present configuration 
would differ significantly from that shown in Fig. 3. For 
comparison purposes, the dashed curve shows the total 
(uplink plus downlink) delay change relationship for a 
single ray traveling only along the axis of the subreflector. 
It can be seen that the total change in the far-field due to 
contributions from the entire aperture is slightly less than 
the delay change for the single on-axis ray. 

From the Sorensen /Rusch Program, it was found that 
the phase delay change was equal to the group delay 


’ 1 he Sorensen/Rusch rrugrani resuits shewed tii£ piMtSC {icisy 
change was the same as group delay change as a function of 
subreflector position. Therefore, Fig. 3 applies to either phase or 
group delay changes. 
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change. Since, by definition, DRVID is group delay 
change minus phase delay change (Ref. 7), the far-field 
DRVID data should show no change as a function of 
subreflector position. However, if there is a multipath 
effect in the far-field, the group delay changes will be 
much larger than phase delay changes and a cyclical 
variation on DRVID data will be observed due to 
subreflector movement. 

III. Experimental Results 

A. VO-2 Tests at DSS 14 

On October 23, 1975, a far-field test was performed 
with the VO-2 spacecraft at DSS 14. At the time of the 
test the VO-2 spacecraft was approximately 13.9 million 
km from Earth. A strong received signal level of -136 dBm 
and a good ranging power-to-noise ratio of 33 dB were 
prevalent at DSS 14 during the far-field ranging tests. 
Other favorable conditions for this test were the short 
round-trip light time (RTLT) of 93 seconds and the 
absence of any solar and ionospheric effects introduced 
into the DRVID data. A short RTLT was desirable for the 
far-field multipath test because a short RTLT is equivalent 
to a short interplanetary medium. A shorter interplanetary 
medium generally has less electron content than a longer 
one, and, therefore, the DRVID data will be less likely to 
change due to electron density changes. No ionospheric or 
solar effects were introduced because the tests at DSS 14 
were done during night-time hours. 

The far-field multipath test consisted of ranging to the 
spacecraft and recording DRVID and received signal level 
data as a function of subreflector positions. For this test, 
the subreflector was moved in 1.27-cm (0.5-in.) increments 
over a total distance of 15.24 cm (6 in.). The test began 
and ended with the subreflector at the nominal setting so 
that one could observe and remove the long-term drift in 
the test data. 

The experimental test results may be seen in Figs. 4-7. 
Figure 4 shows the results of taking 4 hours of DRVID 
data while moving the subreflector and 2 hoius of data 
while keeping the subreflector stationary for comparison 
purposes. The drift in the DRVID data during the 4-hour 
multipath test was about 10 ns. After removing this drift 
by fitting a second-order curve to the experimental data, 
the residual differences at each subreflector position were 
then averaged. A least squares fit of the multipath 
equation was then made to the averaged DRVID points 
by use of the Multipath Computer Program described in 
Ref. 1. The resulting plot shown in Fig. 5 shows relatively 
good correlation between calculated and experimental 


values. The peak-to-peak change in range delay was about 
3 ns. The Multipath Program calculated the multipath 
signal to be -45 dB relative to the primary signal in the 
far-field. This relative multipath signal strength will 
produce a ripple of 0.1 dB in the downlink received signal 
level focus curve. 

Figure 6 shows the downlink received signal level data 
obtained during the far-field multipath tests. A comparison 
of experimental and theoretical results may be seen in Fig. 
7. The experimental points were obtained from averaging 
the received signal level values at each subreflector setting 
and normalizing them to the received signal level value at 
the 0 in. position.^ A least squares parabolic curve was 
then fitted to the experimental data points. For compari- 
son purposes, the peak of the tricone theoretical curve 
(Fig. 2) was aligned with the peak of the exp>erimental 
curve. It can be seen that excellent agreement was 
obtained between theory and experiment. The deviations 
of the data about the experimental curve are consistent 
with the deviations caused by a multipath signal being -45 
dB with respect to the principal wave. 

For purposes of comparison, the results of pre-cals done 
with the zero delay device at DSS 14 are shown in Fig. 8. 
It can be seen that the peak-to-peak variations are 40 ns 
and 3 dB for range delay and received signal level, 
respectively. One explanation as to why a small multipath 
effect occurs in the far-field measurements but a large 
effect occurs in zero delay calibrations is the following: 
The primary waves impinging on the paraboloidal dish 
surface add up in phase in the far-field. Although the 
multipath wave signal strengths at a localized area can be 
as strong as -15 dB relative to that of the primary wave, 
the multipath waves are not the same in phase or 
amplitude over the entire dish surface. Therefore, when 
the multipath components combine in the far-field, they 
do not necessarily add in phase. Since cancellations can 
occur, it is not unreasonable that the net far-field 
multipath signal level can be -45 dB relative to the 
primary wave in the main beam. 

B. VO-1 Tests at DSS 63 

On October 26, 1975, a far-field multipath test was also 
performed at DSS 63. This test was similar to the DSS 14 
test but differed in that the VO-1 spacecraft was used 
instead of VO-2. The RTLT was 146 seconds, which 
indicates that the VO-1 spacecraft was approximately 21.9 

^The subreflector position on the equipment is indicated in inches rather 
than in centimeters. In order for data to be meaningful and be useful 
to the user and also to avoid double conversions, the units should be 
kept in inches as reported. 
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million km from Earth. VO-1 was farther from Earth than 
VO-2 was; therefore, the received signal level was weaker. 
The ranging power-to-noise ratio was 22 dB as compared 
to 33 dB obtained with VO-2. In order to eliminate 
ionospheric and solar effects, the test at DSS 63 also was 
done during night-time hours. 

Although far-field data had already been successfully 
obtained at DSS 14 with the VO-2 spacecraft, the test was 
performed at DSS 63 because the DSS 63 64-m antenna 
was known to have the largest multipath effect (Ref. 1) 
and also had a cone missing. These differences in ground 
station configurations could cause far-field multipath 
effrcts to differ. 

Figure 9 shows the DRVID data during the 3-hour 
period that the subrefiector was being moved. It may be 
seen that there was virtually no drift in the DRVID data 
during the test. Because this was an abbreviated test, 
DRVID stability data with the subreflector stationary 
were not obtained during the same VO-1 pass. However, 
examination of DRVID stability data at DSS 63 for a VO-2 
pass three days later show^ that the DRVID drift 
between 0000 GMT and 0700 GMT was only 10 ns. 
Therefore, it can be reasonably assumed that the drift 
during the DSS 63 VO-1 test was also small. 

Figure 10 shows a comparison of the averaged 
experimental DRVID data at each subreflector setting and 
a best fit theoretical curve which is based on the 
multipath equations (Ref. 1). Although there is only partial 
agreement between theory and experiment, the peak-to- 
peak variation of range delay is only about 5 ns. The 


multipath signal was calculated to be 39 dB weaker than 
the primary signal in the far-field. 

Figure 11 shows the experimental received signal leveb 
during the multipath test. In Fig. 12, which shows the 
normalized plot, each daU point is the result of averaging 
data at each subreflector position. A comparison of a least 
squares parabolic curve (fitted to the experimental data 
points) and the theoretical curve for the tricone configu- 
ration may be seen in Fig. 12. It can be seen that the 
agreement between the experimental parabolic curve and 
the theoretical curve is reasonably good. Although the 
agreement between theory and experiment is not as good 
as that obtained in the DSS 14 VO-2 tests, it should be 
pointed out that at DSS 63, the 64-m antenna has a cone 
missing. In addition, the received signal level with VO-1 
was 10 dB lower than that on the VO-2 tests, and 
therefore, the data were noisier. 

IV. Conclusions 

Good experimental results were obtained ftom the far- 
field multipath tests on the 64-m antennas at DSS 14 and 
DSS 63 and the two Vi^g spacecraft The far-field 
multipath effects were found to be smaU. The variations 
due to subreflector movement were typically less than 5 
ns peak-to-peak on total range delay, and less than 0.1-dB 
cyclic variation on received signal level. 

The excellent results obtained from Viking spacecraft 
can be attributed to a good ranging power-to-noise ratio 
(>^ dB). a short RTLT (<3 min), and the absence of solar 
or ionospheric effects. 
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Fig. 1. 64 hd antenna Surveyor focus curve (from Ref. 2) 



SUBREFLECTOR POSITION RELATIVE TO NOMINAL, in. 


Fig. 2. Far-field focus curve from 64-m Cassegrain antenna 
tricone configuration in the absence of multipath 



Fig. 3. Total uplink plus dmmlink phase or group delay change 
in the far-fwld of a defoeused 64-m Cassegrain antenna in the 
absence of multipath 
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ng. 5. Itaduced ORVID data from DSS 14 VO-2 apacacraft 
multipath test on 1975 GMT Day 296 


EXPERIMENTAL 

BEST FIT PARABOLA 

THEORETICAL TRKONE 

FOCUS CURVE 


•2 -10 1 2 
INDICATED SUBREFUCTOR POSITION, in. 


Fig. 7. Theoretical and experimental results for relative received 
signal level on DSS 14 VO-2 multipath test on 1975 GMT Day 
296 
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Fig. 8. Theoretical •ndnperimentalrasuMs tar muKIpith tests on 
DSS 14 VO-2 Pass 44 pre^al on 1975 GMT Day 296. POUrlsation 
is RCP. 
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Fig. 9. Original DRVID date from DSS 63 VO-1 spacecraft 
muHipath test on 1975 GMT Day 299. Polarization Is RCP. 
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Fig. 10. Reduced DRVID data from DSS 63 VO-1 spacecraft 
multipath test on 1975 GMT Day 299 
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Fig. 11. Recaived •ign*l level during DSS 63 VO-1 multipath 
test on 1975 GMT Day 299. Polarization is RCP. 



INOICATEO SUBREFLECTOR POSITION, m. 


Fig. 12. Theoretical and experimental results for relative re- 
ceived signal level on DSS 63 VO-1 multipath test on 1975 GMT 
Day299 
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Design of a Superconducting Cavity 
Stabilized Maser Oscillator 

W. Higa 

Communications Elements Research Section 


The pioneering works of W. H. Hartwig (Ref. 1), S. R. Stein and J. Tumeaure 
(Ref. 2). J. J. Jimenez and A. Septier (Ref. 3), and others (Ref. 4) have shown the 
possibility of using a superconducting cavity to stabilize a microwave oscillator. 
The achievement of cavity Qs of the order of has made possible the 
realization of frequency standards with performance which could surpass that 
of the hydrogen maser. 

The present study explores the possibility of integrating a superconducting 
cavity with a traveling wave maser to obtain a frequency standard with very 
high spectral purity. 


I. Introduction 

There are two basic methods for using a high Q 
resonator to stabilize a microwave oscillator. In one 
method, a voltage controlled oscillator (VCO) is used in a 
FM detector loop, as shown in Fig. 1, to stabilize the 
oscillator frequency. The second method, shown in Fig. 2, 
utilizes the high Q resonator as a transmission filter and 
stabilizes a microwave oscillator to the cavity/frequency. 
Both methods have been used by the investigators 
mentioned in Refs. 1-4. 

The motivation for the present program was the long 
experience with traveling wave masers (TWMs) and 
closed-cycle refrigerators (CCRs), which technologies 
appeared to blend in very naturally with the requirements 
for a good superconducting cavity stabilized maser 
oscillator (SCSMO). 


In the present article, a discussion is presented of the 
microwave electronics problems related to the SCSMO. In 
future articles, the cryogenics problems will be discussed 
and the statistical properties of the output signal from a 
SCSMO wiU be described. 


II. Desiga of a SCSMO 

The low-noise performance of a TWM would appear to 
make it a good device for integrating with a superconduct- 
ing cavity to achieve a stable oscillator.* In such an 
application the TWM would operate at high signal levels 


*S. B. Stein is investigatiDg the use of parametric amplifiers for reasons 
similar to ours (Ref. 5). 
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and in partial gain saturation in order to provide 
amplitude stabilization. Thus, it is apparent that the 
application of the TWM here is completely different from 
the usual linear-uibaturated mode in a receiving system. 


In the usual electron tube or semiconductor oscillator, 
the mechanism for amplitude stabilization is contained in 
the Van der Pol equation (Ref. 6), which introduces 
damping terms proportional to the square of an exciting 
voltage or current Whereas the Van der Pol equation was 
derived through phenomenological arguments to explain 
experimental observations in electronic oscillators, the 
same saturation effects were already inherent in the 
equations for the solid state maser as originally proposed 
by Bloembergen (Ref. 7). Subsequently, it was observed 
that power saturation in a maser was necessarily 
accompanied by a rise in the spin temperature of the 
amplifying medium, and that the spin temperature 
determined the noise performance of a maser. 


Siegman (Ref. 8) has postulated an equation to account 
for saturation in a TWM; however, for present purposes a 
graphical approach will be adequate. Figure 3 shows a 
typical response for a TWM as a fnnction of distance 
aJong an amplifying structure of great length; typical 
operating points a and b for an amplifier of length L cm 
are shown. The small signal gain C of the TWM is given 
by the slope for small values of x, while the oscillator gain 
Co is given by P^/Fo ^uid is, of course, less than G. The 
response shown in Fig. 3 is for fixed conditions on certain 
parameters, such as pump power, operating temperature, 
and the static magnetic field. Variations in these 
parameters will cause corresponding power fluctuations in 
the oscillator output. Moreover, to the extent that these 
variations also cause frequency and phase fluctuations, 
they will need to be examined in greater detail in a later 
study. 


An important observation can, however, be made at this 
point: the low-noise performance of a TWM can be 
realized in an oscillator provided the input section of the 
maser operates in the linear-unsaturated mode. The signal- 
to-noise ratio then can be made sufficiently high by the 
time the signal sufifers amplitude saturation, toward the 
end of the TWM, that the spectral purity is not degraded. 
In other words, the high spin temperature (which implies 
high noise temperature) toward the end of the structure, 
due to saturation, will not degrade a properly designed 
SCSMO. 


The microwave design problem is then to provide the 
proper coupling coefficients to the superconducting cavity 
so as to achieve the desired gain Gg while maxim izi n g the 
loaded Q. 

Figure 4 shows the equivalent circuit for the SCSMO. 
The transmission cavity has been analyzed in great detail 
by Montgomery, Dicke and Purcell (Ref. 9) and only the 
results need to be given here. 

The transmission loss throu^ the cavity at resonance is 
given by 

7ta>o) - (4^,^2)/(1 + + /5 s) 

and the loaded Q is given by 

Ql ~ + fi2> 

In the above equations, 

7(ci>o) transmission loss at resonance 
fit ~ coupling coefficient at input 
~ coupling coefficient at output 
“ loaded Q of cavity 
Qg unloaded Q of cavity 

An X-band TWM can easily have a gain in excess of 40 
dB; hence, an oscillator gain of 30 dB is reasonable. Thus 
the transmission loss through the cavity should be arotmd 
10~^. Assuming the input coupling can be made equal to 
the output coupling, 

Pi ~ P2 ~ P 

one has, finally, 

Ttwo) - 4/5V(l 2fi) 
and 

Ql ~ QAi + ZP) 

Using the numbers quoted above, 

Ql - (inxa)Qu 

Thus, the high gain in a TWM makes possible a loaded Q 
which approaches the unloaded Q. 

One of the advantages of the servoed oscillator shown in 
Fig. 1 is that the time constant of the loop may be made 
sufficiently long to smooth out the noise in the system. 
With the SCSMO the noise is negligible, and the cavity 
performs a smoothing of any phase instabilities which may 
be inherent in the TWM. For Qs of the order of 10*° and 


JPL DEEP SPACE NETWORK PROGRESS REPORT 42-31 


SI 



oscillations at a frequency of 10“’ Hz, the time constant of 
the cavity is around 1 second The transit time throu^ a 
TWM is something less than a microsecond and the 
superconducting cavity thus averages a signal which 
makes over a million round trips through the TWM. Small 
phase fluctuations throu^ the TWM are nullified and a 
hi^y monochromatic signal should be realizable. 

Ttiineaure and Stein (Ref. 2) have shown that the Q of a 
superconducting cavity can increase by as much as an 
order of magnitude per Kelvin in the region of 4 to 1.5 K. 
It is essential, therefore, to operate the cavity at around 
1.5 K; the TWM may be operated at 4J2 K. It is apparent 
then that the difficult problem to be solved for a 
continuously operating SCSMO is that of cryogenic 
refrigeration. Not only is it necessary to maintain the 
cavity at around 1.5 K with high temperature stability but 
the cavity also needs to be kept free of any mechanical 


motion. These problems will be discussed in a future 
study. 

III. Conclusions 

The superconducting cavity appears to be compatible 
with a traveling wave maser. The real advantage of the 
superconducting cavity stabilized maser oscillator is that 
the cavity can be integrated with the TWM in such a way 
that rigid microwave coimections can be made. Second- 
order perturbations resulting from fluctuations in transmis- 
sion line joints are thus eliminated However, the more 
serious problem is that of the mechanical (dimensional) 
stability of the cavity itself. 

A special closed-cycle refrigerator needs to be devel- 
oped for the SCSMO, and it will be discussed in the next 
report. 
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Fig. 1. Frequency modulation aervoed oscillator. (The supercon- 
ducting cavity is used as a frequency discriminator.) 



ng. 3. Typical power saturation curve for a TWM (see text) 




Fig. 4. Equivalent circuit for a SCSMO (A small part of output 
signal would be used as the stable oscillator output) 


Fig. 2. Regenerative microwave os- 
dUator (The cavity is used as a 
transmission filter.) 
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QVLBI Doppler Demonstrations Conducted 
During Pioneer 11 Encounter and Solar 
Conjunctions 

C. C. Chao 

Tracking and Orbit Determination Section 


During the Jupiter encounter of Pioneer 11 in December 1974, a limited 
amount of simultaneous two-way and three-way doppler was obtained. It was 
demonstrated in this study that, based on a very short arc (two weeks) of 
differenced doppler (quasi very long baseline interferometry or QVLBI doppler), 
the B-plane predictions were as good as 200 km at E-1 day and 400 km at E-3 
days. 

During the two solar conjunctions of Pioneer 10 and 11, which occurred in 
March 1975, several passes of QVLBI doppler were obtained and analyzed. The 
data quality of QVLBI doppler was found to be two to three times better than 
the conventional two-way doppler when the Sun-Earth-probe angle was less 
than 5 degrees. 


I. Introduction 

Studies in Ref. 1 have shown that using difierenced two- 
way and three-way (QVLBI) doppler data rather than the 
conventional two-way doppler data may significantly 
improve the accuracy with which the trajectory of a 
spacecraft can be determined in the presence of process 
noise. Process noise describes the efifect of unmodeled 
spacecraft accelerations due to attitude control gas leaks 
and solar radiation pressure anomalies, etc. This data type 
has been successfully demonstrated using the radio metric 


data of Pioneer 10 (Ref. 2) and that of Mariner 10 (Ref. 3). 
The results of the Pioneer 10 demonstration indicate an 
order of magnitude improvement at Jupiter B-plane in the 
presence of large unmodeled accelerations (10~^ km/sec^) 
due to the four massive Galilean satellites. The results of 
the later demonstration conducted with Mariner 10 
conclusively showed that this new data type can reduce 
the effects of moderate unmodeled spacecraft accelera- 
tions (10*^° km/sec^) by an order of magnitude and reduce 
the effects of the solar corona by a factor of five when the 
Sun-Earth-probe (SEP) angle is less than 5 degrees. 
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Unmodeled spacecraft accelerations and solar plasma 
effects will be two major difficulties in navigation of future 
space missions such as Mariner Jupiter-Satum and the 
Saturn flyby of Pioneer 11. The latter effect is even more 
important because solar conjunction will occur during 
many of the critical periods of planet encounter of 
Mariner Jupiter-Satum and Pioneer 11 missions. The 
order of magnitude increase in data noise due to the solar 
corona makes the orbit determination (OD) extremely 
difficult. During the solar conjunction of Mariner 10, the 
doppler data noise became so large (1.1 Hz) that neither 
the S-band-X-band (S-X) dual-frequency technique nor the 
sequential estimation method was able to calibrate or 
remove this effect. Therefore, the OD strategy relied 
heavily on the solutions based on QVLBl doppler data. By 
comparing short arc solutions with that from conventional 
data, a factor of 4 improvement in ^plane uncertainties 
was found (Ref. 3). 

During the Pioneer 11 encounter of Jupiter, one month 
of simultaneous two-way and three-way doppler data were 
received. The two-way doppler data diuing this period 
were relatively clean, and the accelerations due to 
Galilean satellites had been correctly modeled. Therefore, 
we do not expect significant improvement in data quality 
from the QVLBl doppler data. Consequently, the Pioneer 
encounter demonstration emphasized the orbit determina- 
tion capability of this new data type. 

Later in the spring of 1975, both Pioneer 10 and 
Pioneer 11 were in solar conjimction, and a few passes of 
two-way and three-way doppler were obtained. Because of 
very limited data coverage, data quality improvement was 
demonstrated with no attempt at studying the OD 
improvement. 


il. Pioneer Encounter 

Similarly to previous demonstrations (Refs. 2,3), the 
results of orbit determination based on this new data type 
will be compared with those from conventional data 
types. In the meantime, both of them will be compared 
with the true values determined from both pre- and post- 
encounter data. The simultaneous two-way and three-way 
doppler data coverage was little less than a month during 
encounter, with 18 days before and 10 days after the 
closest approach of Jupiter. Figure 1 shows the data 
distribution of two-way and three-way simultaneous 
doppler for that period. On most of the days, the coverage 
was less than 2 hours, which is relatively poor compared 
with the continuous coverage of conventional two-way 
doppler data. 


During this period, the conventional doppler data was 
relatively quiet. Consequently, we do not anticipate 
significant improvement in data quality from the new data 
type as we Ad in previous demonstrations. Instead, we 
emphasize the actual capability for navigation of this two- 
way and three-way Afferenced doppler. Besides, the 
ephemerides of the four massive Galilean satellites, which 
were turned off in the previous Pioneer 10 demonstration 
(Ref. 2), were included in the orbit determination program 
this time. 

Three Afferent data types are included in this ai^ysis: 
the conventional two-way doppler (F2), the two-way and 
three-way Afferenced (QVLBl) doppler (F3C) and the 
ramped range (ETR). The selected data arc starts from 30 
days before Jupiter encounter, with conventional doppler 
loosely weighted (Table 1). Six ramped range data were 
evenly distributed in this interval, and they are weighted 
at 10 km. Three cases of OD solutions at three Afferent 
data arcs were obtained. In the first case, six state 
parameters and two parameters for the constant frequency 
bias for the DSS 14-43 baseline (AlOOl) and the DSS 43- 
63 baseline (A1002) were estimated. In the second case, 
two more parameters for the relative frequency Aift 
(A2001, A2002) were added to the estimated list. The third 
case is the state only solution without QVLBl doppler and 
wiA two-way doppler weighted normally at 0.045 Hz and 
ramped range (10 km). Table 2 shows all the estimated 
and considered parameters, together wiA a priori and 
nominal values used in the OD analysis. Those a priori 
values are standard values used by the Pioneer 10-11 
navigation team. 

The three data arcs give solutions at 5, 3 and 1 days 
before encounter. The B-plane values shown in Fig. 2 
inAcate an excellent agreement in the B * T component 
between that preActed by QVLBl doppler (solid and 
empty triangles) and the post-encounter best estimate (less 
than 60 km). The B * R component, which is the out-of- 
plane component and usually is difficult to preAct, has a 
spread of 1000 km in the solutions from QVLBl data. This 
is not surprising because at E-5 days, the QVLBl doppler 
only sparry covered an arc of about 10 days. At E)-3 
days, the B * R preAction missed the best estimate by less 
than 400 km. llie best preAction from QVLBl data was 
good to 200 km, which is at E-1 day with only a 16-day 
arc of QVLBl doppler. The preActions based on 
continuous two-way doppler (case 3) at E-5, E-3 and E-1 
day agree with the b^ estimate within 120 km. It is 
estimated that the QVLBl doppler data coverage actually 
obtained during this time span (E)-30 days to E)-l day) is 
less than 15 percent of the possible QVLBl doppler 
covered during the same time interval (based on 7-hour 
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two-way and three-way overlap per 24 hours). Had we had 
all the possible two-way and t^ee-way data during that 
time span, considerable improvement in B-plane accuracy 
would have been expected. 

The difference between the solid and empty triangles 
(Fig. 2) tends to indicate the effect due to the drift 
between frequency standards (new Rb maser 5065A). The 
estimated values for frequency biases are shown in Fig. 3. 
The bias between DSS 14 and DSS 43 was as large as 33.9 
mHz, which is consistent with the estimated value from 
the Mariner 10 QVLBI doppler demonstration (Ref. 3). 
The slope of the frequency ^ift between DSS 14 and DSS 
43 for this period is unusually small at 0.1 mHz per month. 
The slope ^tween DSS 43 and DSS 63 is also small during 
this interval. It is seen from Fig. 3 that the estimated 
values of bias (case 1 solution) differed from the values 
when slopes were estimated at the same time (case 2 
solutions). This is true particularly for the £-5 day 
solution. It may be responsible for the large error in 
B-plane predictions. When both bias and slope are 
estimated from a short arc of data, the two parameters are 
highly correlated and thus degrade the uncertainty of 
estimated biases. By comparing results (covariance) of case 
1 and case 2, the bias may be determined from 2 weeks of 
QVLBI data within a few tenths of a mHz if slopes are not 
estimated and station location has no error, while in case 
2, when slopes are estimated, the uncertainty of estimated 
bias becomes as large as 20 mHz. The above fact strongly 
suggests that better frequency standards such as H-masers 
with long-term drift less than one tenth of a mHz or better 
are highly desirable for QVLBI doppler tracking. If we 
had H-masers as frequency standards, with the constant 
bias properly removed either by estimating in the Orbit 
Determination Program or other independent methods 
such as VLBI or SITT*, the B-plane solution would be 
improved significantly. 

Figure 4 shows the sensitivity of B-plane components to 
various parameters for the three different data arcs. It is as 
expected that the QVLBI data type is much less sensitive 
to errors in the mass and ephemerides of Jupiter than the 
conventional data type. QVLBI doppler is more sensitive 
to errors in station location and frequency standards. 

As a research effort, only the simultaneous two-way and 
three-way doppler data (without differencing) from E-18 
days to E-t-10 days were used together with the six 
ramped range data to estimate the dynamic constants of 


•siTT is the Simultaneous Interference Tracking Technique which, as 
has been demonstrated, can estimate the relative frequency offset good 
to a few parts of a mHz. 


the Jupiter system. It was hoped that the common error 
sources would be implicitly removed from the simulta- 
neous data. Those estimated values for the mass of Jupiter, 
J2, and J4 and the masses of the four Galilean satellites 
have go^ agreement with earlier determinations (Ref. 4) 
and are shown in Table 3. 

ill. Solar Conjunctions of Pioneer 10-11 

The two Pioneer 10-11 spacecraft had solar conjunc- 
tions in late March and early April 1975. Since the 
celestial longitude is approximately the same for both 
spacecraft, the minimum Sun-Earth-probe angle occurs 
approximately at the same time — March 24 for Pioneer 11 
and April 4 for Pioneer 10. This minimum angle reaches a 
little less than 2 degrees for both spacecraft. 

A total of 4 passes of QVLBI doppler were obtained in 
late March during the solar conjunction of Pioneer 11. A 
summary of these data is shown in Table 4. Unfortunately, 
one pass of three-way doppler received at DSS 44 was 
bad, and another pass of three-way doppler at DSS 42 on 
March 23, when the 2-degree minimum angle occurred, 
was questionable becau% the three-way data were found 
to be 60 percent more noisy than those of the two-way 
doppler received at the same time at DSS 14. According 
to earlier observations (Ref. 3), it is unlikely that the 60 
percent increase in data noise is due to the gradient in the 
solar corona. A possible explanation for the noise in the 
three-way data is that during that time, the 26-meter 
antenna of DSS 42 was at only 9 dB, which was very 
marginal to receive useful doppler data. Two good passes 
of QVLBI doppler were received on March 29 and 31. On 
March 29, the SEP angle was about 4 degrees and the data 
noise of both two-way and three-way doppler was around 
0.026 Hz (60-second count time). After differencing the 
simultaneous doppler, the data noise became 0.0098 Hz, a 
factor of 3 improvement (Fig. 5). Two days later, the SEP 
angle was about 5 degrees, and the data noise was 
improved from 0.0177 Hz to 0.0089 Hz after differencing 
(Fig. 6). 

During the Pioneer 10 solar conjimction, two passes of 
QVLBI doppler were received on March 22 and March 
30. On March 22, the SEP angle was about 10 degrees and 
the data noise of two-way doppler received at DSS 12 was 
0.0134 Hz (see Table 5). There is no noticeable improve- 
ment in the differenced data, which have a data noise of 
0.0132 Hz. This is about the normal data noise of two-way 
doppler at 60-second count time. Comparison with the 
two-way doppler data noise (0.0144) during the Pioneer 11 
conjunction, when the SEP angle was 15 degrees, seems to 
indicate that solar corona effects become less important 
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when the SEP angle u greater than 15 degrees. The effect 
due to the ionosphere was not calibrated because, near 
local noon, it is estimated to be less than 0.002 Hz. The 
most interesting pass in this analysis is the one received on 
March 30 with two-way at DSS 14 and three-way at DSS 
43. The SEP angle was at 4 degrees, and the noise seen in 
two-way doppler residuals was an order of magnitude 
higher than that found during Pioneer 11 conjunction (SEP 
4 deg). This increase strongly suggests active solar 
activities, particularly during the first 15 minutes of this 
one-hour pass, as shown in Fig. 7. After differencing, the 
noise is still as large as 0.28 Hz in the first 15 minutes and 

0.092 Hz in the remaining 35 minutes (Table 5). The noise 
in the differenced doppler implies large spatial gradient 
and/or smaller-scale (less than the DSS 14-^3 distance of 
7000 km) inhomogeneities in the solar plasma. This 
particular pass of data may be useful in studying solar 
activities. 

It is also interesting to see the integrated doppler 
residuals due to solar plasma. Figure 8 shows the 
equivalent range residuals after integrating the two-way 
and differenced doppler residuals as shown in Fig. 7. It 
clearly indicates that the 30-meter variation in integrated 
two-way doppler residuals was improved to less than 5 
meters after differencing. The linear drift in the differ- 
enced range residuals seems to reveal the frequency o&et 
between the two station standards. The factor of 3 
improvement in doppler quality and the improvement in 


range quality of nearly an order of magnitude during solar 
conjunction have demonstrated the importance of this 
differenced data type. 

A summary of per pass (l<r) data noise is shown against 
SEP angle in Fig. 9. The improvement in the differenced 
doppler data noise is not as good as that of the Mariner 10 
solar conjunction, which had a factor of 5 or better 
improvement. This difference is not well understood. A 
possible explanation is that the Pioneer 10 and 11 
spacecraft were 8 and 5 AU, respectively, away and thus 
the signal-to-noise ratio was low compared to Mariner 10 
data. Cortsequently, this tended to cause cycle slipping in 
the two-way or tlu-ee-way doppler, which would degrade 
the differenced doppler. 

IV. Conclusion 

The results of this study have demonstrated the 
navigational capability based on a short arc of QVLBI 
doppler during Jupiter encounter and the improved data 
quality of this data type during solar conjunction. With 
only two weeks of QVLBI doppler at 3 days before 
encounter, the prediction at Jupiter B-plane was good to 
400 km. This differenced doppler has improved the data 
noise during solar conjtmction by a factor of 2 to 3. During 
a time of active solar activities, the integrated doppler 
(range) noise was reduced by nearly an order of magnitude 
after differencing. 
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Table 1. Data summary during Pioneer II encounter 


Station 

ID 

Data 

type 

Time of earliest point 

Time of latest point 

Total 

points 

DSN-43 

F2— S 

03 NOV 74 

10:10:32.00 

12 DEC 74 

12:55:32.00 

919 

DSN-63 

F2— S 

03 NOV 74 

17:35:32.00 

12 DEC 74 

22:30:32.00 

831 

DSN-14 

F2— S 

06 NOV 74 

02:09:32.00 

13 DEC 74 

00:59:32.00 

588 

DSN-14 

ETR— S 

07 NOV 74 

02:30:30.00 

01 DEC 74 

03:20:30.00 

3 

DSN-43 

ETR— S 

13 NOV 74 

08:55:30.00 

10 DEC 74 

07:05:30.00 

3 

DSN-63 

F3C— S 

15 NOV 74 

14:37:32.00 

12 DEC 74 

12:59:02.00 

71 

DSN-14 

F3C— S 

15 NOV 74 

22:47:32.00 

12 DEC 74 

21:13:32.00 

135 

DSN-43 

F3C— S 

18 NOV 74 

04:45:32.00 

13 DEC 74 

03:38:32.00 

152 


QVLBI: Data 

Weight 






F3C 

0.005 Hz 






ETR 

10 km 






F2 

0.5 Hz 
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Table 2. Estimated and consider parameters and the a priori values 


Name 

DX 

Sigma 

A priori sig. 

New value 

Previous 

Nominal 


Estimated parameters 


X 

-0.20963859 +003 

0.12816 +005 

0.10000 +005 

-0.2616640180689983 + 008 

-0.26166192 + 008 

-0.26166166 +008 

Y 

0.94837902 +003 

0.52173 +005 

0.10000 +005 

0.2884050131638451 +008 

0.28839553 +008 

0.28839064 + 008 

Z 

-0.29388873+004 

0.16898+006 

0.10000 +005 

0.1040944516061885 +008 

0.10412384 + 008 

0.10413440 + 008 

DX 

0.39066160-004 

0.33225-002 

0.10000+001 

0.6129445346992491 +001 

0.61294063 +001 

0.61294066+001 

DY 

-0.22425842 -003 

0.13698-001 

0.10000+001 

-0.6547890364795064 + 001 

-0.65476661+001 

-0.65476163 + 001 

DZ 

0.62536468 -003 

0.43661-001 

0.10000 +001 

-0.2471630467680742 + 001 

-0.24722558 +001 

-0.24723309 + 001 

AlOOl 

0.13401015-010 

0.30045-010 

0.10000-010 

0.1340101533576107 -010 

0.00000000 

0.00000000 

AlOO?. 

-0.62326979-011 

0.23034 -009 

• 0.10000-010 

- 0.6232697922876701 - 01 1 

0.00000000 

0.00000000 

A2001. 

0.10147428 -020 

0.31975-019 

0.31623 -017 

0.1014742772031248 -020 

0.00000000 

0.00000000 

A200!! 

0.48452651-020 

0.15498 -018 

0.31623-017 

0.4845265095899446 -020 

0.00000000 

0.00000000 


Consider parameters 


GM5 

Mass of Jupiter 

0.25000 +004 

0.12671303+009 

J502 

J2 of Jupiter 

0.44159-003 

0.14730637-001 

J504 

J4 of Jupiter 

0.60000-003 

-0.59322944-003 

501 CM 

) 

0.75000+002 

0.59369039+004 

502CM 

f Masses of the four satellites 

0.75000+002 

0.32476577 +004 

503GM 


0.10000+003 

0.98768103 +004 

504GM 

) 

0.21500+003 

0.71323529 +004 

L0141 


0.10000-003 

0.24311052 +003 

L043> 

Station longitude 

0.10000-003 

0.14898130+003 

Loeo) 


0.10000-003 

0.35575202 +003 

CU14 1 


0.10000-001 

0.52039962 +004 

CU43 > 

Spin axis distance 

0.10000-001 

0.52052508 +004 

CU63) 


0.10000-001 

0.48624504 +004 

DMW5 

\ 

0.70520-008 

0.44149977 -008 

DPS 

1 

0.53292-008 

0.16706886-006 

DQ5 


0.52240-006 

0.26149561-006 

EDW5 

/ Set 111 ephemerides of Jupiter 

0.26653-008 

-0.58308401-007 

DAS 

\ 

0.28634^007 

-0.15887297 -007 

DE5 

) 

0.26127-008 

1 

-0.24894730-008 



Table 3. Summary of estimated values of the dynamic constants of the Jupiter system 



Sampson 

(1921) 

de Sitter 
(1931) 

Pioneer 10, 
Anderson, et al. 

Pioneer 11, simultaneous 
2-way and 3-way doppler 

GM, 



126,713,600 ±2500 

126,714,027 ± 4000 




(1.4720 ± 0.0040) X 10-2 

(1.4754 ± 0.0022) X IO -2 

h 



(-6.5 ±1.5) X 10-* 

(-5.97 ±3.66) X 10-* 


4.497 

3.81 ± 0.40 

4.696 ± 0.06 

4.606 ±0.03 

^^Europfc 

2.530 

2.48 It 0.07 

2.565 ±0.06 

2.519 ±0.13 


7.988 

8.17 t: 0.13 

7.845 ±0.08 

7.878 ±0.48 


4.504 

5.09 ±0.53 

5.603 ±0.17 

5.637 ±0.11 

^The satellite masses are in 

units of the mass of Jupiter X 10"’. 




Table 4. Data summary of Pioneer 11 solar conjunction 




3-way doppler 


2-way doppler 

Differenced doppler 

Date 

DSS 

No. of data 

Data noise, Hz 

DSS 

No. of data 

Data noise, Hz 

No. of data 

Data noise, Hz 

3/1 

44 

24 

Bad 

12 

28 

0.014 

0 


3/23 

42 

149 

0.16 

14 

145 

0.10 

145 

0.16 

3/29 

11 

14 

0.025 

61 

106 

0.026 

14 

0.0098 

3/31 

11 

58 

0.017 

61 

56 


56 

0.0089 


Table 5. Data summary of Pioneer 10 solar conjunction 




2-way doppler 

Difierenced doppler 

Date 

DSS 

No. of 
data 

Data 

noise, 

Hz 

DSS 

No. of 
data 

Data 

noise, 

Hz 

3/22 

12 

36 

0.0134 

44 

36 

0.0132 

3/30 


13 

0.48 

43 

13 

0.286 

i 14 

35 

0.212 

43 

35 

0.092 
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RELATIVE FREQUENCY OFFSET, mHi 



-7 -i -5 -4 -3 -2 -I ALL DATA 


DAYS FROM JUPITER ENCOUNTER 

Fig. 3. Estimated valM«s of frequency offset 
(all relative to DSS 43) 
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RESIDUALS, Hi X I 0*' _ AB • T, km AB • 




TIME, fflin 

Fig. 5. Residuals of 2-way and differenced doppler on 
March 29 (Pioneer 11), when SEP = 4 deg 
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Fig. 6. Residuals of 2-way and differenced doppler on 
March 31 (Pioneer 11), when SEP = 5 deg 



Fig. 7. Residuals of 2-way and differanead doppler during 
Pioneer 10 solar ooniunction (SEP = 4 deg) 
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TIME, min 


Rg. 8. Integrated doppler residuals during Pioneer 10 solar conjunction 
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Threshold Analysis for VLBI Delay and Doppler 

J. W. Layland 

Communications Systems Research Section 


This article considers the problem of estimating the delay difference and doppler 
frequency difference for the arrival of a white noise signal from a distant Quasar 
at two widely separated receiving stations on Earth. The nonlinear Barankin 
bound is reviewed as it applies to the very long baseline interferometry (VLBI) 
problem, and used to evaluate the signal-to-noise ratio threshold for VLBI Esti- 
mates. We conclude by comparing the performance of several sampling strategies 
for VLBI data. 


I. Introduction 

This article considers the problem of estimating the 
delay difference and doppler frequency difference for the 
arrival of a white noise signal from a distant Quasar at 
two widely separated receiving stations on Earth- At both 
receiving stations, the signal is corrupted by additive 
white Gaussian noise, and further, by filtering, and per- 
haps other processing. This additive noise is statistically 
independent of the signal, and the noises at the two 
stations are independent of each other. Our parameter 
estimator is that parameter transformation which maxi- 
mizes the cross correlation between the corrupted signals 
received at the two stations. If the signal-to-noise ratio 
(SNR) in this estimator is sufficiently high, the maximiz- 
ing parameters are well-defined, and the calculation of 
estimator variance is quasi-linear. As the estimator SNR 


decreases, the problem becomes distinctly nonlinear, and 
bounds or approximations are required. The SNR at 
which the nonlinear effects become important is called 
the threshold of the system. 

The bounds derived by Barankin (Ref. 1) for estimation 
error under a wide variety of conditions apply here. 
These bounds include the quasi-linear Cramer-Rao bound 
as a special case. Swerling (Ref. 2) has applied the 
Barankin bounds to the variance of the estimate of time 
delay and doppler shift for a radar signal, and McAulay 
and Seidman (Ref. 3) have applied it to a threshold 
analysis for pulse-position modulation. The VLBI prob- 
lem can be transformed in a systematic way to the 
problem of estimating the parameters of a “known” wide- 
band Gaussian signal in wideband Gaussian noise, so that 
much of this above-cited work can be directly applied. 
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II. The Barankin Bound 

For the approximate application of the Barankin boimd 
to the VLBI problran, it is convenient to follow the nota- 
tion and arguments of Swerling (Ref. 2). A more rigorous 
and exact formulation could be built from the formula- 
tions used by Barankin (Ref. 1) or Kiefer (Ref. 4). 

Let (F(t, i)} be a family of real-valued functions of 
time t and defining parameter vector We wish to esti- 
mate some of the components of | by observing functions 
Z(t) over Ti<t <Ti where 

Z(t) = F(f,f) + x(t) (1) 

and where x(t) is a real- valued Gaussian random process 
with zero-mean and covariance function tpi(s,t). Given 
appropriate conditions, x(t) can be expanded in the 
Karhunen-Loueve Series 

(2) 

•=» vx» 

where the Xi are independent, zero-mean, unit-variance 
Gaussian random variables, and x>> are the eigen- 
values and orthonormal eigenfunctions of the integral 
equations 


G(f,^"||o) = exp|]£x< ~ ~ 

( 7 ) 

This bound may be made arbitrarily tight by appropriate 
choices of K and (ot. ^)- 

The Cramer-Rao bound can be derived from Eq. (6) by 
taking a limit as two of the smoothly approach each 
other at ^o, and their corresponding o* grow equally and 
oppositely to infinity. The resulting zeroth-index deriva- 
tive term can be included in the bound (6) as (Ref. 3) 

(8) 

i:sO IsO 

where 

doO=-^g|7G(f'.rU.)|{- = f-={o 

doi = d*o = -^G(|',|i' lfo)|f ={o te[i,X] 

~ G(^,|"|fo)|*.ieu,jt) (9) 


4>(s) = xj^'*>is,t)4>(t)dt (3) 


We further suppose that for all we can express 




over the observation interval, where 


= rF(t,^)Ut)dt ( 5 ) 


Barankin’s result for the greatest lower boimd to the 
variance of any estimator for ( is 


— ( 6 ) 

ZZG(^i;,^;'|fo)-a*a. 

lal 


Equation 8 is the Cramer-Rao Bound if Oi = 0, for all 
» > 0. Continuing to follow Ref. 3, let A denote the 
column vector of the a,, and N denote the column vector 
of the tii defined by n» = 1, nj = {( — fo for »e[l,K]. If 
D denotes the matrix of {d*i|k,l€0, K}, then the bound 
may be rewritten as 

^ > A'NN'A/A'DA (10) 

The right-hand-side of Eq. (10) is maximized by letting 
A = X.D~'N, for any A, so that the greatest lower bound 
on the variance is 

<7* > A/‘D->N (11) 

We will apply this form of the Barankin Bound to the 
VLBI estimation problem. 

III. Very Long Baseline Interferometry 
Detection and Parameter Estimation 


where is the true r>arameter value, K is any finite 
integer, the K are alternate parameter values, and the a* 
are arbitrary real numbers. In the above notation, the 
fimction G(*, • | •) is defined by 


X sliov/s MM ovs3r\”0\v of tlio features of the VLBI 
problem. The white noise signal from the radio star is 
received at the two stations with some relative delay and 
doppler offset which can be related to Earth’s rotational 
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position and velocity. Local ambient noise, usually much 
stronger than the radio-star signal, is added at each re- 
ceiver, and the resultant is filtered, translated by a local 
oscillator to a manageable frequency and bandwidth, 
and then sampled, hard-clipped, and recorded. The re- 
corded signals are then brought together at leisure, modi- 
fied by an a priori estimate of the relative delay and 
doppler, and cross-correlated. Refined estimates of the 
delay and doppler are e.xtracted from this cross-correlation. 

The input to the cross-correlator from station 1 
may be denoted as = Sign {S(; + Ik,/) + n,(; k,f)}, 

and the input from station 2 may be denoted 
tjj = Sign {S*(; -f d,f„) n-(;,/„)}. Here S*(*,*) Is the 

naturally delayed and doppler-shifted version of S(*, •), 
and k and / are our a priori estimates of this delay and 
doppler, respectively. Because of their origins as white 
noise, and subsequent filtering, S(*, •) and n,(*, •) have 
identical autocorrelation functions, as do S*(*, •), and 
*) The cross-correlation 


where N is the number of samples, and ^(*, •) is the 
normalized autocorrelation (ambiguity) function for the 
filtered signal process. 

The second moment is 

Z(fc,/)- Z(i,g)= 

S S [’’•♦I + ) 

' ‘ (14) 

If we expand the right hand side of Eq. (14), and drop 
those terms with zero expected values, 

Z{Jc, /) • Z(l, g) = L E {rM • »•»« •r’rr'i 

i • 

■f ■ 

+ ri.k*ri»i • m'i • m'i 

+ (15) 


Z(fc,/)=2:.r;.,-i/, (12) 

j 

is evaluated as a function of k and / near its maximum. 
The final estimate of the actual delay and doppler is made 
much finer than the integral steps in k. 

Let us denote by T,, and the noise temperatures of 
the environment at stations 1 and 2, respectively, and by 
T,, and T.j the increment to these noise temperatures 
caused by the radio-star “signal” at the two stations. T,, 
and r.- may differ as a result of differing antenna collect- 
ing areas, efficiencies, etc. Denote as Tj and nij the signal 
and noise terms respectively in X/, and as and m'j the 
signal and noise terms in i/,. By definition of the hard 
limiting, the average values of (r, -I- mj)^ = (x,)’ = 1 and 
(r'i + m'iY = (i/j)= = 1. Furthermore, Er{{E„(rf + m,})=) = 
(2/r) (T../r„.), and E,((£„{r' + m;})=} = (2A) (T^/T^), 
where Ea denotes expected value over the distribution of 
the random variable a, and where we assume T,i«Tni- 


In order to get the most information possible on a per- 
sample basis, the sampling is performed at a rate such 
that X)Xt s: 0 and ~ 0 for ; ^ k. Hence at die sample 
points 

m,^k • mui = P - 

m’j 'm\ — (16) 

For some constants p, p' and any particular f, /, k, I. We 
will assume that Eq. (16) holds at interpolated values 
also. With this assumption, Eq. (15) may be rewritten in a 
statistically-equivalent but simpler form; 

Z{k, f) Z(l, g) = £ 2: • r,., . rj • K 

(17) 


Because the number of samples in the summation (12) 
is large, Z(*, •) will be approximately Gaussian (by the 
central limit theorem), and its distribution specified by its 
first two moments. The first moment is 

i 

= Z 

i 

= • <^{k - d,T- /„) (13) 

7T V i nl 1 « 


Equations (13) and (17) are representative of the detec- 
tion of a known Gaussian signal, (rt) embedded in 
Gaussian noise. The constants p, and P' are 



On a per-sample basis, the noise-to-signal ratio for this 
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detection problem is 


where 


NSRj = /S + )9' + 



It is convenient to consider the detection problem as if 
the interfering noise samples have unit variance. If it is so 
normalized, each of the “signal” samples has energy given 
by E, = 1/NSRj, or 


R = N‘ 

The D matrix of Eqs. (9) to (11) is defined from Eq. (23) as 
D = {d«|l:.l=0. -.K} 

d*.i = exp (R[^(0) - <!>(& - lo) - <^(f" - lo) + m 




With the VLBI parameter estimation problem thus 
transformed to one of parameter estimation for a known 
signal in noise, the formulation of the Barankin bound 
(Eq. 8) is directly applicable. Since the sequence of sam- 
ples of the signal and noise processes at both stations are 
virtually independent from sample-to-sample, we take as 
the orthonormal basis for the signal set, the sequence of 
sample functions themselves. We normalize so that the 
interfering noise samples have unit variance, and hence 
so XI = 1. fo*’ sample indices j. The signal projections 
are a scaled version of the tj. The number of sam- 
ples, N from which the estimates of delay and doppler 
are made, is large, so that central limit theorem applies, 
and 

( 21 ) 

j 

Or using the normalized autocorrelation function 

i: xMi) • /3 Af') N Hi - i') 

T«-T« 

(22) 

Note that here ( denotes a vector parameter which can 
include delay or doppler offsets, or both. 

The G(*, * I •) function for VLBI (Eq. 7) is defined by 

G(f'. ni*) = exp (R- [0(1' - n - Hr - r) 

-Hr'-r) + Hr-rm (S3) 


for k, l€[l, K] indices of alternate f. 

d*.« = d„.t = Ji'^iHr - r:) - nr - fo)|f-=fo 

forkelLK] 

d<„, = R.-|r^*(f'-r)|j=t-=f„ (24) 

Furthermore, we note that is presumed to be the global 
maximum of 0(*) so that the second term of dot is trivially 
equal to zero. If we restrict the set of alternate to local 
maxima of 0(*), the first term of dot, and hence dot itself, 
is zero for all k. With this restriction on D, the Barankin 
bound for a* becomes 

a>-^ + N[Di^N, (25) 

“o.o 

where the subscripts 1 in the second term of Eq. (25) are 
meant to imply that the zero-index row and column have 
been dropped. The first term of Eq. (25) is the Cramer- 
Rao Bound; the second is the nonlinear threshold term. 

IV. A Very Long Baseline Interferometry 
Example 

The wideband quasar noise signal is processed at the 
two receiving stations into a form which is suitable for 
recording. This processing corresponds to constraints 
on available equipment bandwidths, recording rates, data 
storage capacities, etc. The basic constraint is the band- 
width of the low-noise Maser amplifier, which we will 
hypothesize to be 40 MHz centered at 2 or 8 GHz for 
present systems. This is still too wide for any but the 
fastest recorders, and is further reduced by filtering sev- 
eral narrow channels from the passband which span the 
full bandwidth available (bandwidth synthesis); or by 
samphng the fiJl bandwidth into a buffer memory, and 
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then recording at a slower rate once that buffer is filled 
(burst mode). Both sampling modes have side-lobes in the 
delay-estimator which arise at roughly the reciprocal of 
the spanned bandwidth (1(H sec). There are also side- 
lobes in the frequency estimator at roughly the reciprocal 
of the recording or sampling interval (Kh’ Hz). Burst 
sampling will generate additional sidelobes spaced a few 
Hz away from the true frequency, but these should be 
able to be removed by a priori knowledge, as should some 
of the more widely spaced side lobes in time-delay. 

Bandwidth synthesis systems which have been pro- 
posed utilize two channels at the extremes of the avail- 
able bandwidth, plus a third channel placed somewhere 
between to reduce the sidelobe amplitudes. The exact 
sidelobe structure is sensitive both to the position of this 
third channel, and to the bandwidth of the channels. For 
an overall spanned bandwidth of 40 MHz, a reasonably 
tolerable side-lobe structure is achieved if the third chan- 
nel is located 5 MHz from the passband center. While 
this position ,is perhaps not optimum, it offers near 
optimum side-lobe suppression which is far better than 
that provided at many alternative positions, such as at, 
e.g., 6 MHz. Figure 2 shows the effect of varying the 
individual channel widths for a 40-MHz spanned band- 
width with the third channel located at 5 MHz. The 
function plotted is the aggregate autocorrelation function 
for the three recorded channels, which results with opti- 
mum use of all available information. The recording 
bandwidths are 6 MHz, 4 MHz and 2 MHz for Figs. 2a, 
2b, and 2c, respectively. The corresponding Barankin 
lower bound to delay estimate error is shown in Fig. 3, 
lines B, C, and D, respectively. As could be anticipated 
from the sidelobe structure of Fig. 2, the threshold per- 
formance degrades seriously as the individual channel 
bandwidth is narrowed, and would continue to degrade 
if the recording bandwidth were reduced further. At 
Estimator signal-to-noise ratios well above threshold, the 
delay estimate error depends only minimally upon the 
individual channel bandwidths. 

Figure 3 also shows the Barankin lower bound to 
the dday estimate error achieved by a full bandwidth 
40-MHz Burst Mode System (line A). The autocorrelation 
function is the familiar sin (x)/x. Threshold performance 
of the burst mode sampling is better by about 2 to 3 dB 


than that for bandwidth synthesis with 6-MHz charmels, 
and better by about 10 dB than that for bandwidth 
synthesis with 2-MHz charmels. At high SNRs, however, 
the delay estimate error is 2 dB worse for burst mode 
than for bandwidth synthesis. 

The doppler estimator error has not been calculated 
per se, but its shape can be inferred as follows: The 
recording of data is approximately uniform in time over 
some fixed interval. Accordingly the transform of this 
interval is the sin (x)/x function and the doppler estimate 
error is identical to line A of Fig. 3 with the abscissa 
relabeled as appropriate: for example, if the recorded 
sample occupied 40 sec, then ot = 10"® sec on Fig. 3, 
line A is comparable to a/ = 10^ Hz. 


V. System Implications 

One can draw a variety of conclusions from Figure 3, 
depending upon where and how he looks. Above thresh- 
old, for example, a three-fold improvement in delay 
estimate error requires roughly a 10-dB increase in esti- 
mator SNR, or a hundred-fold increase in the number of 
samples processed. This is a potentially expensive way to 
gain accuracy. This improvement in accuracy could per- 
haps be better achieved by other means, such as increas- 
ing the spanned bandwidth, if possible, and operating the 
system at some nominal margin above its threshold. We 
would like furthermore to be able to operate with a per- 
formance curve like Fig. 3A or Fig. 3B with a threshold 
as low as possible to minimize the amount of data needed 
to ensure above-threshold operation. Line 3B corresponds 
to a recording channel bandwidth of 15 percent of the 
spanned bandwidth, and for spanned bandwidths which 
are significantly larger than the 40 MHz considered here, 
it suggests that the use of a combined buist/bandwidth 
synthesis recording is desirable. For example, it does not 
appear reasonable to burst-sample a 200-MHz spanned 
bandwidth with present-day equipment, but 15 percent 
bandwidth channels, each 30 MHz wide, could be easily 
burst-sampled and recorded. Such operation would lead 
to a performance curve of the type of line B, Fig. 3, but 
with the delay-estimate error in the above-threshold 
region reduced five-fold by the increased spasmed band- 
width. 
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Fig. 3. Delay estimate error vs. estimator SNR in dB 
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N76- 18174 


Analysis of Command Detector Signal-to-Noise 

Ratio Estimator 

R. G. Lipes 

Communications Systems Research Section 


We have investigated a specific technique jot making (SNR)''‘ estimates by 
using in-phase channel output averages to estimate signal and quadrature channel 
output averages to estimate noise. We have produced bounds to determine the 
accuracy of this technique when fluctuations of one standard deviation occur. Our 
results show the estimate is relatively independent of actual input signal-to-noise 
ratio (SNR) and can be improved only by irwreasing the number of samples in the 
averages. 


I. Introduction 

The present design of a command detector being de- 
veloped for NASA calls for a signal-to-noise ratio (SNR) 
estimate that will be used for monitoring operations of 
the detector. Naturally, we wish to know for a given 
SNR and ntimber of samples how good the estimate will 
be. The samples are derived from in-phase and quadra- 
ture outputs that are efiPectively data and error inte- 
grators. As is to be expected, the command detector 
circuitry must be as simple as possible. For this reason 
the absolute value of the in-phase and quadrature out- 
puts rather than their squares will be sampled. Further- 
more, since samples of the absolute value of the output 
provide an estimate of the square root of the SNR rather 
than the SNR itself, in this report we will be concerned 
with the relationship between SNR, number of samples. 


and the estimated square root of SNR. We will use the 
ratio of the mean plus standard deviation of the estimate 
to actual (SNR)’'* as a measure of the accuracy of the 
estimate. This ratio will be determined as a function of 
number of samples. 

In Section II, we will model the probability distribution 
of the in-phase and quadrature outputs and develop an 
estimate for the square root of SNR. In Section III we 
will give bounds for the mean and variance of this 
estimate that converge to the exact v^ues as the number 
of samples becomes large. Because these asymptotically 
tight bounds are easily calculated, they are presented 
rather than the exact values whose inte^al representa- 
tions required numerical integration for evaluation. Fin- 
aUy, in Section IV we will discuss results and conclusions. 
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II. Probability Distribution of Samples and 
Development of Estimate 

A. Data Integrator Output 

The output of the in-phase or data channel is inte- 
grated for one symbol time T. The absolute value of the 
resulting integration becomes a single sample so that 
the tth sample Xj can be expressed as 


Xi 



[d(t) + n,(t)] 


dt 


( 1 ) 


where n^(f) is the noise of the quadrature channel also 
asstuned to be white Gaussian with power spectral den- 
sity No/2. Because the samples Yj are taken from non- 
overlapping intervals, they are independent and identi- 
cally distributed. Their probabihty density is obtained 
from Eq. (4) with AT = 0: 

(o. 


o > 0 

( 6 ) 

a <0 


where d(t) is the data assumed to be ±A over one symbol 
time, and n,{t) is the noise of the in-phase channel as- 
siuned to be white Gaussian. Thus, 


C. Development of Estimate 

The output of the in-phase channel before the absolute 
value is taken is 


Xi = 


zAT -I- 


fti*T 

/ n,{t) dt 

Jh 


The random variable 



n,(t) dt 


is zero mean Gaussian with variance 


( 2 ) 



-I- ni{t)]dt = zizAT + Z;,i 


(7) 


The square of the mean of this output is (AT)* while its 
variance (Eq. 3) is u* = NoT/2. Conventionally, the in- 
put SNR is defined as the bit signal energy divided by 
the one-sided noise spectral density or 


/ 1<.+ T rii*T V x 

dt ds n,(t) n,(s) = — g— (3) 


A»T 1 (ATY 
No ~ 2 (N„T/2) 


where No/2 is the power spectral density of n,(t). Since 
the noise is assumed white Gaussian and the samples Xj 
are taken from non-overlapping intervals, they are inde- 
pendent and identically distributed. From Eqs. (2) and 
(3) we have for the probability density 


Thus the signal-to-noise ratio estimator of the command 
detector will estimate the quantity 


r (ATY 1 


AT 1 

a 2*'* 


(8) 


(W)’« 

0, a < 0 

(4) 

where = NoT/2. 

B. Error Integrator Output 

The output of the quadrature or error channel is also 
integrated for one symbol time and the absolute value 
taken to form a single sample Y{: 


In Appendix A we show that the mean of the sample 
Xi (Eq. 1) is given by 

Xi =At[i+2€^^)] (9) 

where the function f.(AT/ti) is defined by Eq. (A-4) of 
Appendix A. This function is sufficiently complicated so 
that forming an unbiased estimate of AT using only the 
{Xi} is not feasible. Nevertheless, for input SNRs greater 
than 3 dB (the design point input SNR for the command 
detector is 10.5 dB), |e|< 0.0042. Consequently we can 
neglect e and employ M samples of the (Xi) to form our 
estimate v of AT : 



Yi = 


rti*T 

i ’ 


n«(t) dt 


(5) 



( 10 ) 
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From Appendix A by setting AT = 0, we find the mean where we note that {Xi) and {Yi} are statistically inde- 
of Yii pendent. 




Thus, an unbiased estimate of <r would be 


III. Expressions for the Mean and 
Variance of Estimate 

The mean and variance of our estimate W are given by 



However, <r is in the denominator of Eq. (8), so A*' will 
not proxade an unbiased estimate of <r~‘. As we show in 
Appendix B as A/— » oo, — » <r'* so we will use A** as 
an estimate of even though it is not unbiased. Our 
estimate W of (SNR)*/’ is, therefore. 


W = TA-‘-^ (13a) 

2<rw = A"* — (v)* (A"‘)* = <r?<ii-i + <ri(A~*)’ + (v)*«ri-i 

(13b) 

The mean of v is and is given by Eq. (9), while the 
variance is 1/M times the variance of Xt; so from Ap- 
pendix A: 




( 12 ) 


(14) 


In Appendix B we derive integral representations for the mean and variance of A~*, but here we will display only easily 
calculated bounds for these quantities that are obtained from the integral representations: 


1 M 1 4M ^ ^ 1 M , 1 4M 2 

<rM + l'^a(M-l)(M + 3)\^*y (Af- l)(M + 3) ■» 


(15a) 


, ^ 1 ( 64M (M + 1) 1 , 8Af 1 M’ 

- cr* \(M« - 4)(M 4- 4) ^ (M -I- 1)(M -f 4) V (M -I- 1)’(M + 2) , 

1 ( 8M> P 4M IV ) 

(r’ t(M*-l)(M + 3) Vir^ ■'■[(M-IHM+S) J f 


where the derivation requires M > 3. Since (1/2) (AT/v)^ is the actual SNR, we have the following bounds on our 
estimate W using Eqs. (9), (13), (14), and (15), 


[1+ 2.K W < (SNR)«[ + (M-ltrM + 3) (l)][* ^ 

( 16 a) 

. 64M(M + 1) 1 , 8Af 1 M’ 8M’ 

‘'*'-t(M’-4)(M + 4)»’'''(M-(-l)(M+4) * (M + 1)’ (M + 2) (M> - 1)(M + 3)[ir ) 

[(M - 1) (M + 3)] ( 7 ) } { 2 M ~ f ^ 

[ mTT (M-l)(M + 3) (7)] ~ ® 2M 


(SNR)’/’ 


M 


4M 


M-f 1 (M-l)(M + 3) 




where e = «([2(SNR)]’/’) is tire function defined by Eq. (A-4) of Appendix A. 
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IV. Discussion of Results and Conclusions 


We can measure tlie accuracy of the estimate by con- 
sidering the ratio 


-f aw 
l§NRp" 


(17) 


where is the upper bound on the mean W of Eq. (16a) 

and aw is the upper bound on the standard deviation of 
Eq. (16b). This ratio should be close to unity if our 
estimate is good and measures roughly how closely the 
estimate approximates (SNR)*^’ when statistical fluctua- 
tions of one standard deviation occur. We have plotted 
this ratio in Fig. 1 as a function of M, the number of 
samples, for input SNRs (as defined by Eq. 8) of 0, 5.25, 
and 10.5 dB. The design point input SNR for the com- 
mand detector is 10.5 dB so the 0 dB curve repre- 
sents performance of the estimate when the input SNR 
is 10.5 dB below design. 


From Fig. 1 we notice immediately two features of the 
estimate. First, the accuracy of the estimate is quite 
insensitive to actual input SNR; the ratios are within 
0.15 dB of each other for A/ > 16 for input SNRs diflfer- 


ing by 10.5 dB. Furthermore, mcreasing input SNR 
above 8 dB has virtually no effect on the estimate. Sec- 
ond, the accuracy is extremely dependent on the number 
of samples: an estimation accurate to within 1.5 dB 
requires about 45 samples while one good to within 
0.5 dB requires more than 500 samples. The number of 
samples M must be increased to improve the estimate, 
and since the standard deviation of the estimate de- 
creases as for M large, an extremely large number 
of samples is required for very accurate estimates. Fig- 
ure 1 gives quantitative support to these concluding 
remarks. 

Finally, we should note that previous analyses (Refs. 
1-7) have arrived at essentially the same conclusions. 
The estimate in this work is for SNR*^* and is obtained 
from averages of absolute values of in-phase and quadra- 
ture channels. In the previous works the estimates were 
for SNR and were obtained from averages of squares of 
the relevant chaimel outputs. Intuitively, we might ex- 
pect this work to agree qualitatively with the other 
analyses, but, in fact, it agrees quantitatively as well. 
For given input SNR and number of samples, averages 
of either absolute values or squares provide estimates of 
approximately the same accuracy (Ref. 7). 
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RATIO f 



Fig. 1. Plots of ratio of upper bound of estimate plus one standard deviation o„. to actual SNR*/2 as a function 
of number of samples M. The three plots are for input SNRs of 0, 5.25, and 103 dB defined by Eq. (8) of the text 
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Appendix A 


In this appendix we will derive the mean and variance 
for a random variable X whose probability density is 
given by 

+ «>0 

0, o < 0 

(A-1) 



Since Q(—x) = 1 — Q(x), we have 

X = (A-3) 

where the function €{y) is defined by 

€(y) = exp [ - 1 y i - Q(y) (AA) 

and vanishes as exp [— i/V2](l/y’) as y becomes large. 


I. Mean of X 

The mean of X is given by 
X = a Px(o) da 

*L. 

<-) 

where we use the function 

grx;=y%xp[-in’]^ 


II. Variance of X 

The second moment of X is 


- /■“ . r l/a-AT\-l da 

X«-|aexp|^ „ ) J(W)»* 

/•“ , r l/o + ATVI da 
+ 1 a=e^[- 2(— 7-j JiWPS <A-5) 

If we substitute /9 = — a in the second integral, we have 


(A-6) 

which is the second moment of a Gaussian random vari- 
able with mean AT and variance a*. The variance of X is 

(A-7) 

where €{AT/a) is given by Eq. (A-4). 
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Appendix B 


In this appendix we wish to give integral representa- 
tions and bounds for the mean and variance of 
where 


where the subscript M indicates Af samples. We will 
utilize the following identity to obtain a representation 
for [ci + ••• + 


A = 



Yi 


I. Integral Representations 

We will give a representation for A~^, where M > L + 1, 
which insures convergence of the integral. From the 
probability density for (Eq. 6 of the text), we have 


=/•••/&. - da^ exp 


X 




(B-1) 


[c -t- Oi + • • • + Ojf ]"* — 

d/8 exp { — /S [c -I- a, -I- • • • -f- a*]) 

(B-2) 

Differentiating L times with respect to c and setting 
c = 0 gives 



[Ci -I- • • • -i- Cm] *'■ = 


1 

(L - 1)1 



exp {- ;8[ai H + oj,]} 

(B-3) 


Substituting (B-3) into (B-1) and interchanging the order of integration gives 


(B.<) 


where 


The mean and variance of are obtained from (B-4) using L = 


dy 

(27t)‘'* 

1 . 2 . 


II. Bounds for Mean and Variance 

The bounds we will obtain depend upon the inequality of Ref. (8), which for our application states: 

(t) <'W exp j < (-iy '* [X + (i* + 8/,r)*/»]-' (B-5) 

Applying this to (B-4) gives 

1 /MX'' /■« 1 

a.- TT ! -(~) ( - ) . [* + (** + dx < A-J- <7^zriM 

/MX*- /2X<*'*‘'*'’ /■* 

X^-^j I i^*[* + (x’ + 8/T-)*'*]-*dx (B-6) 
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A. Bound for the Mean 

Consider (B-6) when L = 2. The integrals involved can 
be done by letting x = osinh (t), so 


X [X + (X* + a»)w*]-^dx = (B-7) 


Thus, (B-6) becomes 




(B-8) 


We can obtain the boimd on A]) by relating A^* to Aj* 
Consider (B-4) for L = 2 and integrate by parts with 
u = (Q(x))“ and dv = exp [ — (l/2)x*M] x dxi 






-I- 


M 2 Af’ 1 2 

— 2 ^ T7 ^ Ai 

a IT 


Using (B-9) with inequality (B-8), we find: 


(B-9) 


1 M ^ 1 4M 


(f) 




ff M + 1 <r (M - 1) (M -f 3) 

1 Af 1 4Af 2 


< a; 


< 


Af -f 1 <r (Af — 1) (Af -f- 3) ir 


(B-10) 


B. Bound for the Variance 


Consider (B-6) when L = 3. Again letting x = a sinh (f) 
simplifies the integrals giving: 


i: 


x» [x + (x* a»)‘/*]-^dx = 


2M 


(Af» - 1) (Af « - 9) o“-« 

(B-11) 


So, (B-6) becomes 

1 8Af‘ 


( 2 

v) 


8Af‘ 


(Af*-l)(Af- 




(B-12) 


We can relate A-J to A^* by using L = 3 in (B-4) and 
integrating by parts with u = x(Q(x))' and do = 
exp [(l/2)x’Af] X dx: 


1 4Af* /2V''*”" — 1 4Af' /2V 




1 Af' — 1 Af — 




(Af - 1)^ 


(B-13) 


Using (B-13) with inequality (B-12) we find a bound for 
the variance <ri-> (Af — 1) = Aj^, — ( A-'_j )' 


oX-.(Af- 1) < 


1 64Af(Af - 1) 
vV' (Af -l-l)(Af»-9) 

, 1 /Af- 1\'_ , 

V ( Af ) 


(B-14) 


Using bound (B-10) for A^', A^.,, we finally obtain 
64Af(Af + 1) 1 




+ 


4) (Af -I- 4) 
8Af 1 


Af' 


(Af -I- 1) (Af -I- 4) «■ (Af-l-l)'(Af -t-2) 
8Af' /2\<»«)» 


r} 


If 8Af' /2\ 
o'\(Af*-l)(Af -t-3)\»j 


(B-15) 
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High-Power Transmitter High-Voltage 
Power Supply Ripple 

E. J. Finnegan 
R. F. Systems Development Section 


This article reports the results of decreasing the high-voltage power supply 
ripple by redesigning the low-pass filter. 


I. Introduction 

During a routine checkout of the high-power transmit- 
ter at DSS 14 located at Goldstone, California, it was 
noticed that the digital fast klystron body detector, a body 
overcurrent sensing device, continuously tripped the 
transmitter high-voltage off. Investigation showed that the 
ripple voltage of the high-voltage power supply had 
increased sufficiently to cause the body peak currents to 
activate the fast body detector. 


II. Investigation 

Measurements were made of the ripple voltage at the 
output of the dc filter (Fig. 1). It was found that the ripple 
voltage was very high, approximately 300 volts; the high- 
voltage specification is .1 percent of the high voltage. This 
high ripple voltage indicated further problems in the filter 
choke (inductor) or capacitor. The capacitor was checked 
and found to be in good condition. The inductor was then 
checked and found to have changed its inductance from 1 
henry to .1 henry. It was theorized that the gap in the iron 
of the inductor had changed resulting in a change in 


inductance. However, the high-voltage insulation isolating 
the inductor coil from the iron core had not broken down. 
It was tested to 80 kV with a high-potential tester. The 
spare inductor was installed and the ripple voltages 
checked again including the inductor input and output. It 
was found that the input ripple was alarmingly high. The 
three phase full wave output of the transformer-rectifier 
should be 4.2 percent instead of the 75 percent found. 
This indicated strongly that the filter system was 
resonating. DSS 63 and DSS 43 were then checked and 
they were found to be operating in the same condition; 
the ripple voltage was 30,000 volts while running the 
klystron at 35,000 volts. The filter in its old configuration 
was a choke (inductor) input filter. It was then changed to 
a pi configuration. Tests were run to detennine the 
optimum value of the input capacitor to be used and .2p.F 
was the determined value. 

III. Conclusion 

A temporary Ejigineering Change Order was written in 
order to change the filter to a pi configuration using a 
.I5/sF capacitor on the input that was immediately 
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available. In the future this will be changed to the .2ftF 
capacitor. Results of the new filter are shown in the 
enclosed photos. Figure 2a shows the ripple voltage 
without any capacity at the input to the filter which is 

15.000 volts (original configuration). The dc voltage is 

20.000 volts. Figure 2b shows the input of the filter at 20 
kV dc with a .15pF capacitor at the input pi configura- 
tion. As can be seen, the ripple voltage is reduced to less 
than 1,000 volts peak-to-peak. 

Figure 3a shows the ripple voltage at the output of the 
filter (choke input). The ripple voltage is greater than 50 
volts peak-to-peak at 20 kV dc. Figure 3b shows the ripple 


voltage at the output of the new filter, pi configuration. As 
can te seen, the ripple is reduced to less than 20 volts 
peak-to-peak. 

Figure 4 shows the ripple at the klystron operating level 
of 60 kV. The ripple at this point is approximately 30 
volts peak-to-peak. The input ripple to the filter was 3,800 
volts; this is within specs for the filter and the output 
ripple is within specs for the transmitter. 

This will reduce the modulation on the klystron output 
in the future, thus transmitting a purer spectrum. 
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Fig. 1. High-voltage power supply and filter configuration 
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Venus Station Automation: 
Communications Link 

A. Zygielbaum 

Communications Systems Research Section 


A highly error-resistant communications system has been developed in 
support of the demonstration of the remote automatic control of the Venus 
Station (DSS 13). This article describes in general the computer communications 
software used to connect the JPL control terminal to the Venus Station, and 
describes in some detail the software used to drive the teletype line which is part 
of the link. 


I. Introduction 

Remote automatic control of the Venus Station (DSS 
13) has recently been demonstrated (Ref. 1). During this 
demonstration, data from several pulsars were acquired by 
a Sigma 5 computer at ]PL. The system configuration at 
the Venus Station is an outgrowth of that used in a 
demonstration of local automatic control of that station 
(Ref. 2). The system conjuration at JPL is similar to that 
used previously to monitor the performance of the Mu-I 
ranging system at DSS 14 (Ref. 3). This article describes in 
gmeral the computer communications software used to 
conirect the JPL control terminal and Venus Station, and 
describes in some detail the hi^y error-resistant software 
used to drive the teletype (TTY) line which is part of the 
link. 

At the Venus Station, an XDS 930 is used as the main 
control computer. It is connected to two XDS 910 
computers by high-speed links known as D/R units (Rek. 


2 and 4). Currently, the only available data path from the 
Venus XDS 930 to JPL is via a 10 character-per-second 
teletype line. Using standard XDS teletype buffers, the 
930 is connected to an XDS 920 at JPL, which, in turn, is 
connected to a Sigma 5 computer via another D/R link. 
The software which drives t^ link has been previously 
described (Ref. 5). 


II. Definition of Task 

The primary task was to implement a simple to use, 
reliable coirununications link between the Venus XDS 930 
computer and a Xerox Sigma 5 computer at JPL. The 
topology of the path between these machines in shown in 
Fig. 1. To facilitate user programming, identical Fortran- 
callable user interfaces were defined for both ends of the 
link. Extensive error checking and retransmission capabil- 
ity were required for reliability. 
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Three software design and development efforts were 
embodied in the primary task. A teletyp>e handler was 
needed to link the XDS 930 at the Venus Station and the 
XDS 920 at JPL. To extend the connection, a “switcher” 
program had to be written for the XDS 920 to interface 
the teletype handler and the D/R software controlling the 
link to the Sigma 5. Finally, Fortran-callable subroutines 
were required for the Sigma 5 to drive the D/R software, 
and provide a user interface identical to that in the Venus 
XDS 930. 


III. User Interface 

The teletype handler appears as the collection of 
Fortran-callable subroutines shown in Table 1. Identical 
software is used at each end of the link. As a result, both 
the Venus XDS 930 user and the XDS 920 “switcher" 
utilize the same subroutine calls to access the teletype 
lines. These same calls are made available in the Sigma 5 
by the transparent combination of the XDS 920 switcher 
and the Sigma 5 interface package. 

Two logical full duplex channels are provided by the 
communications link. Messages conveyed by the link are 
blocked to facilitate retransmission. The user can assign a 
transmission priority of 1 to 15 (1 is highest) to either 
channel. The channel with the highest priority can 
preempt the other channel. If they are given equal 
priorities, data blocks on the channels are interleaved so 
that the traffic on both channels appears to be transmitted 
simultaneously. 

Sending a message is equivalent to transferring a core 
image from one end of the link to the other. The 
respective users must supply the buffers from which these 
data are taken and into which these data are placed. At 
either end. two transmit flags control transmission on both 
channels while two receive flags control reception. 

To supply a buffer and begin transmission, the 
SNDMSG(LOC,CHN,PRY.RTRY,CHCNT) subroutine is 
used. The parameters are: IXX), the first location of the 
buffer; CHN, the channel munber (1 or 2); PRY. the 
channel priority; RTRY, the number of times any 
particular block of the message is retransmitted before a 
line error is declared; and CHCNT, the number of 4-bit 
bytes in the buffer. A buffer is designated for reception 
with the SETRX (CHN.LOC.CNT) subroutine. Once 
again, LOG and CHN are the first location in the buffer 
and the channel number, respectively. The size of the 
buffer is given as CNT computer words to prevent data 
overflows. 


It would be luiacceptable to transmit a message if a 
receive buffer had not been provided (i.e., if SETRX had 
not been executed at the receiver). Coordination of the 
two sides of the link for transmission is provided by 
transmit flags. The mechanism for coordinating flags at 
both ends of the link will be discussed later. 

These flags are accessed by the subroutine TXSTAT 
(CHN.FLG), where CHN is the channel number and FLG 
is a dummy variable into which the flag value is placed. 
The transmit flag normally has the value 1, indicating that 
a receive buffer for that channel is available at the remote 
end of the link and that the channel is free. A value of 0 
implies that SNDMSC has been executed and, therefore, 
the channel is busy. If the transmission is successful, the 
flag remains busy, 0, until a new receive buffer is supplied 
via SETRX. When the receiving user supplies the buffer, 
the transmit flag again becomes 1. If the transmission is 
unsuccessful, that is, if a particular block has been 
retransmitted RTRY times, then the flag is set to -1. No 
action is required by the receiver user as the receiving 
buffer is still considered empty. 

Data reception is coordinated through receive flags. A 
receive flag is accessed by the subroutine RXSTAT 
(CHN.FLG). Again, CHN is the channel number and FLG 
becomes equal to the receive flag. Once a receive buffer 
has been supplied, the receive flag remains 1 until a 
message has been succesfully received. Note that there is 
no indication to the receiving user that a message receipt 
is in progress or that an unsuccessful transmission was 
attempted. If a correct message is received, the receive 
flag is set to 0. It will not be reset to 1 until SETRX is 
called to supply a new buffer. 

A utility subroutine MSGCNT(CHN,CNT) is provided 
to indicate the number of 4-bit bytes received in a 
message. This subroutine will return the count as the 
dummy variable GNT for channel CHN provided that the 
receive flag is 0, that is, if a message has just come in. If 
the flag has any other state, MSGCNT will return -1. 

Sample receiver and transmitter programs are shown in 
Fig. 2. The transmit routine first executes a busy wait until 
the channel is free. Upon exiting this loop, it starts the 
transmission of the buffer IRRAY. Coincidentally, the 
receiver waits for a message to come in. Upon receipt, 
MSGCNT is used to determine the incoming byte count. 
Elxecution of SETRX supplies a new receive buffer and, 
consequently, frees the channel. 

System initialization is facilitated through the TINIT 
subroutine. This subroutine may be called by either user 
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at any time. Its e£Fect is to initialize all flags and to 
terminate all trafiic activities. At both ends of the link, the 
transmit flags are set to 0 and the receive flags are set to 
-1. This receive flag state alerts the remote user to the 
fact that the link has been reinitialized and indicates that 
new receive bufiiers must be supplied as well as any 
interrupted transmissions restarted. When SETRX is 
called, the transmit flag at the opposite side of the link 
will be set to 2, initialization complete. Concurrently, the 
local receive flag will become 1, i.e., no message. A 
typical intialization sequence is shown in Fig. 3. The 
CALL TXSTAT checks for initialization complete. 
RXSTAT is called to detect and protect against a 
simultaneous initialization from the other end. 

IV. The Teletype Link 

Characteristics of the teletype circuit influenced the 
software design. Preliminary experiments indicated that a 
particular teletype circuit including the computer buffers, 
microwave, land lines, and associated equipment has a per 
character error rate of ICH. Two error checking schemes 
are used to ensure reliable communication. First, informa- 
tion to be transmitted is divided into four-bit bytes and 
encoded into the first foiu- bits of the teletype character. 
The fifth, bit of the teletype character is character parity. 
-Second, data are divided into blocks of 120 characters, and 
the characters within a block are bit-wise exclusive-ored 
to form a checksum parity character, or longitudinal 
'' parity check. If an error is detected, the block is 
retransmitted. A block size of 120 characters was selected 
as small enough to minimize the probability of retransmis- 
sion and yet large enough to prevent channel saturation by 
overhead bookkeeping characters. 

Two types of teletype characters are designated, data 
characters and utility characters. Constrained to be of 
even parity, data characters contain text and block 
overhead information. Odd parity characters are used as 
utility signals called “interhandlers” to control the 
interlocked transmission of blocks and synchronize the 
initialization of the system. The use of distinguishing 
parity permits the interleaving of interhandlers and text 
with no ill effect. 

The interhandlers consist of two odd parity characters 
transmitted consecutively twice. The five-bit characters 
are each denoted by two octal digits. Table 2 relates these 
numbers to teletyqje characters. For instance, the interhan- 
dler which notifies the remote teletype program to 
initialize is 37/01/37/01. Any combination of 37 and 01 
occurring within 4 characters will cause the interhandler 
to be recognized (e.g., the message 47/01/37/02 is valid). 


Since the interhandlers are selected to be of at least 
distance 2, this technique yields simple error correction. 
The interhandlers are summarized in Table 3. 

Interhandlers are given odd parity to permit them to be 
inserted into the text of a message. This prevents delay of 
up to one block time (13.7 seconds) in the transmission of 
an interhandler. If an odd parity character is received, its 
validity as an interhandler character is checked. If it is 
valid, it is accepted. If it is not valid and message receipt 
is in progress, an error flag is set to indicate that a 
character with improper parity was received. If a message 
is not in progress, the character is ignored. 

Whenever data or link parameters are transmitted, they 
are encoded into four-bit bytes. The data bits become the 
four most significant bits of a five-bit teletype character, 
with the remaining bit used to set even parity. 

Data are sent in blocks of 120 characters (480 bits). The 
format of a block is shown in Fig. 4. The start-of-block 
interhandler alerts the receiver to imminent data. 
Following immediately, a header information sub-block 
includes a message number, channel number, block 
number, and block character count. The message number 
for each channel is set to 0 during initialization and 
incremented by one for each message of a channel. It is 
transmitted as one character, hence it is modulo 16. The 
channel number identifies whether the block is assigned to 
channel 1 or 2. To order the message, a block number is 
initialized as the number of blocks in a message. It is 
allowed 2 characters (8 bits). As each block of a message is 
sent, the block number is decremented by one. The final 
block of a message, therefore, has a block number equal to 
one. Used in conjunction with the message number, the 
block number helps keep track of the sequence of blocks 
and messages as an error detection tool. The character 
count is 8 bits (two characters) containing the number of 
characters in the block. “Text” contains the data to be 
transmitted. End of text is marked by the interhandler 10/ 
25. It is followed by a checksum derived from a character- 
by-character bit-wise exclusive-or of the text and header 
information sub-block. Finally, the end of block is 
designated by 10/34/10/34. 

After a block is sent, the transmitter awaits a reply from 
the receiver before continuing. This reply or handshake is 
used to synchronize the link and control retransmission. As 
a block is received, an error flag is set (1) if any particular 
text or header sub-hlock character has odd parity, (2) if the 
count of received characters and the block character 
count disagree, (3) if the derived checksum parity and 
transmitted checksum parity disagree, (4) if the block 
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. number is out of sequence, (5) if the message number is 
changed prior to receipt of block number one, or (6) if the 
end of block is received prior to the end of text 
interhandler. If the error flag has not been set, the 
receiver replies to the transmitter with the block- 
received-OK interhandler; otherwise, block-not-okay is 
sent. At the transmitter, in the former case, the message is 
continued or if complete, the transmitter chaimel remains 
busy but inactive imtil the receiver signals that a new 
receive bufier-ts available. In the latter case, a parameter 
set by the user for that channel, designated RTRY, is 
decremented. If it remains >0, the block is retransmitted. 
IF RTRY is 0, the user is notified that the chaimel is in 
error and the channel is freed. 


V. The XDS 920 Switcher 

With the teletype link from the Venus Station XDS 930 
to the XDS 920 at JPL established, the next step is the 
path between the XDS 920 and the Xerox Sigma 5. These 
machines are, fortunately, already interconnected via the 
D/R high-speed link. Software to support the D/R unit is 
available and utilized. As this software will be described in 
a cursory fashion here, the reader is again referred to the 
article by Layland (Ref. 5) for more detail. 

A program to interface between the D/R software and 
the teletype software described earlier resides in the XDS 
920. This program, in combination with the Sigma 5 
software, discussed in the next section, appears transpar- 
ent to the user. Hence, the Sigma 5 user believes he is 
communicating directly with the teletype handlers. The 
main tasks, then, for the XDS 920 “switcher” program are 
to route and reformat message trafiic through the D/R 
link and to extend the teletype transmit and receive flap 
to the Sigma 5. 

Communication with the Sigma 5 is maintained via the 
aforementioned D/R link software. This software is 
capable of transferring blocks of 126 8-bit bytes at rates 
approaching 25 kilobytes per second. Seven logical 
bidirectional channels are provided. Five of these were 
used in this implementation. D/R channels 1 and 3 are 
assigned to teletype channels 1 and 2, respectively, of the 
teletype handler for trafiic from JPL to Coldstone. D/R 
channels 2 and 4 were assigned to teletype channels 1 and 
2, respectively, for Coldstone to JPL trafiic. D/R channel 
5 was used for “signals” to coordinate and control the XDS 
920/Sigma 5 link. 

As far as the user is concerned, the teletype handlers 
can transfer memory blocks of arbitrary size. On the other 
hand, the D/R unit software is limited to transferring 
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blocks of 126 8-bit bytes. Therefore, the switcher and 
Sigma 5 software must cooperate by disassembling 
messages into D/R unit blocks for transmission and by 
reassembling these blocks into core images upon recep- 
tion. To main tain compatibility with the teletype link, the 
Sigma 5 software must be informed of the receive 4-bit 
character coimt for the entire message. This is done by 
inserting the character count into the first two bytes, i.e., 
16 bits of each D/R block. Although this information is 
unnecessarily repeated, the time penalty is negligible. In 
the opposite direction, the switcher makes use of the 
SNDMSC subroutine to transmit trafiic to Coldstone. The 
subroutine must be supplied with channel priority, 
number of retrys and the four-bit character count. This 
information is supplied by the Sigma 5 user and sent to 
the switcher program as the first four bytes of each D/R 
block. 

It should be noted that the teletype channel priority 
scheme is not reflected in the D/R link. This is simply 
because of the high speed of that link. Whenever a full 
message is ready to be transmitted from the switcher to 
the Sigma 5, or vice versa, it is sent on a first-come first- 
serve basis without regard to priority. 

The signals used to coordinate and control the message 
trafiic over the D/R link are given in Table 4. These 
signals allow the Sigma 5 applications program to direct 
the teletype subroutines in the XDS 920 and interrogate 
their status much as if the teletype were directly under its 
control and not separated by two additional interface 
programs. The sequence of events in handling this link is 
identical to that discussed earlier. 

VI. The Sigma 5 End 

An overriding requirement in the communications link 
design was that both ends of the link appear identical to 
the user. The purpose of the Sigma 5 program is to 
establish an interface between the Fortran-callable 
subroutines and the D/R link software. This interface 
provides compatibility with the teletyp>e routines, splits 
arbitrary-length outgoing data messages into D/R blocks, 
and concatenates incoming D/R blocks into arbitrary- 
length messages. 

The environment of the Sigma 5 led to an immediate 
problem. The communications software is a subprogram 
of the user’s program. The D/R link’s software is a 
subroutine of the communications software. While the 
D/R software is interrupt-driven and therefore runs in real 
time, the communication subprograms run only when 
called by the user program. ’The communications software 
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must, however, chop messages into D/R blocks and “feed” 
them to the D/R subprograms. In the inverse, the 
communications software must concatenate D/R blocks 
into messages. To keep the link running, therefore, either 
repeated calls to some “keep running” subroutine or some 
“busy-wait” is required. The former alternative would 
violate the premise that both ends of the link be identical. 
The latter could cause possibly hazardous delays to the 
user program. 

To avoid this problem, the interface subroutines are 
constructed around a FIFO list, the task list. Whenever 
some subroutine has a pending operation, it places a 
particular number in the list and executes a return to the 
calling program. The internal subroutine STACKPUT(N) 
is used for ^is operation. To keep the system running, the 
subroutines TXSTAT and RXSTAT contain calls to 
STACKGET. STACKGET removes, in turn, the topmost 
item in the list and transfers control to the designated 
subroutine. When the subroutine execution is complete, 
control is returned to STAGKGET. This process is 
repeated until the task list is emptied. Gontrol then 
returns, as appropriate, to TXSTAT or RXSTAT. Since 
TXSTAT and RXSTAT are routinely called by a user 
program, and in fact are usually part of some testing loop. 


the use of the task list provides satisfactory operation 
which is transparent to the user. 

The communications link is normally operated from the 
timesharing terminal on the Sigma 5. Because of the 
design of the link, some time is spent in busy-waiting on 
the status of the transmit or receive flags. During these 
periods, it was deemed inappropriate to delay concurrent 
background or batch processing jobs. On the other hand, 
TXSTAT and RXSTAT must be repetitively called, as just 
described, to keep the system going. A subroutine, 
“ENDSLIGE,” was added to the link software to allow 
the user program to return the unused portion of a time- 
share time slice to the operating system. ENDSLIGE will 
return the slice only if the link is not actively receiving or 
transmitting, that is, only if repeated calls to TXSTAT or 
RXSTAT are not necessary. 

VII. Conclusion 

The communications system has now successfully 
operated for an elapsed period approaching 100 hours. It 
is therefore considered “operational” and will continue to 
be used in the station automation experiment as long as 
teletype is the primary digital data link with the Venus 
Deep Space Station. 
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Table 1. Summary of subroutines 


Table 2. Teletype characters 


SNDMSG(LOC.CHN,PRY,RTRY,CNT) 

Causes a data block to be transmitted. 

LOC — Location of the first word of data 

CHN — Channel 1 or 2 

PRY — Priority 1 to 15 

RTRY — Number of block retransmissions 

CNT — N umber of 4-bit bytes of data 

TXSTAT(CHN,F1.C) 

Returns transmit status of channel CHN as one of the 
following values for FLC; 

2 — Channel initialized 
1 — Message received, channel free 
0— Channel busy 

— 1 — Message not received (error), channel free 
RXSTAT(CHN,FLC) 

Returns receive status for channel CHN as one of the fol- 
lowing values for FLC; 

— 1 — No buffer 

0 — Message in 

1 — No message in 

SETRX(CHN,LOC,CNT) 

Assigns a receive buffer starting at LOC, CNT words in 
size, to channel CHN. 

TINIT 

Initializes the link. 

ENDSLICE Sigma 5 ONLY 

Releases tmused time slices to the operating system. 

MSGCNT(CHN,CNT) 

If the receive fiag for CHN is 0, CNT will be assigned a 
value equal to the number of 4-bit bytes just received. If 
the flag is not 0, then CNT will become —1. 


Letters 

character 

Octal 

code 

Figures 

character 

Octal 

code 

Letters 

37 

Figures 

33 

A 

03 

0 

26 

B 

31 

1 

'27 

C 

16 

2 

23 

D 

11 

3 

01 

E 

01 

4 

12 

F 

15 

5 

20 

G 

32 

6 

25 

H 

24 

7 

07 

I 

06 

8 

06 

] 

13 

9 

30 

K 

17 

- 

03 

L 

22 

? 

31 

M 

34 



N 

14 

$ 

11 

0 

30 

! 

15 

P 

26 

& 

32 

Q 

27 

# 

24 

R 

12 

* 

13 

S 

05 

( 

17 

T 

20 

) 

22 

U 

07 


34 

V 

36 


14 

w 

23 

; 

36 

X 

35 

/ 

35 

Y 

25 


21 

Z 

21 

Blank 

00 

Bell 

05 

Space- 

04 

Line feed” 

02 

Carriage return 

10 

“Figures character. 
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Table 3. Teletype link interhandlers 


Interhandler 

Teletype characters 

Order initialiae 

37/01/37/01 

Initialization complete 

37/31/37/31 

Block received OK 

37/26/37/26 

Block received in error 

37/23/37/23 

Start of block 

10/02/10/02 

End of text 

10/25 

End of block 

37/34/37/34 

Channel 1 buffer available 

37/15/37/15 

Channel 2 buffer available 

37/13/37/13 


Table 4. 0/R link signals 


Signal 

Number 

Order initialize 

0 

Initialization complete 

1 

Channel 1 free 

2 

Channel 2 free 

3 

Start message channel 1 

4 

Start message chaimel 2 

5 

End message chaimel 1 

6 

End message channel 2 

7 

Bad teletype data channel 1* 

8 

Bad teletype data channel 2“ 

9 

“920 to Sigma only. 
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JPL 


VENUS STATION 


OPBIATOR 



TELETYPE 

UNK 


Fig. 1. Communications system 


CALL TINIT 

3 CALLSETRXd, LOCI, CNT1) 

CALLSETRX (2, LOC2, CNT2) 

1 CALLTXSTAT (1, FLG) 

IF (FLG-2) I, 2, 1 

2 CALL RXSTAT (1, FLG) 

IF (FLG) 3, 4, 4 

4 CONTINUE 


Fig. 3. Initialization sequence 


TRANSMITTER 


1 CALLTXSTAT(1, IFLG) 

IF(IFLG)3, 1, 2 

2 CALL SNDMSG (KRAY (I), I, I, 2, 200) 


3 TYPE 4 

4 FORMAT (STRANSMBS ION ERRORS) 


RECEIVER 


1 CALL RXSTAT (1, IFLG) 

IF (IFLG) *, 2, 1 

2 CALL MSGCNT (1, CNT) 

CALL SETRX (1, KRAY (1), 2000) 


* THIS INITIALIZATION STATE IS DISCUSSED 
ELSEWHERE 


START OF BLOCK HEADER INFORMATION SUB-BLOCK 


HEADER 


02 

10 

3 



N 

N 

B 

3 

TEXT ' 

^ TEXT - MAXIMUM 120 CHARACTERS ' 

TRAIUR 

10 

25 1 

□ 


34 

10 

34 






r" ^ 


^ 








END OF TEXT END OF BLOCK 


M ^ MESSAGE NUMBER 
C ■= CHANNEL NUMBER 
N > 2-OIGIT BLOCK NUMBER 
B = 2-DIGIT CHARACTER COUNT 
P = LONGnUDlNAL PARITY 


Tig. 2. Sample message transfer programs 


Fig. 4. Block format 
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N76-18177 


Implementation of Wideband Digital 
Recording Equipment in the DSN 

K. R. Kimball 

DSN Data Systems Development Section 


The DSN is implementing wideband digital recording equipment for mission 
support in the 197&-1982 time era. This article describes the development status 
of this equipment and the factors which were of major importance in the design 
approach. Design of equipment is discussed to the functional block diagram 
level. 


I. Introduction 

Beginning in late 1977, wideband digital recording 
equipment will be implemented at the 64-meter Deep 
Space Stations (DSS). This equipment is necessary for 
support of flight project data acquisition and spacecraft 
navigation requirements. Specific^y, a 12-megabit/sec 
re^rding capability is required in support of the Pioneer 
Venus probe entry wind experiment in December 1978; 
the capability to record 24 Mbits/sec is required for the 
DSN to use very long baseline interferometry (VLBI) 
techniques for determination of Earth platform parame- 
ters (UT-1 and polar motion) to the accuracy necessary for 
spacecraft navigation in the 1978-1981 time era. 

Efforts toward implementing such a capability began in 
late 1974. Magnetic tape recording was adopted as the 
most acceptable means to accomplish digital recording at 
the required rates; multi-track recording on 2.5-cm (1- 
inch) wide magnetic tape, using an instrumentation-type 


tape transport, is the general approach used by the 
recording industry to achieve such data rates. 

II. System Design Factors and Influences 

A. Capability Avaiiable From Commercial Sources 

A key aspect of the implementation approach adopted 
by the DSN for wideband digital recording is to 
commercially procure from a single contractor that 
portion of the recording equipment which is to perform 
the primary functions of recording and reproducing user- 
supplied digital data at the rates and error tolerances 
required by the DSN. This approach places the responsi- 
bility for engineering those interfaces which are most 
critical to wideband digital recording in a single organiza- 
tion, and simplifies the specifications for the equipment 
being procured. 

One of the first activities in the design effort was to 
become familiar with the capability that the recording 
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industry had achieved in wideband digital recording. The 
following summarizes the key observations made during 
this efibrt: 

(1) Several industry sources claimed digital recording 
capability in excess of 50 Mbits/sec. 

(2) There were (and are) no standards within the 
industry for wideband digital recording (i.e., formats, 
densities, track widths, track spacing, encoding 
system, etc.) 

(3) Although advertised as such, no truly “off-the-shelf’ 
product line was offered by any industry source; this 
was in part due to the fact that users of such 
equipment had had widely varying requirements. 

(4) Most systems were still undergoing design changes 
and improvements. 

(5) There existed a large variance in the prices for 
equipment which ostensibly performed equivalent 
functions; in general the pricing of the equipment 
appeared high relative to the amount of electronics 
provided. 

(6) The specifications of system bit error rates were 
vague; specifically, techniques for distinguishing 
equipment errors from those errors resulting from 
tape imperfections were not well defined. 

Following the industry survey, a JPL specification for 
25-Mbit/sec recording equipment was generated. The 
technical requirements of this specification were written 
to be generally consistent with the level of capability 
observed during the industry survey; i.e., it was expected 
that several sources would be capable of supplying 
recording equipment in accordance with the JPL specifi- 
cation. 

B. Supporting Research and Technology Activities 

A second major influence on the implementation design 
approach was the related effort conducted by JPL’s 
Communications Systems Research Section. In July 1974, 
a purchase order was issued by this section for the 
procurement of an off-the-shelf 80-Mbps record/repro- 
duce system. The selected manufacturer encountered 
several difBculties in producing a system which met the 
required specifications. The specific technical outgrowths 
of this procurement which had an impact on the 
specification and selection of the DSN equipment were as 
follows: 

(1) A conservative design approach would preclude the 
use of tape track widths less than 25 mils for digital 
recording. 


(2) The procedures and techniques for testing bit error 
rates must be carefully specified to ensure that a 
truly parallel bit error rate is measured, and to 
eliminate the effects of tape flaws. 

(3) The minimization of electrical and mechanical 
adjustments is vital if the system is to be operatio- 
nally acceptable. 

C. Related Efforts Within NASA 

Similar requirements for wideband digital recording 
exist within NASA for several planned and proposed 
applications. The rate requirements for these applications 
fall in the range of 15 to 120 Mbits/sec. Cognizance of 
these related efforts is under NASA’s Goddard Space 
Flight Center (GSFC). 

In. March 1975, a JPL-GSFC wideband recorder 
technology working sub-group was established (as part of a 
larger VLBI working group); the primary objectives of this 
sub-group were to: 

(1) Elxchange information as to development status and 
design approach. 

(2) Assess the degree to which format and interface 
standards could be adopted to permit inter-network 
experiments. 

(3) Assess the degree of operational risk associated with 
the relatively new level of technology present in 
wideband digital recording systems. 

The first meeting of the working group was very 
productive in terms of information exchange, and there 
were several areas of agreement as to desirable interface 
and format characteristics. Although JPL had by this time 
selected a system contractor and the design of the 
equipment was well underway, several aspects of the 
detailed design were influenced by the results of this 
meeting. It is believed that the DSN equipment as 
currently designed could be readily compatible with any 
potential “NASA-standard” system. 

D. Operational and Reliability Considerations 

The most common reservations generally expressed 
about recording equipment center around operability, 
maintainability, and reliability. In an attempt to be 
responsive to these concerns, the following design goals 
were established and have been of major importance in 
the system design: 

(1) An automated pre-pass performance verification of 
the recording equipment should be provided. 


JPL DEEP SPACE NETWORK PROGRESS REPORT 42-31 


97 



(2) Some degree of real-time read-after-write data 
monitoring capability should be provided. 

(3) Tape-changing should be the only real-time operator 
function; the functions of starting, stopping, and 
rewinding tape should be automated. 

(4) A hard-copy log of all recording activity should be 
provided. 

(5) The capability should be provided to receive control 
information from a remote operator, and to report 
system status to station monitoring equipment. 

(6) The recorder transport should accommodate 41-cm 
(16-inch) reels for extended recording time. 

(7) DSN Hi-Rel hardware should be used wherever 
appropriate. 

(8) The recording system MTBF should be >500 hr. 

(9) Modular design should be provided to facilitate 
troubleshooting and repair. 

(10) Mechanical and electrical adjustments should be 
minimized. 

E. Cost Considerations 

The effort to minimize costs while providing capability 
consistent with requirements has had an impact upon the 
design in two significant respects. The first was the 
procmement of the commercial recording equipment in 
two distinct configurations. One configuration is intended 
for use at DSSs, where data recording is the primary 
function, and the other is destined for use at a central 
playback facility. The “record” configuration has full 25- 
Mbit/sec recording capability but limited reproduce 
capability; the “reproduce” configuration has complete 
playback capability but only that record capability 
required for system performance checks. An estimated 20- 
percent reduction in price of the commercial recording 
equipment was achieved by procuring units without the 
full record/reproduce function. 

Equipment costs were also minimized by implementing 
only 20 channels of electronics in recording equipment 
which is designed to accommodate 28 chaimels; the digital 
recording rates required for the DSN application can be 
comfortably accommodated on 20 recorded tracks. 
Although all systems are equipped with 28-track read/ 
write heads and 28 channels of head drivers and reproduce 
amplifiers, those channel electronics which exhibit greater 
complexity and cost (record and reproduce amplifiers, 
encoders, and decoders) have been limited to the 20- 
channel configuration. 
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III. Detailed Design Description 

A. Commercial Wideband Digital Recording Equipment 

As outlined above, two configurations of recording 
equipment are being procured from a systems contractor. 
A functional block Vagram of the record configuration is 
shown in Fig. 1. The input to the recording equipment 
consists of 20 synchronous channels of digital data and a 
common data clock. The nominal data rates per channel 
for the DSN configuration are shoWn in Table 1. 

Each input channel is recorded on a separate track after 
being reformatted by the encode/sync-insertion logic. The 
reformatting process consists of merging additional data 
bits (on a non-replacement basis) with the input data for 
each chaimel. For every 200 user input bits received on a 
given channel, an additional 43 bits are generated auid 
merged to create a 243-bit data frame. The frame content 
is shown in Fig. 2. As a result of the reformatting process, 
the actual bit rate of information recorded on each track is 
21.5 percent higher than the original input bit rate. The 
primary functions of the reformatting are (1) to provide a 
means for removing skew between tracks during playback, 
(2) to limit the low-frequency content of the recorded 
serial bit stream, and (3) to provide a means for error 
detection during reproduction. 

A limited read-after-write capability also exists in the 
record configuration. Two reproduce-decoder chaimels 
are provided; one of these chaimels is hardwired to a 
specific track, while the other can be switched via a 
multiplexer to any of the 20 recorded tracks. This 
multiplexing capability allows real-time monitoring of 
data quality on a track-selective basis. The functions of the 
decoder channels are (1) to eliminate data skew between 
recorded cbaimeb, (2) to remove the previously merged 
“overhead” bits and reproduce the data as it originally 
appeared at the user input, (3) to produce an error output 
signal which will indicate the bit error rate of the 
respective channel, and (4) to provide a signal which 
indicates the synchronization status of each decoder 
channel. 

A block diagram of the reproduce configuration, 
intended for use at a central playback facility, is shown in 
Fig. 3. The individual elements of this configuration are 
equivalent in function to those of the record configuration. 
A full complement of reproduce amplifiers and decoders is 
provided for playback of 20 data tracks. The decoders and 
deskew electronics produce 20 synchronous channels of 
digital data. The content and phase relationship of these 
data is identical to that supplied by the user during the 
record process; the fact that overhead bits were added and 
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subsequently deleted during the record-reproduce process 
is transparent to the user. 

An externally supplied data clock controls the data 
output frequency within a ±2-percent range of the 
nominal playback rates. Data can be reproduced at the 
originally recorded rate or at any of the nominal rates 
below the original rate. 

B. Equipment Configuration at DSN Sites 

A functional block diagram of the planned equipment 
configuration at DSN sites is shown in Fig. 4. The basic 
functional elements of the system are: 

(1) Analog-to-digital (A-D) converters and data 
formatter 

(2) Wideband recording equipment (2 units) 

(3) Recording monitor unit 

(4) Recording controller and log printer 

The primary functions and characteristics of these 
elements are outlined below. 

The A-D converters and data formatter will perform 
the following functions: 

(1) Digitize up to 6 analog inputs; single-bit quantiza- 
tion is planned for very long baseline interferometry 
(VLBI) applications, while 3-bit quantization of a 
single input will be provided for the Pioneer Venus 
wind experiment. 

(2) Derive A-D sampling clocks from the DSN hydro- 
gen maser frequency reference; sampling rates up to 
25 MHz can be provided with less than 3 nanosec- 
onds sampling jitter. 

(3) Provide a utility input port which will allow any 
external digital data stream to be recorded. 

(4) Distribute the digitized data onto 18 parallel data 
lines for input to the recording equipment. 

(5) Generate two redundant channeb of digital house- 
keeping data for input to the recording equipment. 
Housekeeping tracks will include data time tag, 
system configuration data, and auxiliary data re- 
ceived from the recording controller. 

(6) Provide a data generator which will serve as a data 
source for stand-alone performance evaluation of the 
equipment. 

The functional characteristics of the “commercial 
recording equipment” shown in Fig. 4 were discussed in 
detail in Section IIIA. This equipment (transport and 


electronics) will be controlled and monitored via a remote 
control interface connection to the recording controller. 

The function at the recording monitor unit is to 
accumulate, display, and transmit to the. controller data 
which reflects the status and condition of the recording 
process. Specific functional tasks performed are: 

(1) Selection of the multiplexed channel monitored 
during the read-after-write process. 

(2) Accumulation and display of statistics on the 
selected channel, and input of these statistics to the 
controller for verification of prof>er values. 

(3) Monitoring of the housekeeping tracks and extrac- 
tion of the time tag information for display and 
verification of proper content by the controller. 

(4) Reflection of the status and activity of the tape 
transports. 

It is planned for the recording controller, in conjunction 
with the log printer, to perform the following functions: 

(1) Stop, start, and sequence transports and rewind 
tapes when appropriate. 

(2) Monitor real-time data quality. 

(3) Receive auxiliary data (e.g., weather, radiometer, 
antenna position, etc.) through a DSN standard 
interface port and reformat these data for inclusion 
on the recorded housekeeping data tracks. 

(4) Detect system alarm conditions and transmit alarm 
status via the standard interface port. 

(5) Provide a hard-copy log of recording activity and 
performance characteristics. 

(6) Provide an automated pre-p>ass performance verifi- 
cation test. 


IV. Installation Schedule 

The current implementation schedule calls for installa- 
tion and integration testing to be completed as indicated 
below: 

DSS 14: 1 January 1978 
DSS 43: 1 February 1978 
Santiago (STDN station): 1 January 1978 
Guam (STDN station): 1 February 1978 
DSS 63: 1 March 1979 
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The equipment installed at Santiago and Guam is being 
loaned to the Spaceflight Tracking and Data Network 
(STDN) for support of the Pioneer Venus entry probe 
wind eqieriment and will be returned to the DSN in 
January 1979. 


V. Potential for Expanded Recording 
Capability 

An attempt has been made throu^KUit the design 
proc^ to provide a recording system which is general- 
purpose in nature, and to allow for future expansion or 
modification of recording capability. Some of the potential 
requirements for which such equipment could be used 
are: 

(1) Recording of digitized pre-detection telemetry data. 

(2) Recording of radio science data. 

(3) Post-detection recording of multi-megabit telemetry 
data streams. 

Increases in recorded data rate capacity could be 
accommodated by any of the following means: 

(1) Recording data at a tape speed of 305 cm/sec (120 
inches/sec), yielding an aggregate 55.6-Mbit/sec 


capacity; recording time per reel, however, is 
recced to 20 minutes at this tape speed. 

(2) Recording data at a tape speed of 610 cm/sec (240 
inches/sec), yielding an 111-Mbit/sec capability. 
The operational feasibility of this tape speed is 
questionable, as reel changes are required at 10- 
minute intervals (assuming continuous recording on 
41-cm (16-in.) reels). A second drawback is that the 
recording systems being supplied to the DSN will 
have no reproduce capability at 610 cm/sec. Data 
recorded at 610 cm/sec must be played back at 305 
cm/sec. 

(3) Increasing the number of recorded tracks from 20 to 
28. This would result in a 40% increase in capability 
at any tape speed. 

(4) Increasing on-tape bit density from the current 11 
kbits/cm (28 kbits/inch) to 13 kbits/cm (33 kbits/ 
inch). 

All of the above are technically feasible, and there exists 
a potential system rate capability in the 150-Mbit/sec 
range. However, any future consideration of increasing the 
data rate capacity of the system must be heavily 
dependent on (1) the experience gained with the 
equipment in its present configuration, (2) operational 
considerations, and (3) equipment costs. 
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Table 1. Nominal DSN recording rates 


Tape speed, 
on/sec (in./s) 

Input bit rate 
(per trade). 
Mbits/sec 

Aggregate bit 
rate (20 tracks). 
Mbits/sec 

305 (120) 

2.78 

55.6 

152 (60) 

1.39 

27.8 

76 (30) 

0.694 

13.9 

38 (.15) 

0.347 

6.94 
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S-Band-X-Band Directional Coupler 

R. Lay 

R. F. Systems Development Section 


In support of an effort to reduce microwave noise bursts at the DSN tracking 
stations, an S-band-X-band directional coupler has been fabricated. The goal 
was to develop a directional coupler that has a coupling of about 30 dB using 
two different rectangular waveguide sizes. 


I. introduction 

The DSN tracking stations have experienced occasional 
problems with excessive microwave noise bursts appear- 
ing in the low-noise receivers during high-power diplexing 
at S-band. Locating the source of such noise bursts has 
often proved to be quite difiBcult. The S-band-X-band 
directional coupler, shown in Figure 1, has been 
investigated as a possible aid in locating the source of the 
noise bursts. 

It is assumed that the electrical breakdown that 
generates the noise burst has a broad band spectrum. The 
directional coupler couples the X-baiid component of the 
noise burst through WR 112 waveguide to a test receiver 
while, in efiect, completely isolating the test receiver from 
the high-level S-band being transmitted in WR 430 
waveguide. By inserting the directional coupler at various 
locations in the waveguide system, it should be possible to 
isolate the source of noise bursts to a particular 
component. 

II. Design Considerations 

With reference to Figure 2, the directional coupler 
consists of two unequal size rectangular waveguides 


coupled together by means of a series of small coupling 
aperttues on the broad walls between the two guides. The 
coupling theory is based on the diffraction of electromag- 
netic radiation by an aperture, where its size is small 
compared with the wave length. An incident wave in port 
1 couples power into ports 3 and 4. 

By proper adjustment of the aperture position d, the 
radiation in the direction of port 4 can be minimized and 
that in the direction of port 3 enhanced. The amplitude 
equations of the excited fields in the secondary waveguide 
are given in both References 1 and 2 in terms of its 
exciting and excited fields. In the case of a circular 
aperture, they are proportional to the cube of the radius. 
With the aid of a computer program the four amplitude 
equations have been solved for the proper aperture 
position d in order to obtain the maximum radiation into 
port 3 using a single-hole coupler model. The apertme 
position d •= 3.1750 cm (1.25 inches) was chosen, and a 
circular aperture of radius Tq = 0.436563 cm (0.171873 
inches). Since the single-hole coupler model has a narrow 
bandwidth, many apertures were used in order to obtain 
better directivity over a larger bandwidth. There is a 
drawback in couplers using different sizes of coupling 
apertures because the ratio of the coupling coefficients is 
not constant over the frequency band. Thus, equal size 
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apertures were used, and equal spacing between the 
apertures was used for simplicity. Fifteen circular 
apertures were used, and the spacing between them, D <= 
1.0320 cm (0.80 inches) was chosen. 

III. Tests 

A waveguide coupler test sample has been fabricated. 
Low-power tests indicate that both good coupling and 
directivity were found over several bands of frequencies. 
Using a waveguide taper from X-band to S-ban(i it was 
possible to propagate an X-band signal inside the WR 430 


waveguide of the coupler. With reference to Figure 2, the 
output of the excited waveguide (WR 112) was plotted on 
the X-Y plotter over, a bandwidth as sho'> n in Figure 3. 
Port 4 monitors the forward direction i.i the radiated 
power from port 1 and port 3 monitors the reverse 
direction of the radiated power from port 1. The 
directivity is obtained by subtracting one from the other. 

With the test equipment available, a desirable fre- 
quency range would be the 8407-8417 MHz, which would 
afford about 30-dB coupling and 22-dB directivity over the 
10-MHz bandwidth (Figure 3). 
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Heat Transfer Criteria of a Tubuiar Solar Collector- 
The Effect of Reversing the Flow Pattern 
on Collector Performance 

F. L. Lansing 

DSN Engineering Section 


A new evacuated tubular solar collector has been selected for further investiga- 
tion for the application of Goldstone energy conservation projects. This article 
presents briefly the exact heat transfer analysis of this two-pass-flow collector in an 
effort to determine the difference in performance characteristics with two different 
flow patterns. The results from this analysis are not only the determination of the 
collected heat rate and the temperature profiles at each cross section but also the 
prediction of the maximum attainable fluid temperature at zero flow rate in both 
cases. 


I. introduction 

Part of the updated technology in energy conservation 
measures, using solar-assisted equipment, is the develop- 
ment of a new evacuated tubular two-pass-flow solar 
collector. This is a prospective candidate collector with a 
relatively higher, thermal performance than the flat-plate 
type. The collector has been selected by the DSN Engi- 
neering Section for further theoretical and practical inves- 
tigation for application in Goldstone energy conservation 
projects. 

Limited laboratory measurements were the only data 
available for the performance of this collector. The latest 


attempt (Ref. 1) at NASA Lewis Research Center using a 
solar simulator and an analogous one-pass-flow tubular 
collector was found insufficient to provide the relevant 
performance characteristics of the actual two-pass-flow 
tubular coUector. The conunon analysis procedure (Ref. 
2) of assuming an overall, one-directional heat loss to the 
surroundings in a flat-plate collector expression is no 
longer adequate for exact performance analysis with two 
different kinds of flow patterns. This suggested the need 
for exact heat transfer and flow analyses as an essential 
step in the determination of performance characteristics 
over wide ranges of irradiancies, fluid flow rates, inlet 
fluid temperatures and ambient conditions. This article 
presents the steps and results of these analyses. 
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II. Collector Description and Flow Patterns 

Tlie tubular collector, as shown schematically in Figs, 
la and 2a, is composed of three all-glass concentric 
tubes, namely, a feeder tube, an absorber tube and a 
cover tube. The annulus space between the absorber and 
cover tubes -is evacuated to such a low pressure 
torr) that convection and conduction losses are negligible. 
The absorber tube surface is coated with a selective ma- 
terial which minimizes outward long-wave radiation 
losses. In flow pattern (1), as shown in Fig. 1, the fluid, 
normally water, starts from the inlet section of the feeder 
tube and absorbs some of the useful solar energy along its 
path, thus raising its temperature. At the closed end of the 
collector, the fluid reverses its direction and passes in the 
annulus spacing between the feeder and absorber tubes. 
The temperature distribution along one streamline is as 
shown in Fig. lb. 

In flow pattern (2), as shown in Fig. 2, the fluid path is 
reversed from the above. It starts from one end of the 
collector through the annulus spacing, turns around at the 
closed end, and then leaves from the other end from the 
feeder tube. Both flow patterns (I) and (2) alternate, in 
each collector module, as shown in Fig. 3. 

For more irradiancy augmentation, the set of collectors 
is mounted, with lateral spaces separating them from each 
other with a highly reflective back reflector, as shown in 
Fig. 4. This increases the overall irradiancy (Refs. 1 and 
2) by as much as 64 percent at normal solar incidence. 


dQr, = the long-wave radiation exchange between the 
cover and absorber tubes. (This is absorbed by 
the cover, thus raising its temperature.) 

dQr, = the heat transfer to fluid (2) in the annulus pass 
from the absorber tube surface. 

dQ- = the augmented radiation and convection loss 
between the cover tube and ambient air and sky. 

dQ„ = the heat transfer to fluid (I) in the feeder tube. 

dQa = the sensible heat gain to fluid (2) in the aimulus 
pass. 

For both flow patterns, at steady state, the heat balance 
of the absorber tube alone can be expressed as 

dQ, - dQ, - dg. = 0 (1) 

and the heat balance of the cover tube alone can be 
written as 

dQ, + dQ, -dQ, = 0 (2) 

By eliminating the cover tube temperature, T,(x), and 
the absorber tube temperature, T,(x), from the heat bal- 
ance equations for fluids (I) and (2), the problem is re- 
duced to the solution of the two simultaneous differential 
equations, as follows: 


III. Analysis 

For a segment of the collector of length dx at steady- 
state conditions, the heat flux is divided as shown in 
Fig. 5, where 

dQi = the total irradiation on the cover tube from all 
sides (including irradiancies from the back re- 
flector and reflections from adjacent collectors). 

dQ, = the energy absorbed by the absorber tube (in- 
cluding absorption from multiple reflections 
within the evacuated space). 

dQ, = the total outward reflection loss from all sides of 
the cover tube. 

dQt = the energy absorbed by the cover tube (includ- 
ing absorption from multiple reflections within 
the evacuated space). 


dUx) 

dx 

dT,{x) 

dx 


+ K,T,{x) = K,T,{x) 

+ K,T,{x) = K,T,{x) + K, 


and 


dT,{x) 

dx 

dT,{x) 

dx 


-I- K,T,{x) = K,T,(x) 

+ K,7,[x) = K,T,{x) -f K, 


I 


for flow 
pattern (I) 


for flow 
pattern (2) 


( 3 ) 


( 4 ) 


where K„ K,, K,, and K, are constants detennined from 
the steady-state thermal analysis, and T,{x), T,(x) are the 
fluid temperatures in passes (1) and (2), respectively, at a 
distance x from the inlet section. 
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IV. Boundary Conditions 

(1) r,(x) = Ts(x)atx = L for both flow patterns 
(l)and(2) 

(2) r,(x) = r,(0)atx = 0 for flow pattern (1) 

(3) Tx(x) = TxiQ) at X = 0 for flow pattern (2) 

The temperature distribution T,(x), Tx(x) of the fluid 
for pattern (1) or T,(x), T 2 (x) for pattern (2) can now be 
determined. 

For flow pattern (1), and at the inlet section (x = 0), 
the temperatiure of the fluid leaving the collector 72(0) is 
given by 


From Eqs. (6) and (8), it can be proved fliat the roots s, 
and si are related to the roots Si,S 2 by 

Si = Si 

Tx = - Si > (9) 

K,(s,S 2 ) + (s, + Si) = 0 ) 

Substituting Eq. (9) in Eq. (7) will show that the expres- 
sions, Eq. (5) and Eq. (7), for the leaving fluid tempera- 
ture in both fluid patterns are identical. 

The results from this analysis are the determination of 
the collected heat rate and also the temperature distribu- 
tion of the fluid as it flows in each pass in addition to the 
absorber and cover temperatiu-es. 


T.(0) = r.(0) + |[ri(0)--^^]|- 

[ 1 - exp [(s, - S 2 )L] 

Si 

1 exp [(s, - S 2 )L] 

where Si and Si are the roots of the quadratic auxiliary 
equation: 



six - (Ka - Ki)s,.2 - K.(K3 - K.) = 0 (6) 


For flow pattern (2) and at the inlet section (x = 0), 
the temperature of the fluid leaving the collector 7,(0) is 
given by 


V. Maximum Fluid Temperature for Zero Flow 
Rate 

The constants K„ Kj. Ks, K„ s,, Sa, sj, and si are found to 
approach an infinity value as the mass flow rate ap- 
proaches zero. However, their quotients have finite values 
independent of the fluid rate. Accordingly, it is always 
possible to find finite values for the ratios Si/st, T,/sl, 
K,/(Ki — Ki), J,/Ki, and s,/K,. Using L’Hopital’s rule for 
indefinite quantities, the maximum fluid temperature, 
72(0), for both flow patterns will be the limiting value, 
from Eq. (5) or Eq. (7), as the mass flow rate approaches 
zero. This is given by 

[72(0)] „„ = 7,(0) + [t,(0) - 

(10) 


T,(0) = T,(0)-|[r,(0)-.;^;^]i 


1 - exp [(s, - S 2 )L] 


(1 — -^exp [(s, - S 2 )U] + -^(exp [(s, - s,)L] - 1) 
. ^2 


(7) 


VI. Numerical Example 

The following numerical example will show the dif- 
ference in performance between the two flow patterns of 
the tubular solar collector. The data were arbitrarily 
abstracted to be as close as possible to actual running 
conditions. However, the conclusions may be generalized 
at any other relevant conditions. 

Irradiation with normal incidence to back reflector 


where s, and Sx are the roots of the quadratic auxiliary 
pnnatiODl 

six + (Ka - K ,)?,.2 - K,(Ka - K.) = 0 (8) 


Inlet water temperature 

7,(0) in pattern (1) '1 70 

or 7,(0) in pattern (2) J 

Ambient temperature 7o 30 °C 


110 


JPL DEEP SPACE NETWORK PROGRESS REPORT 42-31 



Sky temperature 

4'>C 

Irradiation intensity 

0.75 kW/m’ 

Tube length 

1.067 m 

Feeder tube diameter 

0.029 m 

Absorber tube diameter 

0.041m 

Cover tube diameter 

0.51m 

Glass refraction index 

1.526 

Cover tube transmissivity 

0.91 

Cover tube emissivity 

0.90 

Absorber tube absorptivity 

0.85 

Absorber tube emissivity 

0.10 

Mass flow rate 

5 kg/hr 

Area of back reflector/ 
single collector 

0.10645 m’* 

Distance between two 
consecutive cover tubes 

0.05 m 

Distance between cover 
tubes and back reflector 

0.05 m 

Wind speed 

11 km/hr 

Water specific heat 

11.634 X 10-‘kWhAgK 

A. Collector Efficiency 


The relevant constants K„ 
culated and given as 

Ki, K}, and K, were cal- 

K. = 0.88532 m-‘ 

K 3 = 0.90242 m-’ 

K, = 2.50333 m-’ 

K, = 13.44064 K/m 

The roots of Eq. ( 6 ) were given by 

s, = -1- 0.13189 m-’ 

= — 0.11479 m"’ 

and those of Eq. ( 8 ) were given by 

T^= - 0.13189 m-’ 

57= -f- 0.11479 m-’ 


Substituting in Eqs. (5) and (7), the leaving temperature 
T 2 ( 0 ) in flow pattern (1) or T,(0) in flow pattern (2) was 
foimd to be the same in both cases and equal to 350.96 K 
{77.96‘C). 
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The instantaneous collector efficiency, based on the 
irradiation on the back reflector, was then calculated as 

(77.96 - 70) X 5 X 11.634 X 10-* 
collector efficiency 0.75 X 0.10645 

= 58% 

for both flow patterns. 

B. Temperature Distribution 

For this example, the absolute temperature of the fluid 
in the feeder tube T,(x), the fluid in the aimulus spacing 
Ti(x), the absorber surface temperature T,{x) and the 
cover surface temperature Ts(x) were given by 

flow pattern ( 1 ) 

T,(x) = 786.026 - 177.572 exp (0.13189x) 

- 265.454 exp ( -0.11479x) K 

T 2 (x) = 786.026 - 204.025 exp (0.13189x) 

- 231.036 exp ( -0.11479x) K 

T,(x) = 5.369 + 0.993 T,(x) K 
T,(x) = 290.47 + 0.026 T,(x)K 
flow pattern ( 2 ) 

,f,(x) = 786.026 - 231.036 exp (-0.13189X) 

- 204.025 exp (0.11479x) K 

%{x) = 786.026 - 285.454 exp(-0.13189x) 

- 177.572 exp (0.11479x) K 

f.(x) = 5.369 + 0.993 T^ix) K 

r,(x) = 290.47 + 0.026 T,{x) K 

The temperature profile was plotted in Fig. 6 for each 
flow pattern. It is obvious that in the case of flow pattern 
( 2 ), the fluid in the annulus spacing at the closed end 
(x = L) will reach a higher temperature than that in flow 
pattern (1-). This is due to the augmenting heat transfer 
from solar irradiancy and from the relatively hotter fluid 
in the feeder tube. 

C. Maximum Temperature for Zero Flow Rate 

Since the running conditions in this example were ini- 
tially characterized by a slow flow rate (5 kg/hr), it was 
reasonable to assume, without great loss of accuracy, that 
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the convective and radiative heat transfer coefficients do 
not change appreciably by further decrease of the flow 
rate. Accordingly, the ratio s./Ki approaches —0.12965 
and Kt/(K., — K,) approaches 786.026 K as the mass flow 
rate approaches zero. The limit of the temperature dif- 
ference [T2(0) — T,(0)] fromEq. (10) approaches 57.44 *C 
as the mass flow rate approaches zero. This means that it 
is possible for the previous operating conditions and both 
flow patterns to reach a leaving fluid temperature of 
127.44°C by one collector with an extremely slow flow 
rate. 


VII. Performance Comparison 

The results of the above steady-state analysis can be 
summarized as follows: 

(1) The working fluid with different flow patterns (1) 
and (2) in the tubular solar collector yields the same 


value of temperature difference only across the 
open end of the collector. This result is a useful 
tool in simplifying the collector module perfor- 
mance by adopting one analytical expression for both 
flow patterns between inlet and leaving fluid tem- 
peratures. 

(2) The design of selective coating and its durability 
should take into consideration the larger thermal 
gradient and the higher absorber tube temperahure 
in flow pattern (2) compared with that in flow pat- 
tern (1). The relatively faster degradation of the 
coating of flow pattern (2) compared with that of 
flow pattern (1) will affect the maintenance costs. 

Having established an analytic expression for collector 
performance enables the study of future topics. Examples 
of these topics are the effects of introducing different 
types of working fluids or setting different collector di- 
mensions. 
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AMBIENT TEMPERATUIE Tq 
SKY TEMBEBATUBE T. 


T,(x) FEEDEB TUBE FLUID 
T2(x) fluid in pass (2) 
TM ABSOBBEB TUBE TEMPEBATUBE 
' Tj(x) COVEB TUBE TEMPEBATUBE 


OUTLET 


INLE' 



Fig. la. Collector configuration (flow pattern 1) 


AMBIENT TEMPERATURE: Tq 
SKY TCMPERATURE: T 


T, (x) FEEDER TUBE FLUID 
Tj (x) FLUID IN PASS (2) 
T^ (x) ABSORBER TUBE HMPERATURE 
Tj (x) COVER TUBE HMPERATURE 


INLET 


OUTLET 



Fig. 2a. Collector configuration (flow pattern 2) 



Fig. lb. Temperature distribution along one 
collector (flow pattern 1) 



O DISTANCE ALONG TUBE x L 


Fig. 2b. Temperature distribution along one 
collector (flow pattern 2) 
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SOLAR RADIATION 



SOLAR IRRAOIANCY 



Fig. 4. Arrangement of two collectors with a flat back 
reflector 




Fig. 6. Temperature distribution along a streamline with a 
flow pattern (1) and (2) 
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Tracking Operations During the Viking 2 
Launch Phase 

J. A. Wackley 

Network Operations Office 


Vising 2 launch phase operational procedures, and in particular the DSN 
initial acquisition procedures, were very intensively considered and conserva- 
tively designed to accommodate even the most unfavorable of launch 
possibilities. These procedures were successfully implemented and strongly 
contributed to the highly successful launch of Viking 2. 


I. Introduction 

The Viking 2 spacecraft was launched from the Air 
Force Eastern Test Range (AFETR) at 18:38:59.96 
Greenwich Mean Time (GMT) on 9 September, 1975, at a 
launch azimuth of 96.507 degrees. Like its sister 
spacecraft. Viking 1, the mission of Viking 2 is to study 
the planet Mars via direct measurements in its atmo- 
sphere, on its surface, and in orbit aroimd it. After being 
placed in a parking orbit by the combination Titan III- 
E/Centaur D-IT launch vehicle. Viking 2 was injected 
into a trans-Mars, heliocentric transfer orbit over southern 
Africa. The resulting near-Earth portion of the orbit was 
such that, within the Deep Space Network (DSN), the 
Australian complex of tracking stations was first to view 
the spacepraft. The Deep Space Station (DSS) selected to 
perform the initial acquisition was DSS 42, with DSS 44 
serving as a backup. 

In the following sections, the prelaunch tracking 
operations plaiming wiU be reviewed and an analysis of 
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the subsequent laimch phase tracking operations at DSS 
42 and DSS 44 will be presented. 


II. Review of the Viking Acquisition Strategy 

Since Viking 2 was similar to Viking 1 in both the 
trajectory it was to follow and in the characteristics of its 
radio frequency subsystem, it was not necessary to make 
any major revisions to the acquisition plan used during the 
Viking 1 launch (described in detail in Reference 1). The 
main points of this plan were: 

(1) Antenna to be driven in computer mode imtil after 
confirmation of two-way. In order to accomplish 
this, tracking predictions based within two seconds 
of the actual liftoff time were needed by the 
stations. Thus, with this necessity in mind, the new 
Polynomial CoefiBcient Tape-Network Operations 
Control Center Prediction software (PCT-PREIDK) 
prediction generation system was carefully exercised 
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in prelaunch testing and an elaborate PCX genera- 
tion plan was developed (see Reference 1). By the 
day of launch, the throughput time of the system 
had been reduced from approximately 41 minutes to 
less than 30 minutes. 

(2) Uplink acquisition to commence at approximately 
rise plus three minutes in order to collect as much 
two-way doppler data as possible before the 
interruption in data to be caused by the uplink 
transfer to DSS 44. 

(3) Uplink acquisition sweep to encompass XA ± 110 
Hz by use of the voltage-controlled oscillator (VCO). 
The extent of the sweep was chosen based upon the 
combination of the trajectory uncertainties and the 
spacecraft receiver imcertainties provided by the 
Orbiter Performance Analysis Group (OPAG) these 
being (at VCO level): 

30- XA TRAJECTORY ~ .05 Hz 
3<r XA S/C RECEIVER - 15 Hz 
3o- XA S/C “RANDOM WALK” - 7.5 Hz 
3<r XA S/C VCO TEMPERATURE - 5.6”C 
and 

axA/ar - 9 Hzz-c 

which when combined give a total 3<r uncertainty of 
approximately 53 Hz. Thus the sweep was designed 
to encompass approximately twice the 3or values 
about the XA. Figure 1 details the XA frequency 
versus time. 

(4) Uplink tuning rate to be 3 Hz/sec (VCO). This rate 
was considered to be close to the limit of the 
accurate manual tuning capability and when re- 
ceived at the spacecraft would result in an efifective 
tuning rate of about 280 Hz/sec, approximately the 
geometric mean of the allowable tuning rates. 

Additionally, the open window launch trajectory for 9 
September resulted in maximum angular and frequency 
rates of: 

d/dt I HA) 2E .06 degree/second 
d/dt ID2) 2 119 Hz/second (S-band) 
d/dt IXA) 2: .6 Hz/second (VCO) 
where 

HA = Local (station) hour angle 

D2 = Two-way doppler frequency 

XA = Spacecraft receiver best lock with doppler 
accounted for 


While these rates are higher than those encountered 
during the Viking 1 launch, they are much lower than 
those encountered in previous parking orbit ascent 
trajectories. Figure 2 is a stereographic illustration of the 
launch pass over DSS 42. The EUevation An^e versus time 
is shown in Figure 3, which also serves as a time line for 
the important tracking events at DSS 42. 

III. Postflight Analysis of the Viking 2 
Launch Phase 
A. Tracking Predictions 

1. The PCT-PREDIK Prediction Generation System. 
The PCT-PREDIK System functioned smoothly and 
efficiently during the Viking 2 launch phase. Because the 
Viking Flight Path Analysis Group (FPAG) was able to 
deliver the first PCT at launch minus two hours 30 
minutes. Stations 42 and 44 had more than enough time to 
generate a drive tape well before launch using predicts set 
B09D. During the launch countdown, changes made to the 
predicted frequency were small enough (XMTREF was 
changed by 3 Hz) that it was decided not to generate the 
planned frequency update predicts (predicts set B09F). 

As was the case in the Viking 1 launch, the remaining 
throughput time problems and concerns were allayed 
when launch occurred within a fraction of a second of the 
expected time. This left only one predicts set to be 
generated between launch and spacecraft rise; this set 
updated the frequencies and gave the stations text predicts 
with a GMT time field (all previous text predicts had been 
generated in time from launch (TFL)). 

2. Prediction accuracy. During the early portion of the 
DSS 42 launch pass, the radiometric data, when difiFer- 
enced with the preflight nominal predicts set B09E by the 
Network Operations Control Center (NOCC) pseudoresid- 
ual program yielded the following residuals: 

AHA ~ -.07 degree 

AD2 ~ -29 Hz (S-band) 

AXA - 14 Hz (VCO) 

These can be compared to the 3<r uncertainties supplied 
by the Viking project: 

AHA ~ .002 degree 

AZ>2 - 10 Hz (S-band) 

AXA ~ 55 Hz (VCO; total frequency/trajectory 
unC'ertamty ) 

The residuals, in general, exceeded the 3a uncertainties 
(which must be considered miniscide when compared to 
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previous mission trajectory uncertainties). The hour an^e 
residua] had no impact on the acquisition of the downlink 
and caused no degradation of received signal level even 
thou^ it S&bstantially exceeded the So- uncertainty. 
Additionally, the magnitude of this angle residual is nearly 
equal to t^ hour an^e residual of the Viking 1 laimch 
which places it among the smallest early launch pass angle 
residuals yet achieved. 

The two way doppler residuals shown in Figure 4 
started quite large (approximately -300 Hz), but gradually 
decreased to approximately -29 Hz (or three times the 
expected 3cr value) as the apparent motion of the 
spacecraft approached sidereal rate. 

The difierence between the measured and predicted 
best lock frequencies, AXA, fell well within the total 
trajectory/frequency So* uncertainty. Though larger than 
the AXA for the Viking 1 launch (AXA -0.3 Hz) the 
magnitude of this AXA fell easily within the boundaries of 
the prescribed uplink sweep and caused no problem in the 
uplink acquisition. 

B. One-Way Acquisition at DSS 42 

Acquisition of the Viking 2 one-way downlink at DSS 42 
was reported at 19:27:01 GMT, 25 seconds prior to the 
expected spacecraft rise time of 19:27:26 GMT. This again 
indicates a possible error in the DSS 42 horizon mask used 
in the prediction software. (This problem was first noted 
during the Viking 1 launch. See Reference 1, Section VI- 

B. ) 

The downlink acquisition is illustrated in Figure 5. As 
can be seen, it appears that the receiver was being tuned 
through the region near the predicted one-way doppler 
prior to the expected spacecraft rise time. The signal was 
apparently detected at approximately 19:26:41 GMT with 
Receiver 5 (coimected to the S-band Cassegrain Mono- 
pulse (SCM) Anteima) reported in lock at 19:27:01 GMT. 
Both monitor and tracking data indicate, however, that 
Receiver 6 (connected to the S-Band Acquisition Aid 
(SAA) Antenna) may have sustained lock from as early as 
19:26:41 GMT. 

C. Two-Way Acquisition at DSS 42 

DSS 42 was instructed to perform the following uplink 
acquisition sweep designed according to the specifications 
described in Section U: 

TRANSMITTER ON: 19:30:40 GMT 

START SWEEP: 19:31:00 GMT 

STARTING FREQUENCY: 22.035090 MHz (VCO) 

SWEEP RATE: 180 Hz/min (VCO) 


END SWEEP: 19:32:30 GMT 

ENDING FREQUENCY: 22.035360 MHz (VCO) 

SWEEP DURATION: 90 seconds 

A comparison of the instructed sweep with the actual 
uplink acquisition sweep is depicted in Figure 6. As is 
shown, the sweep began approximately seven seconds 
later than planned. At 19:32:02 GMT, the switch to the 
two-way coherent mode occurred. The two-way downlink 
was very quickly acquired by DSS 42 with the receiver 
back in lock at 19:32:12 GMT. However, as the receiver 
was being locked to the downlink, tuning slowed and 
almost stopped for several seconds. The tuning rate for the 
remainder of the sweep was somewhat slower than during 
the initial portion of the sweep, causing the ramp to take 
fifteen seconds longer than originaUy planned. 

When the doppler extractor was switched from the 
SAA antenna receiver to the SCM antenna receiver, it was 
found that that receiver was in lock on a sideband located 
approximately 10 kHz from the main carrier (see Table 1). 
Receiver lock on the sideband was broken approximately 
two minutes later and the carrier reacquired at 19:34:44 
GMT. 

D. Angle Tracking 

In following the an^e strategy reviewed in Section II, 
the antenna at DSS 42 was initially computer-driven using 
the preflight nominal predicts set B09D, generated at 
laimch minus two hours. The drive mode was changed to 
autotrack at 19:35:51 GMT following completion of the 
uplink acquisition. 

As was the case during the Viking 1 initial pass, it was 
necessary to transfer the uplink to DSS 44 in order to 
allow DSS 42 to switch maser No. 1 into the antenna 
microwave subsystem. In preparation for this reconfigura- 
tion, the drive mode at DSS 42 was changed back to 
computer mode at 20:17:00 GMT. Autotracking was 
resumed five minutes later, following the uplink transfer 
back to DSS 42. After the foilure of maser No. 1 at 
23:11:16 GMT, DSS 42 returned to computer drive for the 
remainder of the pass. 

E. Ranging 

The acquisition of range data at DSS 42 began at 
22:25:02 GMT and continued with generally good results 
throu^ the acquisition of seven range points. Using the 
Pseudo-Differenced Range Versus Integrated Doppler 
(DRVID) technique the range points were evaluated with 
the following results: 
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Acquisitions 

Pseudo-DR 

1/2 

4 

2/3 

-11 

3/4 

376836 

4/5 

' * 

5/6 

-1 

6/7 

352235 


*Unable to compare due to erroneous 
reset of doppler counter. 


As can be seen, two acquisitions, numbers four and seven 
of the seven completed, resulted in erroneous range 
points. Acquisition number four was started at 23:10:02 
GMT and completed two minutes later, during which time 
maser No. 1 failed causing the acquisition to be erroneous. 
The seventh and final acquisition began at 00:15:02 GMT. 
During this acquisition the antenna drove approximately 
two degrees off point causing a bad acquisition. 


The pseudo-DRVID computed for the remaining five 
points indicates that they were of good quality and usable 
for orbit determination. 


(3) Bioshield separation needed to be observed in the 
attitude control system telemetry. 

(4) Two-way doppler and telemetry needed to be 
continu^ without major interruption. 

To resolve these difficulties, it was decided that both 
stations would be configured in the maser by-pas$ mode 
(i.e., the maser would be completely out of the antenna 
microwave system) during the period of high signal level 
and that the S-band Cassegrain Monopulse (SCM) antenna 
would be used for transmission during the roil turn. Thus, 
an uplink transfer would be effected between DSS 42 and 
DSS 44 at approximately launch plus 100 minutes in order 
to provide continuous telemetry and alleviate the need to 
reacquire the uplink while DSS 42 reconfigured the maser 
and transmitter. When it again had the uplink, DSS 42 
would be configured for normal cruise support (i.e., maser 
in, SCM aptenna transmitting) without loss of data to the 
Viking project. 

Following this plan, DSS 44 acquired the one-way 
Viking 2 downlink at 19:29:18 GMT (as was the case with 
DSS 42, this time was considerably earlier than the 
predicted spacecraft rise time). The acquisition is shown 
in Figure 7. Autotracking of the spacecraft began at 
19:33:38 GMT, after completion of the DSS 42 uplink 
acquisition sweep and confirmation of good three-way 
downlink. 


F. Launch Pass Activities at DSS 44 

Since it was designated to serve as a backup station, 
DSS 44 played a somewhat passive role in the Viking 2 
launch phase operations. Its availability however allowed 
some unique (for a launch phase) configurations to be used 
at DSS 42. 

In prelaunch planning sessions, the acquisition planning 
group was presented with the following constraints and 
requirements (among others): 

(1) The downlink signal level at spacecraft rise was 
expected to be strong enough (~ -80 dBm for the 
SCM antenna) to cause saturation of the maser. 

(2) A planned roll turn (at launch plus 150 minutes) of 
the spacecraft to execute a sky map could cause 
unfavorable antenna angles. This would result in the 
possible loss of two-way lock if the S-band Acquisi- 
tion Aid antenna was transmitting. 


The uplink was handed over to DSS 44 by means of a 
tracking synthesizer frequency (TSF) transfer (which does 
not require tuning of the uplink) from 20:18:02 to 20:21:02 
GMT. Telemetry and tracking data continued uninter- 
rupted during this period. 

DSS 44 was reconfigured to normal cruise configuration 
at 20:27:40 GMT and continued tracking in the three-way 
mode until spacecraft set at 00:28:20 GMT. 

IV. Summary 

The acquisition strategy used diving the Viking 1 
launch was again successfully implemented during the 
Viking 2 launch. This strategy had been carefully plarmed 
to accommodate even the most unfavorable launch 
possibilities. The culmination of this careful plaiming and 
extensive training was the successful acquisition of Viking 
2 and completion of the Viking launch phase. 
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Table 1. Doppler residuals 


GMT 

Residuals 

Conunents 

19:30:30 

-275.331 

Final good one-way residual 

19:31:00 

-15029.743 

Tuning one-way doppler flagged 
two-way 

19:31:20 

-12163.310 

Tuning 

19:31:40 

-7095.466 

Tuning 

19:32:02 

-1650.137 

Receiver out of lock — switch to 
coherent mode 

19:32:12 

-298.039 

Good two-way residual 

19:32:20 

-291.406 

Good two-way residual 

19:32:30 

-286.623 

Good two-way residual 

19:32:34 

-10026.788 

Switch doppler extractor — SCM 
RCVR on sideband 

19:33:00 

-10019.764 

Receiver on sideband 

19:33:30 

-10007.272 

Receiver on sideband 

19:34:00 

-9993.609 

Receiver on sideband 

19:34:36 

-11827.154 

Receiver out of lock 

19:34:44 

-228.606 

Receiver in lock on carrier 

19:35:00 

-203.622 

Good two-way residual 


ORIGINAL PAGE IS 
OP POOR QUALIIY 
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XA, MHi (22.03XXXX) 



Fig. 1. Best lock frequency at DSS 42, Viking 2 launch, September 9, 1975 
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Fig. 2. DSS 42 Viking 2 launch, September 9, 1975 
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20 : 00:00 


21:00:00 22:00:W 23:00:00 00:00:00 

GMT, 9 SEPTEMBER 1975 

Fig. 3. Elevation at OSS 42, Viking 2 launch, September 9, 1975 
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20 : 00:00 

GMT, 9 SEPTEMBER 1975 

Fig. 4. Doppler deviation from preflight nominal predicts, DSS 42, predicts set B09E 
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GMT, 9 SEPTEMBER 1975 


Fig. 5. Initial downlink acquisition at DSS 42, Viking 2 launch 

















Network Loading Visibility for 
Management 

D. M. Enari and C. A. Holritz 

Network Operations Office 


The responsibilities of the Deep Space Network (DSN) Operations Scheduling 
Group (NOSG) are outlined. The long-range scheduling portion of the NOSG is 
explained in detail. Examples of the DSN Forecast, published quarterly by the 
NOSG, are also included. 


The Deep Space Network (DSN) Operations Schedul- 
ing Croup (NOSG) has the responsibility of forecasting, 
allocating and scheduling the utilization of the DSN 
facilities in response to requirements from all of its users 
such as flight projects, radio science, and DSN develop- 
ment. In addition, requirements which are internal to the 
DSN such as maintenance, engineering, training and 
facility updates must be forecast, allocated and scheduled. 
The NOSG is also responsible for providing management 
with information which will give them visibility concern- 
ing mission support, facility utilization and Network 
loading. 

The ofiicial NOSG publications which comply with 
these responsibilities are as follows: 

1. Seven-Day Schedule — Detailed schedule for the 
coming week. Published weekly. 

2. Mid-Range Schedule — Time-sensitive schedule of 
facility activities for each day for the next 8 weeks. 


Provides the basis for the Seven-Day Schedule. 
Ihiblished monthly. 

3. Forecast — Gross planning of Network allocation for 
the next 3 years. Provides the basis for the Mid- 
Range Schedule. Published quarterly. 

To provide management with the assessment required, a 
method for displaying the data is available. It consists of a 
set of Network Requirements and Allocation Charts, 
which display the requirements placed on the Network by 
users, the capability of the Network and the platmed 
allocation of that capability to the various users in 
accordance with established priority guidelines. The total 
requirements, allocations and capability are displayed at 
the bottom of the chart just above the calendar. Major 
milestones are displayed across the top of the chart as an 
aid to correlating requirements (loading) with events. 
Figure I is an example of this chart. Note that it covers 
two fiscal years, 1975 and FY 1976, in 4-week 
increments, for the 26- and 64-meter stations at the 
Goldstone longitude (DSSs 11, 12 and 14). A complete set 
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of charts comprises 9 charts; 3 for each of the three 
longitudes at which the DSN has Deep Space Stations 
(DSSs), namely the Goldstone longitude (DSSs 11, 12 and 
14), the Australian longitude (DSSs 42, 43 and 44) and the 
Spanish longitude (DSSs 61, 62 and 63). The three charts 
for each longitude cover the years of interest; one prior 
year, the current year and two future years. 

For the purpose of long-range forecasting, the following 
guidelines are used in establishing the capabilities of the 
Network. The available station hours, as defined by the 
budget, are divided into two basic categories. Seventy 
percent are allocated for DSN user support and thirty are 
allocated for facility internal activities. Also, for long- 
range forecasting, a “pass” is considered to be an average 
of 12 hours in duration. To determine the capability in 
“passes per 4 weeks,” the station hours are multiplied by 
0.7, divided by 12, rounded to the nearest whole number 
and then multiplied by 4. Thus, a station which is stafiFed 
for 120 hours per week (3 shifts) would be allocated . 84 
hours per week for user support, which will provide a 
capability of 7 passes per week or 28 passes per 4-week 
period. 

For “quick-look” display and for convenience in 
plaiming support of future requirements on the Network, 
a set of “residual” graphs is made. Figure 2 is an example 
of this chart. These graphs show the difference between 
the total requirements for each anteima size and the total 
capability at each longitude. The time scale shown in Fig. 


1 is used. When the requirements exceed the capability, 
the difference is plotted as a positive ( + ) value and 
indicates an “overload” and a conflict. If the requirement 
is less than the capability, the difference is plotted as a 
negative (-) value, representing an “imderload” and no ! 

conflict. The zero line represents the capability. The scale 
is in increments of 8 passes per 4-week period. This value 
was chosen for convenience because, to a first-order 
approximation, 8 passes represents the difference in 
capability resulting from changing the staffing at a station 
by one shift over a 4-week period. 

Because of the method of display and the scales that 
have been chosen, these charts and graphs provide [ 

management with the visibility of Network loading and 
provide the basis on which decisions can be made 
concerning budgeting, station staffing and negotiations 
among Network users. 

The DSN has some users whose requirements are of a 
nature which makes it possible to use the charts to 
“smooth” the loading on the Network. This can be done 
because their requirements, unlike those of flight projects, , 

are not sensitive to time of day or specific days. Radio 
science, DSN development. Pioneer 6-9 coverage and 
certain types of tests fall into this category. The detailed 
scheduling of this support is usually accomplished during 
the eight-week and seven-day scheduling process in a 
maimer which meets the broad requirements and also aids 
in normalizing the loading on the Network. 
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DSN Research and Technology Support 

E. B. Jackson 

R. F. Systems Development Section 


The activities at the Venus Station (DSS 13) and the Microwave Test Facility 
(MTF), both operated by the Development Support Group, during the period Oct. 
16-Dec. 7, 1975 are discussed and progress noted. Continuing testing of the 
remote-controlled automated station is noted as well as routine pulsar 
observations.. Automatic stability-reliability testing of the station maser- 
receiver-noise adding radiometer combination is described along with the data 
collected while so doing. Comparative measurements on the production version 
of the dichroic plate installed on the 64-m antennas are described along with 
mention of testing of a National Bureau of Standards (NBS) radiometer. 
Engineering measurements on the microwave power transmission test setup are 
described and a progress report of the X-Band Radar, Transmit-Receive; K-Band 
Receive (XKR) feedcone rehabilitation is given. Special testing of the new design 
feedcone for the Unified S-Band stations is mentioned, and measurements on 
interference received from the navigational equipment aboard military or 
commercial aircraft are described. Reporting of routine maintenance and 
support of clock synchronization transmissions and various radio science 
experiments is also included. 


During the period Oct. 16-Dec. 7, 1975, the Develop- 
ment Support Croup, in operating the Venus Station (DSS 
13) and the Microwave Test Facility (MTF) supported 
Various programs as discussed below. 

I. Station Automation 

Demonstration of a remotely operated automated 
station is planned using DSS 13 as the test station. 
Although a successful demonstration has been conducted 


(DSN Progress Report 42-30, pp. 214-221), refinement of 
software and reliability testing is continuing. 

The antenna waveguide switch controller was modified, 
and computer monitoring of the position of the polariza- 
tion selector and ambient load selector switches was 
provided. System testing and automated pulsar tracking 
continued. A total of 45 hours of station support time, 
including 17 hours of automated pulsar tracking, was 
provided during this period. 
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II. Pulsar Observations 

In support of the Radio Science Elxperiment “Pulsar 
Rotation Constancy,” DSS 13 provided 56-3/4 hours of 
observations during which the emissions from the pulsars 
tabulated in Table 1 were recorded. These data, recorded 
at 2388 MHz, left-circular polarization (LCP), are used to 
determine precise pul^-to-pulse spacing, pulse shape, and 
pulse power content of the signals emitted by these 
pulsars. 

III. Maser-Receiver-Noise Adding Radiometer 
(NAR) Reliability-Stability Testing 

Stability and reliability testing of the DSS 13 receiving 
system is conducted automatically during non-operationai 
station periods. The 26-m antenna is pre-positioned in 
elevation and azimuth, and the NAR automatically records 
total system temperature as a function of time. The 
rotation of Earth sweeps the fixed antenna beam across 
the sky, resulting in generation of a radio brightness 
temperature sky map in addition to the data on stability 
and reliability of the system. During this period 391-1/4 
hours of such data were automatically recorded with the 
antenna at 360 degrees azimuth and progressively 
positioned from 53.3 degrees to 52.7 degrees in elevation. 
Testing is done at 2295 MHz, using right-circular 
polarization (RCP) on the 26-m antenna. 

IV. Dichroic Plate Measurements 

The production versions of the dichroic plates used on 
the 64-m antennas contribute more system temperature 
than design and prototype testing indicate they should. 
Taking advantage of the removal of the DSS 14 Dichroic 
Plate for refurbishment of the XKR feedcone, DSS 13 
positioned this plate atop the X-Band Low-Noise Antenna 
Measurement (XLA) feedcone which is operational for on- 
the-ground testing at DSS 13. A spot frequency measure- 
ment resulted in a system temperature of 20 K vice 16.5 K 
without the plate. Testing, using an angled flat sheet of 
aluminum and an absorber batten, indicated that a 
reflected signal entering from the side was contributing 
part of this additional system temperature. 

Extensive additional measurements were conducted 
comparing the sensitivity of the plate to the operating 
frequency. The paint was then removed from the plate 
and “after” measurements conducted. Some frequency 
sensitivity was observed, but more testing of the 
prototype dichroic plate is necessary before final conclu- 
sions can be drawn. The plate was repainted and returned 
to DSS 14 for reinstallation. 
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V. National Bureau of Standards (NBS) 
Radiometer Testing 

In support of Antenna Microwave Advanced Engineer- 
ing and X and K-Band Propagation Calibration, an NBS 
design radiometer was tested at DSS 13. This testing, 
conducted at 2278.5 MHz, was to evaluate the stability, 
resolution, and overall usability of the NBS-designed 
radiometer. DSS 13 provided 25 hours of station support 
of which 22 hours were actual observing using the new 
radiometer. The dynamic range of the radiometer was less 
than desired, and further testing is planned after redesign 
and rebuilding. 

VI. Microwave Power Transmission 

Continuing with system testing (DSN Progress Report. 
42-30, pp. 214-221), measurements have been made of the 
sensitivity of the rectenna efiiciency to angle of incidence 
of the arriving rays. Additionally, testing of system 
performance as a function of operating frequency and 
output power has been performed. By moving the 26-m 
subreflector, data have also been taken of the recovered 
power level as a function of transmitting antenna focusing. 
This last series of measurements indicates that additional 
efficiency could be gained by focusing the antenna on the 
rectenna instead of at infinity as was the case during 
preceding system testing. 

System performance was demonstrated to a number of 
visiting groups including a contract motion picture film 
crew from the United States Information Agency. DSS 13 
has provided 15-3/4 hours of system support, with the 
transmitter operating at various powers up to 250 kW 
during this period. 

VII. X-Band Radar 

In preparation for an ’ additional series of radar 
observations, the XKR feedcone was removed from DSS 
14 and brought to DSS 13 for extensive refurbishment of 
the radar system, including feedhorn, waveguide, 
waveguide switches, buffer amplifier, and protective 
circuits. 

All waveguides, including the feedhorn, were removed. 
The voltage standing wave ratio (VSWR) of each piece of 
guide, as well as groups of pieces, was measured. All 
guides were then chemically cleaned, flanges lapped, and 
VSWR measurements again made. “Runs” of waveguide, 
as installed in the cone, were cormected and VSWR 
measurements again made. The waveguide runs were 
“tuned” for minimum VSWR by careful deformation of 
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the waveguide walls and the system was reinstalled, along 
with the cleaned feedhom. 

The buffer amplifier was completely bench-tested, and 
the prototype crowbar circuit (DSN Progress Report 42- 
30, pp. 214-221) was replaced by the production version. 
Response time measurements were made on all of the 
protective circuits (arc detectors, reflected power detec- 
tors, etc.), and an inhibit circuit was added to the 
traveling-wave tube (TWT) power supplies to prevent the 
power supplies from automatically restarting when 
shortcircuited by the crowbar circuit. 

At the request of the Viking project, arrangements were 
made with Varian for accelerated delivery of a repaired 
klystron which was installed into the XKR feedcone on the 
DSS 14 antenna on Dec. 4, and system testing continued 
with two klystrons. This testing is still in progress. 


VIII. Unified S-Band Feedcone Testing 

Negotiations had previously been completed with 
Goddard Space Flight Center (CSFC) for dual-carrier 
testing at the Microwave Test Facility (MTF) of the new 
typ>e feedcone to be installed into Spaceflight Tracking 
and Data Network (STDN) stations. 

Using a test plan prepared by the feedcone contractor, 
dual-carrier testing at total power levels up to 24 kW was 
performed. The MTF provided dual transmitters of 
variable frequency and power as well as the necessary 
exciters, frequency sources, and a maser-receiving system 
with which to detect the intermodulation products 
generated by dual-carrier operation. 

Testing was conducted at various powers up to 24 kW 
and transmitter operating frequencies from 2080 to 2120 
MHz. System performance was good, with intermodula- 
tion products weaker than -141 dBm. Testing at MTF has 
been completed with further testing to take place at 
Goldstone STDN. 


IX. DSN Receiver Interference Susceptibility 
Measurements 

Among the possible sources of interference to flight 
operations is the radiation from commercial or military 
aircraft navigational gear, in particular the second 
harmonic from the Distance Measuring Equipment 
(DME). Several occurrences of interference seem likely to 
be from this source. 


Tests were performed at DSS 13 utilizing the NASA 
shuttle aircraft to fly patterns around DSS 13 and operate 
the aircraft DME on each of the possible channels. When 
the aircraft was in the 26-m antenna beam, at ranges of 
approximately 2 km, interference was observed and 
recorded from each channel of the DME. Signal levels of 
approximately -138 dBm for the interference were 
recorded. Having the aircraft fly toward the antenna, 
remaining in the beam, was particularly useful as it 
allowed relatively long-term observation of the interfering 
signal on both a spectrum analyzer and a chart recorder. 


X. Antenna Maintenance 

The development activities which take place at DSS 13, 
particularly the development of various computer- 
controlled automation schemes, place unusual stress on the 
anteimas. In particular, the azimuth drive gear boxes on 
the 26-m anteima have failed in the past In order to 
forestall, if possible, future catastrophic failure, these gear 
boxes are being modified and refurbished. The last of the 
four gear boxes was removed from the 26-m antenna, 
replaced with a reworked spare, and sent to the vendor 
for modification and refurbishment. 

Scheduled periodic checking revealed that the 9-m 
antenna was unbalanced. Removal of 907 kg (2000 Ib) of 
balancing weight from the elevation axis brought the 
system into balance for smoother tracking in the elevation 
axis. 


XI. Planetary Radio Astronomy 

In support of the Planetary Radio Astronomy experi- 
ment, DSS 13 measures radiation received at 2295 MHz 
from the planet Jupiter and various radio calibration 
sources. These measurements utilize the 26-m antenna, the 
maser-receiving syrtem, and the Noise Adding Radiometer 
(NAR). During this period, observations were made of the 
calibration sources tabulated in Table 2 in addition to 
measurements of Jupiter itself. A total of 40 hours of 
observations were made. 


XII. Platform Parameters, Very Long Baseline 
Interferometry (VLBI) 

In support of this development project, DSS 13, in 
conjunction with DSS 43, provided 16 hours of VLBI 
observations. During these 16 hours, 97 sources were 
observed, and the received data recorded onto magnetic 
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tape utilizing a modified TV video recorder. These 
observations were made at 2290 MHz. 

XIII. Clock Synchronization Transmissions 

Although some troubles have been encountered with 
this system, three transmissions were made to DSS 42-43 


and three to DSS 61-63. DifBculties encountered include 
failure of the programmed oscillator due to a power 
supply failure and a changed capacitor value in the 
computer interface. Additionally, the circulating water 
pump began leaking so badly that replacement was 
necessary. Repairs have been accomplished and the 
system is now operational. 
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Table 1. Pulsars observed at OSS 13 


0355 +54 
0525 +21 
0736-40 
0823 +26 

0833 -45 
1133 + 16 
1237 +25 
1604-00 

1642-03 1911-04 

1706-16 1929+10 

1749-28 1933 + 16 

1818 -04 2045-16 

2121+51 

2218+47 

Table 2. 

Radio calibration sources observed at OSS 13 

3C17 

3C273 

3C309.1 

NGC 7027 

3C48 

3C279 

3C348 

Virgo A 

3C123 

3C286 

3C353 
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N76-18181 


An Experiment in Dynamic Modeling for a 
Complete Solar-Powered Energy System 

C. L. Hamilton 
TDA Planning Office 


Completion of a prototype dynamic model simulating the performance of a 
solar-powered energy system is described. A set of hypothetical components is 
specified, and the outcome of test analyses of the resulting system is outlined. 
On the basis of this exercise, it appears that the dynamic modeling technique 
will constitute a useful and convenient tool for analyzing performance of time- 
dependent systems. 


I. Introduction 

A previous report (Ref. 1) has described in general ' 
terms a technique under development to model the time- 
dependent performance of solar-powered energy systems. 
In that article, a hypothetical system intended to provide 
all the energy needed by the Coldstone Space Communi- 
cations Complex was outlined. Two of the three computer 
programs needed to embody the model of that system 
were discussed. The third program has been finished, and 
some experimental runs have been made to examine the 
modeled performance of the entire system. 

The Coldstone Energy Conservation Project is no longer 
directed toward development of a full energy system for 
the Complex. For that reason the whole-system model has 
not been refined past the prototype stage, and the 
representative programs (SUN, WIND, and SENSMOD2) 
are implemented without the range of user options that 
would be included in a fully operational - version. 


Construction and manipulation of the prototype have 
constituted a useful and informative experiment that has 
contributed to the definition of systematic processes for 
creating dynamic system models. This report will 
summarize the completion of the prototype model and 
results of the exploratory performance analyses that were 
carried out with it. 

II. Characterization of the Test System 

In Ref. 1 the components of a hypothetical baseline 
energy system were identified, and information flow 
between the computational modules representing them 
was outlined. These are illustrated in Figs. 1 and 2. Before 
the model can be used to carry out a performance 
analysis, key characteristics of each component must be 
specified. A component’s key characteristics are those 
parameters which are needed in the corresponding 
computational module to calculate the outgoing informa- 
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tion from the incoming information. In general, character- 
ization can be carried on at several levels of detail, 
depending on information available and the goals of the 
analysis. A heat engine, for example, can be specified 
simply in terms of size and overall efficiency for purposes 
of, say, a quick-scan evaluation. If more precise system 
analysis is desired, the specification can be made more 
detailed, including variation of efficiency with working 
fluid and rejection temperatures. Even a complete ste(>- 
by-step calculation of cycle performance, treating the 
engine as a subsystem in itself, may be incorporated if 
desired. 

The following paragraphs contain a listing of key 
characteristics for the hypothetical system components 
put together for the purpose of exercising the dynamic 
model. Note that these are gross component characteris- 
tics derived from a first guess at some performance 
specifications that would be required to make a whole 
energy system that is economically feasible. They do not 
necessarily correspond to physical characteristics of actual 
existing or emerging components. Analysis at this level 
corresponds to examination of a system's operation in 
terms of subsystem design goals. 

A. Solar Subsystem 

Program SUN, described in Ref. 1, embodies the model 
of a subsystem consisting of solar collectors, heat storage, 
and a heat engine for generation of electricity. 

For purposes of exploratory analysis, 200,000 m^ of fixed 
flat collectors made mostly of glass were postulated. They 
were assumed to face south, tipped up at 35.4 degrees 
from horizontal. Heat loss rate was taken to be about 350 
watts per square meter at a collector temperature 300°C 
above ambient; loss was assumed to be predominantly 
radiative. An organic fluid was specified for heat removal. 
The fluid flow rate was controlled to maintain the outlet 
temperature between 300°C and 340°C (flow was stopped 
if the fluid temperature was below 300°C). Acceptsmce of 
incident radiation was assumed to decrease drastically at 
angles of incidence greater than about 50 degrees from 
normal. 

Heat storage was accomplished with a set of narrow 
tanks holding a total of 3.8 X 10^ kg of the organic fluid. 
The tanks were assumed to be well enough insulated so 
that heat loss from them was negligible over the time span 
of a few days. Fluid in them was assumed stratified into 
two well-defined zones. Immersed in each tank was a heat 
transfer tube through which the heat engine working fluid 
was circulated. Flow to the engine was maintained at a 
constant rate as long as enough heat remained in the high- 


temperature zone to last through the next integration 
interval. When there was insufficient heat stored, at the 
high temperature (roughly 300°C to 340°C), fluid flow to 
the engine was stopped. The total storage unit was sized 
to hold one day’s output from the collectors under 
conditions of maximum insolation, estimated on the basis 
of annual average figures. 

A Rankine cycle engine was included to generate the 
DC output from the solar subsystem. It was to use the 
same organic liquid as its working fluid. It was postulated 
to run at 90% of Carnot efficiency with a mechanical 
efficiency of 80% and generating efficiency of 90%. Total 
fluid flow to the engine was set at 14,200 kg/h, estimated 
to use up in a 24-hour period the whole output of the 
collectors on a day with maximum insolation. 

B. Wind Subsystem 

WIND is the program used to simulate the perform- 
ance of a wind-driven generator. Wind turbine output was 
taken as proportional to the cube of wind speed whenever 
that fell between 8 and 32 km/h. No generation occurred 
when wind speed was outside that range. Maximum 
output rate (at 32 km/h) was 4000 kW. 

C. Conversion and Storage Subsystem, SENSM0D2 

Using the outputs from SUN and WIND, SENSMOD2 
models the remainder of the energy system, including 
hydrogen generation and storage, direct DC to AC 
conversion, dual-fuel engine generators, load consisting of 
waste heat utilization and electrical demand, and the 
chosen dispatch strategy. Characterization of the compo- 
nents represented by SENSMOD2 was done in very basic 
terms. Elfficiency of the electrolysis unit was set at 80%. A 
maximum capacity for hydrogen storage was not defined; 
the amount of hydrogen in storage was monitored, and the 
gas was defined to be available for use long as enough 
remained to satisfy the next hour’s demand. The engine- 
generator was specified to be 33% efficient using either 
hydrogen or diesel fuel. DC to AC conversion was 
assumed . to proceed at 85% efficiency using a converter 
with 1250-kW capacity. Diesel fuel in storage was 
monitored, and whenever the amount there fell below 
2,000,000 kWh a “load” of 200,000 kWh was added to it 
A base waste beat load of 600 kW that varied ±10% with 
ambient temperature was postulated. Electrical load was 
approximated very simply, using a cumulative density 
function that was linear between a base defined as half the 
average load and a peak defined as 1.5 times average load. 

For the prototype Dispatch Module the following 
strategy was devised. All power from the wind turbine was 
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used for electrolysis. Power from the solar subsystem was 
routed through the DC to AC converter to the extent that 
it was required and available; any solar power not used 
directly was consumed by the electrolysis unit. The 
engine^enerators were required to operate to meet waste 
heat demand, producing an equivalent base amount of 
electricity. Electrical demand above that base was met 
with direct power from the converter to the extent it was 
available, up to converter capacity. Demand exceeding 
that level was to be satished by additional engine- 
generator operation. Hydrogen would be used to fuel the 
generators as long as there was enough; then the engines 
would switch to operation on diesel fuel. 

SENSMOD2 is a program analogous in design to SUN 
and WIND. None of the components of the conversion- 
storage subsystem required modeling with variable step 
size integration, and the simulation could be run with 
update at hourly intervals. While SENSMOD2 was being 
put together, a requirement for ordering computational 
modules in the program was discovered. The relationships 
simulating the hydrOgen-oxygen subsystem, the engine- 
generator, and fuel storage all occur in the DOWN- 
STREAM subprogram, and each of those computational 
modules requires information from one of the other two. 
Thus the Hydrogen-Oxygen, Engine-Generator, and Fuel 
Storage Modules had to appear in the appropriate 
sequence in the program. TTiis ordering requirement 
places some restriction on flexibility in manipulating the 
system model as embodied in the program, as it represents 
one more piece of information not related to the system 
imder analysis that the user has to remember. It may be 
that more skillful design of the Dispatch Module will 
avoid this drawback in future eflForts. 

The computational module representing the dispatch 
function differs from those described in Ref. 1. A “normal” 
computational module contains relationships defining 
component performance. Those relationships with no 
need for information from other computational modules 
(components) are contained in the UPSTREAM subpro- 
gram, while the equations using parameters transferred 
from another module appear in subprogram DOWN- 
STREAM. Between the two, CROSS-COUPLE carries out 
the transfer of variables from module to module. While 
the Dispatch computational module is made up of a 
relatively complex set of logical operations, its real 
function remains one of transferring information. The 
pieces of information have been manipulated before 
transfer, but not used to calculate new variables. Thus the 
Dispatch Module belongs strictly in the CROSS-COUPLE 
subprogram. 


III. Performance of the Test System 

To test the performance of the model and programs, 
operation of the hypothetical system summarized above 
was simulated over the course of one month, January. For 
input to SUN, a table of hourly insolation values was 
generated using the ASHRAE model (Ref. 2) with 
clearness number 1.05; each day’s values were modified by 
a probabilistically determined factor to include the effect 
of unclear days. Ambient temperature inputs were from a 
file containing a composite of experimental hourly values 
measured at Coldstone over about three years. Wind 
speed inputs to WIND were drawn from a cumulative 
distribution based on three years’ measurements at 
Coldstone. Output from SUN and WIND was stored as 
hourly average values for input into SENSMOD2, which 
also used the ambient temperature file. 

A run of SENSMOD2 was made, simulating the system 
as specified. Several measures of system jjerformance were 
output, among them the amount of hydrogen in storage at 
the end of each day and daily totals for solar-generated 
electricity, wind-generated electricity, hydrogen con- 
sumed, and diesel fuel consumed. Production of electricity 
from wind turned out to be a small fraction of the output 
from the test system (less than 10%). Figure 3 summarizes 
the salient features of system operation for the month of 
January. ’The top line represents daily totals for heat 
gathered by the solar collectors; it was generated by 
program SUN. Solar subsystem daily total production of 
electricity varied between 58,500 and 198,000 kWh. 
Hydrogen in storage, starting from an initial value of 
200,000 kWh, decreased rapidly while consumption of 
diesel fuel (graphed as a cumulative value) rose to a total 
of over 640,000 kWh by month’s end. 

Inclusion of waste heat utilization in the system 
required that the engine-generators nm ail the time to 
meet that demand, consuming a great deal of hydrogen in 
the process. Direct use of the solar-generated electricity 
would seem to offer a significantly more efficient strategy. 
As a first parameter . variation, then, the provision for 
waste heat utilization was removed, electrical demand was 
increased to compensate, and the capacity of the DC to 
AC converter was doubled to 2500 kW. Figure 4 shows 
the effect of these changes on cumulative diesel demand 
and hydrogen storage level (heat collected and solar 
subsystem output remain unchanged). Contrary to expec- 
tations, lifting the requirement for compulsary base 
operation of the generators did not lower diesel consump- 
tion (which is the measure of overall system performance 
chosen here) but rather increased it substantially. 
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Examination of the hourly electrical output from the 
solar subsystem revealed that the heat engine there 
exhausted the heat in storage in an average of 14 hours, 
and there was no solar power production the rest of the 
day. In an efiFort to better match power generation with 
the load distribution, which was spread over all hours, the 
output rate of the heat engine in the solar subsystem was 
decreased by halving the flow rate of working fluid to the 
engine (from 14,200 to 7100 kg/h). SUN was re-run with 
that modiflcation, and the resulting data input for a third 
run of SENSMOD2. Hydrogen accumulation and fuel 
demand observed under these conditions are plotted in 
Pig. 5. This last modification significantly reduced diesel 
fuel consumption. 


IV. Summary 

A dynamic model simulating the performance of a 
complete energy system has been formulated and 
embodied in a set of computer programs. Those programs 


are functionally complete but do not contain all the 
features providing user convenience that would be present 
in fully implemented versions. Characteristic parameters 
were specified for a set of hypothetical components, and 
the programs were used to carry out some exploratory 
analyses of the resulting system. Parameter variations 
reflecting system and component modifications were easily 
and quickly made. The test runs described provided 
considerable insight into the behavior of a fairly complex 
solar-powered system and allowed rapid refinement of the 
system leading to improved performance (by one measure 
at least). 

Additional work in dynamic modeling involves a 
program that produces instantaneous and integrated 
performance measures for solar collectors alone, using 
experimentally measured solar radiation inputs at short 
time intervals. Collector characterization packages are 
undergoing continual refinement. Also in progress is a 
detailed analysis of a solar heating and cooling system to 
be designed using commercially available components. 
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Fig. 2. Computational modules and data flow for hypothetical system 
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Fig. 5. Performance of hypothetical system without waste 
heat utilization after heat engine rating cut in half 
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